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ABSTRACT 


DEVELOPMENT OF A SIMPLE, SELF-CONTAINED 
FLIGHT TEST DATA ACQUISITION SYSTEM 


I This report describes work done under a continuing program to 

i develop a simple, self-contained flight test data acquialtlon system. 

' In the past. Instrumenting an r.lrplane for flight testing has taken 

a great deal of time and money. With recent advances In sensor and 
microprocessor technology, a simple, low-cost system could be developed 
which would be applicable to general aviation airplanes. 

This system was conceived to obtain performance and stability 
characteristics of airplanes. The design criteria for the system 
were that It be ersy to Install, self-contained, and simple; that 
It require no speclal/difflcult flight techniques; and that It be 
applicable to general aviation airplanes and low In cost. 

The system developed meets these criteria for doing lontltudinal 
and lateral stability analysis. The package consists of three modules. 
These are 1) microprocessor controller and data acquisition module, 

2) transducer module, and 3) power supply module. The system Is easy 
to Install and occupies space In the cabin or baggage compartment of 
the airplane. All transducers are contained In these modules except 
the total pressure tube, static pressure air temperature transducer, 
and control position transducers. 

The data reduction technique used was the NASA-developed MMLE 
program. This has been placed on a microcomputer, and all data reduc- 
tion Is done on the microcomputer. This greatly reduces the cost of 
the data reduction. Also, when compared with the analogue recording 
techniques, still being used, there has been a large Improvement in 
the accuracy of results. 

The flight testing program undertaken has proven both the flight 
testing hardware and the data reduction method to be applicable to 
the current field of general aviation airplanes. 

This report describes the instrumentation system developed, 
the data reduction method used, and Important results of the flight 
test program. 
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1. INTRODUCTION 


1 


This rsport dsserlbss work conplstsd during the second phsse 
of s continuing progrra sponsored by the NASA Dryden Flight Resesrch 
Center.* This progrsm wss scconpllshed during the period Jsnuery 
21, 1979, through Februsry 15, 1981. The progrsm encompasses the 
development of s simple, self-contslned flight test date acquisition 
system. To date the program has consisted of two phases: 

PHASE I 

• A literature survey of flight testing methods (presented 
In Reference 1). 

• The development and testing of a proof-K>f -concept system 
capable of longitudinal stability analysis (presented In 
Reference 2). 

PHASE II 

• Development and testing of a system capable of longi- 
tudinal and lateral stability analysis. 

This report describes in detail the system concepts selected, 
as wnll as results of the flight test program used to show the 
validity of these concepts, as of the completion of Phase II. 

The purpose of this project, and the design criteria developed 
are contained In Chapter 2. The literature survey (Reference 1 ) 
has been used as a primary data base for establishment of these 

*Funding provided under NASA Grant NSG 4019 (FRL/CRINC 4070). 
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criteria. Other Inpute have com from talka with peraonual in the 
general aviation induatry and of HASA Dryden Flight Raaaarch Center. 

Chapter 3 deaeribea the hardware aeleeted and Mnufactured to 
Met the design requireMnta. 

The instruMntation package eaploye tranaducere to allow 
both longitudinal and lateral etability analyeia of general aviation 
type airplanee, although it can eaaily encompass most other typea 
of airplanes. Due to the nature of the data reduction Mthod 
utilized » a minimum number of high-accuracy tranaducere are required. 
Data from the tranaducera are recorded using an on-board micro- 
procesaor and digital cassette recorder. This has proven a simple » 
reliable method to obtain accurate flight data. 

The system has been designed to allow it to be placed in the 
aircraft with a minimum amount of aircraft modification. A recharge- 
able battery pack was selected for airborne power to reduce the 
number of airplane modifications required. This has allowed total 
isolation from the aircraft electrical syatemi’', which simplifies 
installation, enhances safety, and eliminates many electrical noise 
problems in the transducer signals. Tbii;; transducers are all con- 
tained in one module, except for the following: 

- total pressure probe, 

- temperature probe, 

- static pressure probe, and 

- control position transducers. 

A minimum of installation is also required for these devices, as 
they are literally "sticky- taped” to the airframe. 
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PrtMnttd In Chnptnr 4 la an avaluatlon procadura uaad to 
•alaet a ground-basad data raductlon conputar. Dua to tha axtanalva 
mathaaatical procadura uaad for data analyala, a po%Mrful, hlgh-Iaval 
languaga adcrocoaputar la raqulrad. Tha avaluatlon nathod uaad la 
daacrlbad harat aa wall aa tha eomputar aalactad and uaad for thla 
prograa. 

Tha total flight taatlng procaaa la Includad aa Chaptar 5. 

Dlacuaaad ara tha varloua computar programa and oparatlng tachnlquaa 
davalopad. Tha haart of thla ayatan la tha Modlflad MaxlBum Llkallhood 
Eatlsatlon (MMLE) mathod which haa baac uaad for data raductlon. Tha 
natheraat^ca of thla tachnlqua are Includad aa Section 5.6. 

The flight teat program used for system development is included 
In Chapter 6. Tests have been performed using the ICU-FRL* Cessna 172 
airplane (shown In Figure 1.1). The type of flight test maneuver 
required Is discussed, and .results of the actual flight testing are 
presented. 

A flight test program was conducted at Cessna Aircraft to evaluate 
the spin properties of their model 172 airplane. The data management 
portion of the KU-FRL system was used In conjunction with Cessna-supplied 
transducers for data acquisition and analysis. This program Is described 
In Chapter 7. 

Conclusions to be drawn as a result of the work carried out 
under this program, and recommendations for further work are Included 
In Chap’jers 8 and 9. 

’ft 

KU-FRL “ University of Kansas Flight Research Laboratory. 


3 





Figure 1.1 Experimental configuration of the test airplane. 


References, and reports describing this project are presented 
in Chapter 10. 

Appendix A includes descriptions of ail programs required for 
system operation. 

There appears to be some confusion over the many reference axes 
systems used in airplane analysis. Included In the list of symbols 
is Figure v. 1, which explicitly defines the axes system used in this 
report. Appendix B is included to allow conversion of results in 
this report to several other standard axes definitions. 
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hhlrkitltitli'klt 


The system constructed under this program has proven under 
flight test the validity of the concept selected for longitudinal 
and lateral stability analysis. Throughout this flight test program, 
using the Cessna 172 airplane, the acquisition package performed 
reliably. 
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preceding page blank not filmed 

2. PURPOSE OF PROJECT 

Flight testing has always required a high degree of complex 
Instrumentation to get accurate results. This, In the past, and 
still evident today, has taken a great deal of time and money to 
equip each Individual flight test article. Traditional systems 
are placed on aircraft on an Individual basis, utilizing what Is 
available at that time, coupled with the specific requirements of 
a particular program. This has never really led to ideal or 

totally thought-out systems, and normally results In high costs 
or In too much time being required for instrumenting the airplane. 

With the accurate Instrumentation available today, and with the 
recent advances In microcomputer technology, It was seen that an 
accurate, multipurpose data acquisition system could be developed. 

The system described here has been developed to do Just that« 

The basis for design of this system is as laid out here. 

EASE OF INSTALLATION - This has been a major design consideration. 
If possible, NO permanent modification should be done to the airplane. 
The system must be universally easy to install and should require a 
minimum of ir^^tallatlon time and no special procedures. This factor 
includes calibration of the system installed on the airplane. 

SELF-CONTAINED - The system should be totally self contained. 

This should include all data sources, data recording methods, power 
requirements and data reduction techniques. 
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SIMPLE - Th« aystam must be slrople in concept end eeay to uae. 
The need for complex Inatruntentetion, difficult celibretion, end 
epecielized operator knowledge muet be kept to a mimlmun. 

FLIGHT TESTING > The system should not require any specialized 
piloting techniques to obtain accurate results. 

CLASS OF AIRCRAFT - The system to be developed is primarily 
applicable to the general aviation type airplane. This criterion 
does not restrict the methods and thecrles, but it does define the 
requirements for the transducer ranges and accuracies. 

RESULTS - The system Is aimed at stability and performance 
parameter identification, but it must permit adaptation to other 
test requirements. 

COSTS - The system should meet all of the above requirements, 
yet reduce the expenditure required for the instrumentation system 
as compared with current methods. 

. . . - . 4nKn«c*AA**** _ - - _ _ 

The system described in this report has been developed to prove 
that the concepts selected meet the above design requirements. 
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3. INSTRUMENTATION SYSTEM 


The syetem described and constructed under this phase meets 
the objectives stated In Chapter 2. The instrumentation system 
can be broken up Into four parts: 

1) Data Management; 

2) Transducers; 

3) Power Supply; 

4) Pilot Control. 

The package Is shown In the block diagram of Figure 3.1. The 
system Is used In two forms: airborne for recording of flight data 

(Figures 3.2 and 3.3), and the ground-based portion for data trans- 
fer to the data reduction computer (Figures 3.4 and 3.5). 

Installation of this system Is straightforward and requires 
no permanent modifications to the airplane. The major modules are 
shown installed in the KU-FRL's Cessna 172 in Figures 3.6 and 3.7. 
The other components are shown installed on the airplane as they 
are described in Section 3.2. It is seen that the major modules 
are essentially strapped into the cabin compartment. The transducer 
module does require a more rigid attachment and is, therefore, held 
firmly in place by clamping it to the seat tracks. 

Following is a detailed description of the instrumentation 
system, as well as the trade-offs considered in its design. 
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Figure 3.1 Overall system blocR diagram 
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Figure 3.2 Major components of the airborne system 
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Figure 3.4 Data transfer system 



Figure 3.5 Block diagram of data transfer system 


i3 













ORIGINAL PACr 

black and white photograph 



Figure 3.6 Battery and computer module installation 


Figure 3.7 Transducer module installation 
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3.1 D ata Manaawtnt 


It waa dacldad to uaa a mlcroprocaaaor controllad data managamant 
syatan. Using a conraerclally available computer simplified the design 
task, as well as reducing overall system complexity and cost. Also 
with this type of controller, versatility is easily achieved, especially 
if programs are stored on cassette tape rather than in computer hard- 
ware. The trade-offs considered prior to selecting this data manage- 
ment system were those of analog vs digital data storage, and airborne 
recording vs telemetry. Following is discussion of these trade-offs 
and a detailed description of the system constructed. 

In the past, most on-board systems made use of analog recording, 
due primarily to the high cost and complexity of digital systems. 

In recent years, however, progress has been made in the digital field, 
resulting in small, inexpensive, and reliable digital devices, most 
available as solid-state Integrated circuits. The recent advances 
in digital electronics technology have reduced both the complexity 
and cost. Coupling this with the lower likelihood of error in 
digital systems, it was decided to use a totally digital system for 
this package^ 

In the past, telemetry has been a much-used means of transmitting 
data to be recorded on the ground. Telemetry has an Important place 
in aircraft flight testing, specifically in high-risk operations 
(such as flutter testing, spin testing, etc.). Its major disadvantages 
are the requirement of a ground station, and the associated high cost 
and complexity. Telemetry, however, has been primarily used in the 
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past, due to the large elea, complexity, and inaccuracies of the 
older recording media. Many improvements have been made in this 
regard with the Introduction of small, reliable cartridge and 
cassette recording systems. This improvement is largely attributed 
to the recent advances in solid state electronics technology. For 
this system, on-board recording, making use of a digital cassette 
recorder, has been chosen. 

The heart of the unit constructed is a Rockwell AIM 65 micro- 
computer. This is coupled through a Rockwell expansion interface 
to the other components. The other tvo auiijor components of the 
airborne package's recording system are the Datel MDAS-16 multi- 
plexer and analog-to-dlgital converter, the TEAC MT2-02 digital 

h 

cassette tape transport, and RS232 Interfacing port. These are 
shown in the block diagram of Figure 3.8. 

The AIM-65 is an interactive single board computer using an 
8-blt 6502 microprocessor. Contained on the computer board is 4K 
bytes of memory, as well as a monitor and symbolic assembler. 

(An 8K BASIC programming ROM is also available for this computer.) 
A 20 character display, 20 column thermal printer and alphanumeric 
keyboard allow the user to interact with the computer. Two appli- 
cation connectors Increase the computer's versatility. One allows 
interfacing to audio cassette recorder and other computer terminals. 
The second allows adding an expansion interface which facilitates 

* 

RS232 - serial interfacing standard. 

ROM ■ read only memory. 
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Figure 3.8 Computer module block, diagram 


additional features to be adapted to the standard computer. 

These features provide an easy-to-use Data Management Controller. 

' The user Is able to easily program the computer using the 
symbolic assembler and monitor functions provided. Programs 
presently are stored on the audio cassette recorder. Using the 
additional ROM slots on the computer, or the addition of a ROM 
board on the expansion Interface would allow regularly used 
programs to be permanently placed In the system. 

The AIM-65 Is co» d, through the expansion board, by use of 
the MDAS-16, to the transducer package. The MDAS-16 is a 16-channel 
multiplexer coupled with a 12-blt analog-to-dlgltal converter. 

This unit has the capability of addressing channels as desired 
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Fi> 5 ure 3.9 MDAS-16 data acquisition card 


UHIGINA. i/v. 

ULACK AND VVHITt PhOfCXiKAFM 

(either rendonly or eequentielly) , using « microprocessor controller 
Voltage input ranges can be selected (-) volt to -^3 volt was chosen 
for this syrtem). Tlie unit has a 50 KHz through-put rate with 20 u 
sec access time per channel. The MDAS-lb is shown in Figures 3.9 
and 3.10. The MDAS-16 does require calibration. This procedure 
is described in detail in Reference 3. 

The other major component of the data acquisition system 
is the TEAC tape transport (sec Figure 3.11). This unit is a 
low-cost magnetic tape unit designed specifically for digital 
applications. It makes use of standard audio type cassette tapes 
for data storage. All interfacing required is Included in the 
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pAckag*. Input raqulrtmcnts arc TTL'^-coRipaelbla; and tha Capa unit 
raqulraa only control algnala, provldad by tha AIM-65 microcomputer, 
and parallel data Input. All detailed control functlona required 
by the tape unit are handled on board by the unit for both recording 
and playback. Only aimple control signale are required to initiate 
the various functions. 

The data management system is also used for data playback and 
transfer to the ground-based data reduction computer. It was decided 
to use the same recorder and computer system for playback of data 
and in-flight recording. This avoids possible problems due to mis- 
match of tape drives and also reduces overall system costs. An 
Interface system compatible with stvandard computer RS232 ports was 
designed and constructed. A hard wire connection, or use of a modem 
through the telephone can thus be utilized. This type of interface 
allows data transfer to virtually any computer. A program on the 
AIM-65 controls the TEAC tape transport and sends the data over 
the line to the other computer. Once all the data are on the 
other computer, the Rockwell system is no longer required in the 
data reduction proces>s. (See Chapter 3 for a complete description 
of the data reduction process.) 

3.2 Transducers 

At the outset of this ^ ogram, it was decided to keep the 
number of containers in the total system to a minimum. Thus, most 

transducers, as well as their required signal conditioning, are 
* 

TTL - Transistor-Transistor-Logic: Electrical standard. 
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contained In the transducer module. This module contalna all 
filtering, all voltage regulation, and the transducer pallet. 

It was not possible to place all transducers In this module, as 
measurements such as control positions and outside air conditions 
were required to be measured. Following is a description of the 
methods used for selecting transducers required, as well as 
descriptions of the actual equipment selected. 

The primary Input to aid In the selection of the parameters 
to be measured was the literature describing the data reduction 
methods to be used (References 4-20). The transducers discussed 
in the references above are summarised in Table 3.1. Discussion 
with personnel at NASA, Dryden Flight Research Center, was the 
secondary input for transducer selection. The transducers selected 
allow optimal use of the data reduction technique considered 
(basically a maximum likelihood parameter estimation method; see 
Chapter 5 for a detailed description of this method). 

The literature (References 4-17 and 20) was also used as the 
primary reference for selection of transducer accuracies required. 
The results are summarized in Table 3.2. The transducer ranges 
were selected after discussion with the general aviation manu- 
facturers (the secondary reference), and consideration of the 
performance characteristics of this class of airplane. 

The ranges and accuracies required for the various transducers 
selected are summarized in Table 3.3. 
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Table 3.1 Transducers used in various flight test programs 
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TABLE 3.2 Transducer accuracias used In various flight test programs 



DELFT 
ref. 4-11 

ECKHOLD 
& WELLS 
ref 20 

KLEIN 
ref 17 

ILIFF & 
MAINE 

ref 12-16 

SELECTED 

A 

X 

.001 g 

.002 g 

.005 g 


.002 g 

A 

y 


.02 g 

or 


.002 g 

A 

z 


.02 g 

2 Z 


.002 g 

0 


1/2 “ 

.2® or 


.5® 



1/2 ® 

2 % 


.5* 

P 


.15®/sec 

.2“/sec 


.5* /sec 

q 

.02®/se< 

.15“/sec 

or 

(U 

iH 

(0 

U 

(0 

.5“ /sec 

r 


.15*/sec 

2 % 

rH 

rH 

0 

•5* /sec 





(4-1 

o 




.4" 

.2’ 

►e 

tH 

.5" 





O 




.4" 

or 


.5“ 

«R 


.4" 

2 Z 


.5” 

T 


2<* F 



2* F 


.1 m 
160 ft 

10 ft 



10 ft 



5 knots 

2 knots 


2 knots 
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Table 3.3 Transducer Accuracy and Range Used 


IIIIIIIIIR^ 

Sansor 

HHBH 

Range 


longitudinal acceleration 

.002 g 

g 

Ay 

lateral acceleration 

.002 g 

^0.5 g 


normal acceleration 

.002 g 

-1.5 g to 4 g 

e 

pitch angle 

0.5" 

±30" 


roll angle 

0.5" 

±30" 

p 

pitch rate 

0.5*/sec 

±50" /sac 

q 

roll rate 

0.5"/aec 

±50*/sec 

r 

yaw rate 

0.5*/sec 

±50"/sec 


elevator position 

0.5" 



aileron position 

0.5" 



rudder position 

0.5" 


T 

temperature * 

2"F 

-65 to +120"F 


static pressure * 

10 feet 

0 to 25K feet 


dynamic pressure * 

2 knots 

40 to 150 knots 

it 

Indicates transducers used to define initial and final conditions. 


During a specific maneuver, T, Pg and need only be measured 
at the start and finish to define the initial and final conditions. 

The other 11 channels require measurement throughout the maneuver 
to determine the dynamic characteristics and analyze stability and 
performance properties of the airplane. 

To select the data acquisition rate required, the following 
factors must be considered: 

- Minimum rate must be higher than the undamped natural frequency 
of the airplane to be tested. 

- Minimum rate must be high enough to avoid time skewing of 
the data points. 

- Minimum rate must be as low as possible to allow economy in 
the recording media and data reduction process. 
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In data analysis » to obtain reasonable representations of the 
frequency response, an acquisition rate of at least five times the 
undamped natural frequency should be used (Reference 21, Volume 1, 

Chapter 6). In the class of aircraft considered for this Instru- 
mentation system, the natural frequencies are of the following order 
(from Reference 22): 

w 0.5 - 1.0 Hz 

"SP 

0 ) 0.01 - 0.03 Hz 

u 0.25 - 0.60 Hz. 

"d 

Therefore, the maximum frequency (w ) requires an acquisition rate of 

SP 

1.0 X 5 ■■ 5 samples/sec. 

This Is the minimum data requirement. 

From References 12 and 14 and discussion with the authors It 
was determined that an acquisition rate of 100/sec Is required to 
avoid time skewing problems. From the practical applications of 
the maximum likelihood estimation method, this rate (100/sec) also 
results In an excess of data that unnecessarily Increases the com- 
putation time and costs. 

Using a computer-controlled acquisition system allows scanning 

of the transducers as rapidly as possible (20 u sec/channel, 220 y 

sec total ) , and then waiting until the next data point is required 
* 

10 samples /sec was chosen for the KU-FRL system, as this then definitely 
meets the minimum data requirement. This rate also seems to be somewhat 

of an acceptable industry standard. 

** 

Values for the KU-FRL system. 
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(0.1 sec later ). These data are temporarily stored In memory and 
then output to the TEAC tape. This technique allows a high scanning 

it 

rate to avoid time skewing between channels (equivalent to 4545/sec 

it 

and a low overall acquisition rate (10/sec ) to provide economy 
and still satisfy the minimum data requ ..rement . 

The transducers were primarily mounted on one pallet. This is 
shown in Figure 3.12. It was possible to include most transducers 
on one pallet with the exception of the 

- pitot tube, 

- temperature probe, 

- static cone, and 

- control position transducers. 



LATERAL 

ACCELEROMETER 


■STATIC AND 
DYNAMIC PRESSURE 
TRANSDUCER 


•LONGITUDINAL ACCELEROMETER 
NORFLXL ACCELEROMETER 


ROLL AND PITCH ATTITUDE GYRO 


FORWARD 


Figure 3.12 Transducer pallet 

it 

Values for the KU-FRL system. 
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The pallet* contained within the tranaducer module* was mounted 
as closa to the center of gravity of the airplane as possible. 

In this flight test program the transducer module has been clamped 
to the seat tracks of the Cessna 172* In the copilots 's position. 

_____ iffthli'trk'kliltit _____ 

Following are descriptions of the Individual transducers used 
In this program. 

3.2.1 Accelerometers 

The accelerometers used In this package are of the force feedback 
(or closed loop) type. This type of accelerometer derives its measure- 
ment from determining the force required to maintain a mass at a zero 
location. This technique reduces the erroKs caused by mass displacement 
and also does not rely on springs (and their associated inaccuracies) 
as do the displacement (or open loop) type accelerometers. The dis- 
advantage to the force feedback accelerometer is its relatively high 
cost. 


It is essential to note that linear (as opposed to vibration) 
accelerometers be used for this type of package. 

The accelerometers chosen are manufactured by Schaevitz Engi- 
neering. Their specifications are shown in Figure 3.13. These ac- 
celerometers are intended for the measurement of linear accelerations 
such as required for guidance control systems, or vehicle ride analysis. 
Both a precision sensor and electronics are integrated into the ac- 
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Spteifieationt L8B Linear 
•t 20*C 

Nominal Nominal 

Natural Output 

Ranga Fraquancy Impadanca 

g Hz kilohma 



±50.0 200 5 


Input Vollaga ^ISV OC nominal 

Input Currant 10 mA OC maximum (0mA OC avaraga) 

Full-Ranga Opan< 

Circuit Output Vollaga -5.0V OC 

Damping Ratio 0.0 typical (0.3 to 1.0 on requast) 

Unaarily (Notaa 1 A 2) ±0.05% ol full scale output 

Mystarasis (Note 2) 0.02% ol full scale 

Resolution (Note 2) 0.0005% of full scale 

CrosS'Axis Sensitivity ±0.002 g par g up to ±10 g ranga, 

(Note 3) Inclusiva 

±0.005 g par g over ±10 g ranga 

Bias Lass than 0.1% of lull seals 

Sanaltiva Asia to 
Cass Alignmant ±1* 

Nolsa Output SmV rms maximum 

Oparatlng Tamparatura -40'C to -t-SS’C 

Storaga Tamparatura -S5'’C to ±105‘C 

Tharmal Cosffieiant 

ol Sanailivlly 0.02% par *C 

Tharmal CosflielanI 

of Blaa 0.002% par *C 

Shock Survival 100 g - 11 ms 

Weight 3 oz. 



Figure 3.13 Schaevitz Engineering LSB series accelerometers 
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celerometer cat*. Interfacing la relatively alnple, requiring only a 
DC Input voltage, and then a measurement of the DC voltage output, 
corresponding to the acceleration sensed. 

3.2.2 Filtering 

The response characteristics of the accelerometers were such 
that they picked up the aircraft vibration caused by the engine. 

The graph of Figure 3.14 shows the airframe vibration charac- 
teristics (measured using the accelerometers as transducers, and 
observing the output on an oscilloscope) as a function of engine 
speed. It la obvious from these curves that the vibration Is 
caused by the engine and Is a function of the engine speed. Also 
of note Is the fact that all the vibration is at a frequency above 
40 Hz. 


A low pass filter with a cutoff frequency at 10 Hz would elim- 
inate this vibration from the measurement signal. Using a two-pole, 
active filter with a response as shown in Figure 3.15 virtually elim- 
inated this unwanted vibrational noise, yet leaves the desired measure- 
ment (occurring In the order of 1 Hz) essentially unchanged. (The 
measurements of the A^ accelerometer are presented as filtered and 
unflltered measurements in Figure 3.16 to show this.) 

In general, as was the case with this instrumentation package 
and the Cessna installation (see Chapter 7) , only the accelerometers 
required any filtering. 


29 


A (g's) OUTPUT VOLTAGE 

« (with input of 12 volts) AMPLITUDE (volts) 



Figure 3.1A Measured airframe vibration 



Figure 3.15 Measured filter frequency response 
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FREQUEMCY (Hz) 


On* drawback of filtering aignala ia tha introduction of a 
phaaa shift dua to tha filter. To counter this problaa, all signals 
should be filtered the same aaomt» thus eliminating the problems 
of the phase shift. 

3.2.3 Attitude Gyro 

Both roll attitude and pitch attitude are obtained from a 
Humphrey VG-24 vertical gyroscope. Full specifications for this 
gyro are shown in Figure 3.17. This is a DC gyro, with potentiom- 
eters for determining the measurement (28 volt DC need for the 
motor, ±5 volt DC used for potentiometer excitation). This gyro 
has operated reliably during both phases of this program. 

3.2.4 Rate Gyros 

A three-axis DC/DC rate gyro package was used for roll-, pitch-, 
and yaw-rate measurement. The advantage of using a three-axis 
package rather than three separate gyros is that alignment for orthog- 
onality upon installation Is eliminated. Of course, failure of a 
single gyro will require the entire package to be removed for repair. 

The gyros selected are of the displacement type (or open circuit). 
Closed circuit (or integrating gyros) will provide better accuracy; 
however, cost of these is approximately 10 times higher. The accu- 
racy of a good quality displacement t3rpe gyro will meet the require- 
ment (see Tables 3.2 and 3.3), especially considering the type of 
power input used (free of oscillations or any high frequencies). 
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5 ninutae aaxiBuo 
uithin O.S* in 9 ainutta 

2 to 9 */ninuta aftar 3 ain. 
uithio O.S* of trua varcical 
OaS*/aia. noainal; raatad on 
t3 l/Z* Scoreby 6 ain. run 
altamating 


SAND AND DUST 

mucus 

EXPLOSION PROOP 


RADIO NOISE INTSRPERCNCE 
SERVICE LIPE 
SMELT LITE 
INfUUTED RESISTANCE 


VCICHT 

SEALING 


Humphrey VG-24 vertical gyroscope 


Figure 3.17 
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The gyros selected are manufactured by Northrop. Voltage input 
required is 28 volts DC» and output voltage Is from ->5 to volts DC. 
The gyros are Northrop G5 submlnlature rate sensors. The gyro package 
specifications are Included in Figure 3.18. 

3.2.5 Control Position Transducers 

Linear displacement transducers manufactured by Space-Age 
Control, Inc., were used to measure elevator position. This trans- 
ducer is depicted in Figure 3.19. Due to Che small size of this 
unit, it was decided to place it externally on the airframe. These 
transducers are Installed as shown in Figures 3.20. 

A novel technique for attaching the control position transducer 
(as well as the total pressure tube and temperature probe) has been 
used. Double-sided foam tape attaches the external devices onto the 
airframe. The mounting technique is depicted in Figure 3.21. The 
mounting method was first tested in the KU-FRL subsonic wind tunnel 
for wind speeds up to 119 mph. The tests in the tunnel were run 
for periods of up to 4 hours, with no degradation in rigidity of the 
mount (see Reference 23), The method has proven to give excellent 
results in the flight test program. The tape used is 3M number 4265 
neoprene foam, the properties of which are included in the table 
on Figure 3.21. 

It was anticipated that the mounting locations for the control 
position transducers would result in non-linear calibration curves. 
However, the calibration curves appeared to have a linear character 
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ORIGINAL PAGE IS 
OF POOR QUALITY 




Walght 

Outlln* dlaanaioaa 
Pe«ar Input 
Input volcaga llslts 
Full*acale output 
Output Ittpadauca 
Output load raalatanea 
Rlppla 

Zaco rata satting 
Input ranga 

(roll/pltch/yaw) 

Maxlaua Input rata 
Output: voltaga 

(at ovarranga llalta) 
Output atabillty 

(Input voltaga varlatlona) 
Rapaatablllty 
Thra ahold 
Raaolutlon 
Hyacarlaia 

Oparatlng taaparatura 
Taaparatura aaaaltlvlty 


2.0 lb. (nax) 

3.7S X 3.7S X 2.13 In 
13 w. (mx)(31 vde) 
28±3vde 
±S vdc 

3000 ohaa (nax) 

SOOK ohna (nonlnal) 

23 nv.paak-paak(Bax) 

±1/2 2 FS 
30*/aac 

600*/aac 
±7 vdc 

1/2 X FS. 

1 X FS. 

0.01 */aac 
0.01 */aac 
0.1 */aac 
0 - 160 *F 

Zaro output IZ FS/100*F 
Seala Factor 3Z F3/100*F 


Warn up tlna 
Motor aecalaratlon tlna 
Clabla datlactlon angla 
Aecalaratlon sanaltlvlty 
Llnaar 
Angular 
Llnaarity 

Sarvlea Ufa 
Inaulation raaiatance 
Daaptng ratio 
Natural fraquancy 
Envlronnanta 
Shock 
Vibration 

Storage taaparatura 
Radio Intarfaranca 


Figure 3.18 Northrop 3 


axis rate sensor 



10 aln 

30 aec (nax) 

±2* typical 

0.05 */aac/g - 
0.08 */aac/rad/aae 
1/2 Z FS. 0-1/2 acala 
2 -Z FS, 1/2 - FS 
100 hr (typical 16000 hr) 

10 magohaa (min), 30 vdc 
0.3 to 0.9 
33 H± (min) 

250g peak sawtooth, 3 raaac 
0.1 g /H±, 20-2000 Hz 
-65 - 200 *F 
MIL-I-8161D 
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ORIGINAL' PAGP 
BLACK AND WHITE PHOTOGRAPH 


TECHNICAL INFORMATION 



Potentiometer life, 1-turn units 1,000,000 cycles 


3-turn units. 600,000 cycles 
or 900 hours at rated power 
Cable static tension at zero extension 10-16oz 
Op. temp , 85*^F to ■♦ 255°F 
Resistance, 1000 Q - Other esistances avail- 
able on special order are 
5 . 10, 20, 50, 100, 200, 
BOO.g^SK, 10K, & 20K. 
45K {3 turn only), 50K and 
100K (1 turn only) 
Standard pots are used unless otherwise specified. 
Specials are available on special order only. 


Standard Special 

Resistance t3% ±1% 

Linearity (3 turn) ±0.25% iO.20% 

Linearity (1 turn) tO.5% ±0.35% 


Max. current at 155°F (ambient) is 31.6 milliamps 
Max. voltage across coil is 31.6 volts 
Power rating, 1.0 watts at 155°F derated to 
0.0 watts at 255°F 

Insulation resistance, 1000 megohms min. at 

500 VDC 

Dielectric strength, 1000 volts RMS min. at 

60 CPS 



SAC Linear Displacement Trans- 
ducers (LDT) consist of an extension 
cable, spirally wound on a spring- 
loaded rewind drum, which is coupled 
to a precision, wire-wound, rotary 
potentiometer. The cable end is at- 
tached to the object whose movements 
are to be monitored. As the cable is 
extended or retracted, the cable drum 
rotates the potentiometer wiper, vary- 
ing the voltage at the wiper tap (No. 
2) of the potentiometer. The voltage 
may be measured to reflect the posi- 
tion, direction, or rate of motion 
of the object attached to the cable. 


Figure 3.19 Space Age Controls linear displacement transducer 
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WING/HORIZONTAL STAB . /VERTICAL STAB 
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Figure 3.20(c) Rudder and elevator control position transducers 



NOTEl- 

-llghcly sand surface of airplane 
-clean with Isopropyl alcohol 

-surface must be room temperature during attatchment 
-fair with duct tape 


3M #4265 -DOUBLE COATED NEOPRENE FOAM T>J*E 


Adhesive A-20 Firm Acrylic 

Thickness 3/64 In. 

Tensile 60 psl 

Static Shear 66 psl 
Temp max. 225 *F 

Temp min -20 *F 


Figure 3.21 Mounting technique for external devices 
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(linear regression correlation coefficient of between 0.9976 and 0.9998) 
for the mounting locations used. 

3.2.6 Static and Dynamic Pressure Transducer 

A BSD Instruments Company 2504 series transducer (see Figure 
3.22) was used for the static and dynamic pressure measurement. 

This device Includes Its o«m signal conditioning and converts 
the pressures to electrical signals utilizing semiconductor pressure 
transducers. Semiconductor transducers are largely affected by the 
ambient temperature; the B&D unit allows for this by heating the case 
and maintaining a constant temperature. 

The pitot tube was designed and constructed according to 
Reference 24, (See Figure 3.23.) The pitot tube is attached to 
the underside of the wing (see Figures 3.24) using the foam tape 
method shown in Figure 3.21. The pitot tube allows a high angularity 
of the flow and still provides true readings. The distance from the 
wing is such that the tube is out of the boundary layer and thus 
provides a true total pressure reading as long as the pitot tube 
axis Is close to the direction of airflow (±15°). The tube is 
mounted along the wing, halfway between the propeller arc and the 
wing tip (see Figure 3.25). This location minimizes flow effects 
due to the propeller slip stream and the wing tip vortices. 

For the accurate measurement of static pressure, a trailing 
static cone is recommended (see Reference 25). Initial flights 
showed difficulty in deployment of the static cone after takeoff. 
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Figure 3.22 B&D Instruments 2054 pressure transducer 
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Figure 3.24(b) Pitot tube mounted on airplane 
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Figure 3.25 Pitot tube mounting location 

The cone has not been used and is not essential if only stability 
analysis is performed. The airplane static system is sufficient for 
stability analysis; however, a more accurate method would be required 
for any performance testing. 

3.2.7 Temperature Transducer 

An Analog Devices Company Semiconductor temperature transducer 
was used tor measurement of air temperature. Specifications are 
shown in Figure 3.26. The transducer is mounted in a probe, as 
shown in Figure 3.27. The temperature probe is mounted the same 
way as the pitot tube, using the double-sided tape method. The 
temperature probe is shown mounted on the airplane in Figures 3.28. 
The location of the probe is identical to that of the pitot tube, but 
on the opposite wing. 
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Mo<l«l A0S90H 

AbsoluC* tUxiaua laciots 
forward voltag* 444 v 

Ravarva voltaga -20 v 

Braakdown voltaga(to caaa) t200 v 
lacad caaparaeura ranga -53*C to 4-lS0*C 
Storaga caaparatura -63*C Co -fl55*C 

Laad caaparatura (aoldarlng)4'300*C 
FowarSupply 

Oparacing volcaga ranga +4v to +30v 
Output 

Noaloal curraat (42S*C) 298.2 pA 

Noalnal taap. coaffielant luA/*C 
Calibration arror (42S*C) tO.S*C aax 


Abaoluca arror (racad ranga) 

No axtarnal adjuataant tl,7*C aax 
+25*C ealib arror - 0 tl.O'C aax 

Nonllnaarlcy £3.0*C aax 

Rapaacabllicy tO.l'C aax 

Long tara drift ±0.1*C/aonch aax 

Currant noiaa 40 pA/Hx 

Powar aupply rajaction 
44v <Va <4-5v 0.5uA/v 

+5v <Va <+15v 0.2UA/V 

+I5v <Va <+30v O.luA/v 

Caaa laolacion 10** ohaa 

Effacciva ahunc capacitancalOOpF 
Turn on ciaa 20 ua 

Bavaraa biaa laakaga 

(ravaraa voltaga >10 v) 10 pA 


Figure 3.26 Temperature transducer specifications 
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Figure 3.27 Temperature probe 
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original h.w., 

bLACK AND WHITE PHOlOGRAPh 



Figure 3.28(a) Temperature probe mounted on airplane 
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Th« cr«ndsuc«rt ••l«ct«d h«v« thovni th«c th« b««ic decisions 


regarding specific transducers, ranges end accuracies were correct. 

They have ell proved reliable, with no failures encountered; end 
none required any specialized signal conditioning or difficult 
calibration procedures. 

3.3 Power Supply 

There were two options considered for supplying power to this 
instruments^ ion syster>: 

1) Tap off the aircraft electrical system, or 

2) Carry a separate battery package on the flight. 

Considering option one, using the aircraft power system, offered 
several advantages. These were reduction in size of the instrumentation 
system, and no limited usage time due to battery rundown. It was real- 
ized, however, that there are several voltage standards on the current 
general aviation fleet. This would therefore require either a complex 
voltage control system or several systems to account for the various 
voltages available in the airplanes to be considered. Coupled with 
this is the high cost of voltage conversion systems. Also, modification 
would then be required to the airplane's electrical system to Install 
the instrumentation package. 

It was decided to explore the second option. A suitable re- 
chargeable battery was found, manufactured by Eagle-Plcher. These 
lead acid batteries are sealed, rechargeable, and maintenance free. 

A typical discharge curve is shown in Figure 3.29. These batteries 
when used in a deep cyclic regime (i.e., removing 50-100% of the 


TTnCAL cm. VOLTA6C 


battery ratad capacity prior to racharga) havt a racharga time of 
12 to 20 houra. Thay hava an axpactad lifatima of 100 to 150 complata 
charga/diacharga cyclea, with longar Ufa axpactancy whan laaa than 
lOOZ dapth of diacharga ia usad. Thaae battariaa can alao be used 
in any position. The c> . t of thaae batteries is such that several 
battery packs could be purchased for leas than the price of one regu- 
lated voltage divider required if the airplane electrical system were 
used. 


TYStCAl VOLTAOf CHAKACTIRISTICS AT VARIOUS RATIS OP DfSCHA^OI (70»P) 



DISCHAR6I TIMI 


D rt aww liw OlMiMft* Rata at Varlaai Raltariat Mattlaly RataR CaaaaMy (C) fey faatar tlMKai 
far aBaa«ata~TtM rata at wlilali aa aifiit aaiaara kaar hat*ary matt ba RltchartaR ta ytatR a 
atafal taa baart aaaait .OMC ar .OtR ■ t A.H. s ,77 amparaa. 


Figure 3.29 Battery discharge curve 
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Another advantage of a battery ayatera is stability of the voltage 
supplied to the systen. This advantage stens from two conditions. 

One is the fact that no external loads are on the power supply; thus, 
the power being used is steady and unchanging. Second is the fact 
that no ripple or noise will be in the power supplied. With a ship- 
supplied system, ripple will be present in the voltage system due to 
the means of supplying power (from the generator or alternator sys- 
tems). This steady voltage supply, and the lack of ripple when the 
batteries are used, results in transducers being able to normally 
exceed their advertised specifications. 

The voltages required for the complete on-board data acquisition 
system are shown in Table 3.4. Batteries were selected to match the 
power requirements at the various voltages. The '.ring schematic, 
as well as the specifications of the batteries selected, are shown in 
Figure 3.30. The batteries allow a minimum of 3 hours running time 
between recharge. (The 12-volt battery supplying the TEAC tape drive 
is discharged first.) 

The biggest disadvantage when batteries are used is that or 
weight. The battery module, complete, weighs 60.5 lbs. This Is 
the heaviest component in the entire system (see the ‘’ble of 
Figure 3.2). Total weight of the entire Instrumentation system 
including all cables is 132.4 lbs. This system weight is not a 
problem for the majority of general aviation airplanes. 
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Table 3.4 Power Requirements fo Data Acquisition Package 


BATTERY VOLTAGE 

REGULATED VOLTAGE 

REQUIREMENT 

+36 

+28 

Heater (Pg,Fjj) 

Gyro motors(6,(^,p,q,r) 

+24 

+15 

Accelerometers (Aj^, Ay. Ajj) 
Filters, MDAS-16 

+12 

TEAC tape drive 

+12 

+5» 5 

Pg and Pp reference voltage 

+5 

Potentiometers 

AIM 65 computer 
Temperature transducer 

-12 

-5 

Potentiometers 

-24 

-15 

Accelerometers (Aj^,Ay , A^^) 
Filters, MDAS-16 



* maximum current requirement 




NOMINAL CAPACIl’Y 

HRR? 

:nsion5 

(INCHES) 


BATTERY 

NUMBER 

NOMINAL 

VOLTAGE 

20 HR 

10 HR 

5 HR 

1 HR 

LENGTH 

WIDTH 

HEIGHT 

TO 

TERMINAL 

I7EIGHT 
(LB) _ 

CF12V20 

12 





6.51 

4.91 

6.53 

6.75 

16.2 

CF12V15 

12 

li.O 

14.5 

13.0 

KBI 

7.22 

3.34 

6.50 

6.75 

12.8 

CF12V8 

12 

8.Q 

7.7 

7.0 

5.0 

6.00 

4.00 

3.75 

3.97 

7.;: 

CF12V1.5 

12 

-±11 

1.4 

1.3 

0.9 

7.02 

1.33 

2.40 

2«60 

1.9 


Figure 3.30 Battery module schematic and specifications 


47 































3.4 Pilot Control 


The pilot controls the Instrumentation system using a box which 
can be placed on the seat beside him (see Figure 3.31). Thu control 
box performs essentially the same £. motion as the ground keyboard, 
the switches on the box replacing the keys (which are really Just 
momentary contact switches). The controls are described below. 

3.4.1 System Control 

This consists of three switches. 

First Is the ’’INITIALIZE" tape switch. This is a momentary 
contact switch which is used only after insertion of a fresh data 
tape. Thl' function prepares the data cassette to accept data. 

Second, the "RUN/STBY" toggle switch is used to control when- 
data is being recorded. In the STBY position the system is non- 
active. In the RUN position, data is recorded. There are two 
of these switches, one of which is located on the pilot control 
wheel and the other, on the pilot control box. 

Third is the "REWIND" switch. This is used at the end of a 
cassette or flight. Activation of this switch places an "end" 
mark on the data tape and rewinds the tape. 

3.4.2 Transducer Readout Control 

A high-impedance analog voltmeter is provided to the pilot 

so that he can observe a particular transducer as he requires. 

The meter’s installation is shown in Figure 3.32. A rotary switch 
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mm 






TRANSDUCER READOUT CONTROL 


INITIALIZE 


RUN/STBY 


REWIND 


Transducer monitor 





(on th« pilot control console) controls the slgnel which is observed 


This feeture is also used to verify that all transducers are 
operating correctly prior to a test flight. 
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4. GROUND COMPUTER SYSTEM 


The MMLE data reduction process described In this report 
requires a powerful computer capable of being programmed In a high 
level language. Phase I (Reference 2 ) pointed out the requirement 
for a computer system operating under a compiled language. This re- 
quirement Is due to the lengthy execution times associated with In- 
terpretive languages. (A Hewlett Packard 9825 was used In Phase 1, 
programmed In Interpretive Basic.) This chapter presents a bench- 
mark process which has been used to evaluate the capability of the 
computer systems to perform the data reduction tasks. Also, a de- 
scription of the selected computer Is presented. 

A two-step evaluation process was used. The first program 
In this process, the INTEGER SPEED ROUTINE, Is short a ad easy to 
Implement and gives a ball-park speed estimate. Secondly, the 
FLOATING POINT SPEED ROUTINE Is a lengthier program, more closely 
resembling the operations performed in the MMLE process. These 
programs are described below. 

4.1 Integer Speed Routine 

This Is a short, easy-to-lmplement program giving a rough bench- 
mark of the operating speed of computers. The Idea for this routine 
was originally conceived in Reference 25. A listing Is presented 
in Table 4.1. The program does not realistically reflect the MMLE 
data reduction process, but it can be easily implemented in virtually 

*1IMLE - Modified Maximum Likelihood Estimation (see Chapter 5). 
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Table 4.1 


Integer Speed Routine 
(Fortran Listing) 


10 


DO 100 M-5,10000,2 

20 


I-M/2 

30 


DO 200 K-3,I,2 

40 


J-(M/K)*K 

50 


IF(J.EQ.M) GO TO 100 

60 

200 

CONTINUE 

70 


PRINT, M 

80 

100 

CONTINUE 

90 


STOP 

100 


END 


any language on most computers in little time. This increases the 
ease with which a benchmark can be run and gives a ball park estimate 
of a computer's speed performance. The results of this speed compari- 
son are presented In Table 4.2. 

For evaluation of this data, It was assumed that once through 
this program was equal to two Iterations of the MMLE routine. There- 
fore, to obtalr the desired data reduction time through MMLE of 
5-20 minutes, the Integer Speed Routine needs to run at 2-8 minutes 
on an acceptable computer. From Table 4.2 it is seen that all accept- 
able computers had both a compiler, which compiled down to machine 
code, and a floating point hardware package. Also, all acceptable 
computers were either using 16-bit microprocessors or could be con- 
sidered main frame machines. It was obvious that current 8-bit 
microcomputers would not be capable of performing the data reduction 
task in any reasonable time frame. This is evident by the fact 
that the AIM-65 (using a 6502, 8-bit microprocessor) could not meet 
the speed requirements even in assembly language. 

This study narrowed the number of acceptable machines considerably. 
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16-BlT MICRO 8-BIT MICRO PROCESSOR 


Table 4.2 Integer Speed Comparison 


MACHIHI 

LAMCUACE 



s 

s 

s 

£ 

HRStMlNiSEC 

ACCEPTABLE I 

AIM 65 

BASIC (PRINTER OUTPUT) 

■ 

■ 

■ 

4: 14: 44 

■ 


ASSEMBLY (LED OUTPUT) 




0: 21: 36 



ASSEMBLY (PRINTER OUTPUT) 




0: 31: 40 

1 

TRS SO 

LEVEL I BASIC 




7: 12: 27 



LEVEL II BASIC 




6: 31: 10 

H 


ASSEMBLY 




0: 21: 55 



FORTRAN 




0: 54: 18 



MODEL U BASIC 



^^y 

3: 15: 00 

■ 

APPLE 11 

INTEGER BASIC 




2: 24: 31 



FLOATING POINT BASIC 

y 

y 

B 

3: 56; 23 

■ 

TERAC 8510 

PASCAL (COMPILE TO P CODE) 

■ 

yy 

■| 

0: 30: 35 


TEKTRONIX (4052) 

BASIC 




1: 23: 00 


HP 9825 

BASIC 




1: 41: 17 


HP 1000 

FORTRAN RTE IV B (CRT OUTPUT) 




0: 01: 23 

/ 


” (NO OUTPUT) 



y^l 

0: 00: 48 

/ 


FORTRAN RTEM (CRT OUTPUT) 




0: 00: 57 

/ 


” (NO OUTPUT) 




0: 00: 44 

/ 

IBM SERIES I 

FORTRAN (NO OUTPUT) 




0: 01: 30 

/ 


" (PRINTER OUTPUT) 




0: 04: 30 

j 

POP 11/34* 

FORTRAN (RSX 11 M) (CRT OUTPUT) 




0: 07: 10 

j 


" (PRINTER OUTPUT) 



y^l 

0: 11: 20 


MINC 11/23 

FORTRAN RTll-IV (CRT OUTPUT) 




0: 03: 36 

/ 


" (DISC OUTPUT) 




0: 03: 29 

/ 


FORTRAN RTll-IV PLUS (CRT OUTPUT) 




0: 03: 10 

/ 


" (DISC OUTPUT) 

y 

HI 

BB 

0: 03: 00 

/ 

HONEYWELL 60/66 

FORTRAN 

■ 

n 

n 

0; 00: 44 

/ 


PL/1 


mm 

BB 

0: 02: 13 

/ 

CDC CYBER 70 

FORTRAN (NON OPTIMIZED) 


D 

BB 

0: 00: 39 

/ 


FORTRAN (OPTIMIZED) 


By 

BB 

0: 00: 37 

/ 

IBM 370-148 

PL/1 (OPTIMIZED) 

H 

y 

B 

0: 01: 19 

/ 


Tha POP 11/34 was oporaclng in a mulcl-user mode. Its performance la asclmated to be 
approxinatel/ 2-3 times faster than the 11/23 series computer In single-user mode. 


4.2 Floating Point Speed Routine 

To more closely resemble the MMLE data reduction process, yet 

still use a simple-to-impleraent program, the routine shovm in Table 4.3 1 

\ 

was developed. The program is made up of floating point matrix i 

mathematics, which is what MMLE primarily contains. j 

3 
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Tabl« 4.3 Floacing Point Speed Routine 
(Fortran Listing) 


10 

20 

30 

40 

50 

60 

70 

80 

90 

100 

110 ' 

120 

130 

140 200 
150 
160 
170 
180 
190 

200 300 

210 

220 

230 400 

240 

250 

260 

270 

280 100 
290 
300 
310 


REAL A(20,20) .8(20,20) ,C(20,20) ,E(20.20) ,T(20,20) ,0(20,20) ,F 
INTEGER I,J,K,M 
PRINT, ’’START” 

F-. 098625 
DO 400 M-1,40 

DO 200 1-1,20 
DO 200 J-1,20 

E(I,J)-0 

A(I,J)-F*I*J 

B(I,J)-F*I 

C(I,J)-F*J 

D(I,J)-F 

T(I,J)-0 

CONTINUE 
DO 300 1-1,20 
DO 300 J-1,20 
DO 300 K-1,20 

T(i,J)-T(I,J)+(A(I,K)*B(K.J)) 

E(I,J)-E(I,J)+(E(I,K)*D(K,J)) 

CONTINUE 
DO 400 1-1,20 
DO 400 J-1,20 

E(I,J)-E(I,J)+T(I,J) 

PRINT. ”E-” 

DO 100 1-1,20 
DO 100 J-1,20 

PRINT, E(I,J) 

CONTINUE 
PRINT, M 
STOP 
END. 


The program approximates one iteration of the ^^E method. 

This is indicated by the 48 minute run time on the Hewlett Packard 
9825, which requires approximately 50 minutes to perform one iteration 
of the MMLE program. In order to deem a computer acceptable for the 
MMLE process, it must be able to complete the floating point speed 
routine in the order of 1-5 minutes. It was decided that an MMLE 
execution time of 5-20 minutes would be acceptable (assuming 5 iter- 
ations) . 
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The results of this test are* presented in Table 4.4. It seen 
that the 16-blt machines tested, operating in compiled Fortran, meet 
the speed requirement. 


Table 4.4 Floating Point Speed Comparison 


MACHINE 


MIN: SECS 

HP9825 

(BASIC) 

48:15 

HONEYWELL 60/66 


0:20.6 

HP 1000 

(NO OUTPUT) 

1:08.7 


(DISC OUTPUT) 

2:07 

IBM SERIES 1 

(DISC OUTPUT) 

0:58 

MINC 11/03 

(DISC OUTPUT) 

5:35 

MINC 11/23 

(DISC OUTPUT) 

4:00 


4.3 Description of System 

The results of the benchmark evaluation left several computers 
that were deemed acceptable. To select the best machine for the 
KU-FRL requirements, the following factors were also considered: 

- Memory expansion capability 

- Floating Point Hardware available 

- RS232 ports /IEEE 488 ports installed 

- CRT Graphics capability 

- Hard/Flexible disc storage 

- Programming languages available 

- Users group existing 

- Delivery 

- Cost 


•k 

Industry interfacing standards 
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Evaluating the acceptable computera* the DEC MINC 11/03 computer 
waa aelected aa beat meeting the requlrementa. A deacrlptlon followa. 

The MINC 11/03 la ahown In Figure 4*1. The block diagram of 
Figure 4.2 ahowa the baalc featurea and acme of the optlona available. 

The computer uaea a 16~blt DEC LSI 11/03 proceaaor, capable of 
addressing 64K bytes of memory, and contains a floating point hardware 
package, 4 RS232 ports, and an IEEE 488 port. 

Data and program storage Is handled using the dual RX02 flexible 
disc drives. These use 8" flexible discs, capable of holding 300 K 
bytes of information each. 

Computer and program Interaction Is handled using the DEC-VT 105 
graphics terminal. This permits Inputting and outputting of data, 
as well as allowing graphical representation of the flight test 
results. 

The RS232 ports are used for Input and ouptut of the data. 

Four are provided. One Is used for the VT 105 terminal, two are 
configured to allow data transfer from the Rockwell AIM-65, and 
one Is used to control a hard copy printer. 

The IEEE 488 port allows ease of Interfacing to many Industry 
standard components. Planned future use of this port Is for a 
hard copy plotter for analysis and report quality plots of flight 
test data. 

* 

DEC • Digital Equipment Corporation. 
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Th« atandard MINC comaa with BASIC languaga aoftvara. Tha 
KU-FRL packaga haa tha RTll-FORTRAN IV aoftvara option. Thla 
version of FORTRAN allows compiling programs to machine lavel» 
which was datarminad necessary to perform tha data reduction cask 
as indicated in Sections 4.1 and 4.2. 

Tha MINC computer has been found capable of performing the 
function intended. The MMLE process takes approximately 20 minutes 
for 5 iterations, which is close to the prediction of Section 4.2. 

It is recomoended that a hard copy nrinter and plotter be 
added to the standard MINC to make it a complete data analysis 
system. 
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5 . DATA REDUCTION METHOD 


This chapter describes the data analysis procedures used for 
longitudinal and lateral stability analysis. The overall method 
Is best depicted via the flow chart shown in Figure 5.1 

For this phase the system described In Chapter 3 was used for 
airborne data acquisition. The KU-FRL's DEC-MINC 11/03 microcom- 
puter was used for all further data processing. This computer Is 
described in detail In Chapter 4. Segmenting the various data re- 
duction programs Into the blocks as shown In Figure 5.1 allowed 
effective data analysis. 

This section describes theoretical aspects of the computer 
programs used. Flow charts and program listings are Included as 
Appendix A. 

5.1 Data Acquisition 

This program is used as part of the airborne data recording 
system. It is written on the AIM-65 in machine level language to 
allow rapid execution. The program controls the MDAS-16 module, 
as well as the TEAC cassette recorder. (See Appendix A.l for flow- 
chart and listing.) 

The orogram has three control inputs, which are located on 
the pilot control console. The first is the " INITIALIZE’* tape 
button. This is used for getting the data cassette ready to record 
the signals. It is used only once per data cassette. This command 
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-DISPLAYED ON 
GRAPHICS CRT 


-CORRECT FOR 
CENTER OF GRAVITY 
LOCATION 

-REMOVE DATA 
ERRORS 


-DISPLAYED ON 
GRAPHICS CRT 

-HARD COPY PLOTS 
BY TRANSFER TO 
HEWLETT PACKARD 
COMPUTER 


Figure 5.1 Data processing flow chart 
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rewinds the tepe (If required), edvences the tape pest the beginning 
of the tape hole, and then writes a beginnlng-of-tape file mark. 

Second, the "RUN/STBY" toggle switch la used to control the 
recording of data. Placing this switch In the RUN position begins 
the data recordlr.g process. The computer then sends control to the 
MDAS to sample the P^ and T channels. These are sampled 10 times 
each and then output to the TEAC cassette drive In one block. The 

» Ay, Ajj, 0, <^ , P, 
q* Z'* <^£» • These data are temporar^y stored In computer 

memory. After a total of 0.1 seconds has elapsed, the computer then 
samples these channels again and also temporarily stores them In 
memory. After 10 such time points are In the computer memory, the 
AIM-65 outputs this to the TEAC In one block; and the process con- 
tinues until the "RUN/STBY" switch Is placed In STBY. Then the 
computer samples the P^, Pg and T channels again and outputs these 
to the cape. After this the system idles, waiting for the next command. 

To reduce the possibility of error, the hignest order bits on 
the measurement channels, are recorded twice. This is easily done, 
as the aaalog-to-digltal conversion comes out as a 12-bit word, 
available as a tri-state output. The AIM-65 operates on the basis 
of 3-bit words; therefore, the 4 highest order bits of the data are 
recorded twice, resulting in two 8-bit words. These higher order 
bits are compared on readback as a means of error checking. 

Third, the "REWIND” switch causes an end-of-tape mark to be 
written on the cassette and then rewinds the tape back to the start. 


program then runs through the other channels (A^ 
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This program also keeps track of the run number, which is output 
at the beginning of each run to the cassette. 

5.2 Transfer Data to Ground Based Computer 

This operation requires both the AlM-65 Ltid MINC 11/03 and a 
program for each to allow the two to be coupled. A standard RS232 
serial Interface Is available on both computers. The data can be 
transferred across telephone lines If desired. 

The data, before being transferred from the AIM-'bS, Is checked 
for errors. This Is done by comparing the high order bits, which 
have been recorded twice. A running total of any errors is kept 
and printed out by the AIM-65 on Its display printer. Errors have 
not been significant in number, and therefore no correction is made. 
All errors to date have been caused by poor quality data cassettes. 
Using the qualified cassettes (see Table 5.1 and Reference 27) 
no data errors have been found In the flight data. 


Table 5.1 Qualified Data Cassettes 


Manufacturer 


Part N< 

3M 

Scotch 

834A/1-300 

TDK 

Data Cassette 

HR-850 90C 

MAXELL 

Data Cassette 

M-90 

BASF 

Digital Power Typing Cassetv«i 

52346 


(Qualified as per Reference 27) 
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The AIM-65 program la shown In Appendix A. 2; and the MINC 11/03 
program. In Appendix A. 3. These programs are used to transfer the 
flight data from the TEAC cassette tape to the MINC 11/03 disc. In 
this mode the AIM-65 keyboard Is used for controlling the date transfer 
off of Itself. The MINC 11/03 program loads the transferred data Into 
Its memory and then transfers this date to the data disc. 

5.3 Engineering Conversion 

The first step In the actual data analysis procedure Is con- 
verting the raw data bits Into their corresponding engineering 
units. The process Involved first converts the bit pattern of 
each measurement to the voltage representation. (See Reference 3 
for detailed expl<?>natlon of this process.) Then, utilizing the 
particular transducer calibration curve, this voltage representation 
is converted to the units of the actual motion measured. Resulting 
from this, then, is the transducer measurement In the correct en- 
gii leering unit. (See Appendix A. 4 for program listing.) 

This two-step process is presently required due to the cali- 
bration process utilized in this phasfu Currently, transducers 
are excited using known inputs; and the transducer response is 
measured using a voltmeter. A suggested improvement in this 
process is to bypass the voltmeter, using the digital recording 
system in the calibration process. This improvement is planned 
to be implemented upon construction of the calibration rig sug- 
gested in Chapter 9. 
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5*4 Quick-Look Plots 


The next step In the data analsys procedure Is making the quick- 
look plots. The program of Appendix A. 3 is used to do this. Basically 
this program plots the transducer outputs (uncorrected for C.G. loca- 
tion, etc.) on the graphics CRT. This Is a rapid means of determining 
the portion of the recorded data that has the proper aircraft modes 
excited and is thus suitable for further analysis. Operator Inter- 
action has been minimized to reduce the overall time required for 
this step. 

5.5 Detailed Engineering Conversion 

This program Is used to do a rigorous conversion of the data 
Into the form required for the MMLE technique. Accounted for in 
this procedure are accurate transducer calibrations and instrument 
position corrections. 

The first step in the instrument position correction process 
is to account for the misalignment between the transducers and the 
aircraft body axis. Secondly a correction must be applied to correct 
for the distance from the transducer center of gravity to the air- 
plane center of gravity. The following equations are used. (See 
Reference 28 for a more rigorous presentation.) 

Pg - Pj^ cos (9.J.) + r^ sin (0^.) [5.1(a)] 

•"b “ “Pm 
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\ \ ^ si" <®l) + f - (pq - r) I - (pr + q) 

Ay - Ay - (pq + r) ^ + (p2 + r^) |- - (qr - p) 7 [5.1(b)]* 

A^ ■ Aj. sin (9 ) + k cos (0 ) + (pr - q) | + (qr + p) | - (p^ + q2) | 
where 

B indicates Body axis at airplane center of gravity 
M indicates as Measured by transducer 

I indicates as Installed wrt Body axis at airplane center of gravity 
This step also involved checking for and correcting any obvious 
data errors. If any filtering of unwanted noise is required, it 
would also be done at this stage; however, none has been needed 
to date. The quick-look plots are used as the major aid in this 
process. 

A program listing is contained in Appendix A. 6. 

5.6 Modified Maximum Likelihood Estimator 


The flight data were processed through the Modified Maximum 
Likelihood Estimator (MMLE) developed by NASA (see Referencfis 12-16) . 
This technique has been used by NASA for over 12 years. A simplified 
program (NASA Dryden "BONES" version of MMLE) has been placed on 
the. MING 11/03 computer. The actual program listings are included 
in Appendix A, 7. Described here is the theory used in this tech- 
nique, and some of the assumptions made for the KU-FRL version. 

'k • m m 

Where p, q, and r are required, these are determined by digitally 
differentiating the p, q, and r measurements. 
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5.6.1 Parameter Eatlmatlon 


The MMLE estimator Is an iterative process that determines 
the coefficients of a given set of linear equations describing the 
motion of the aircraft. It does this by comparing the difference 
between actual In-f light measured responses of various states, 
and the predicted responses of these states using an estimate of 
the coefficients. The actual measured control Input Is used as 
the Input for the estimating procedure. The estimated coefficients 
are. updated each Iteration, using the differences as determined 
above. The flow chart below shows the MMLE concept. 
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5.6.2 Mathematical Modal 

The mathematical model used to describe the airplane is 
derived from the small perturbation equations of motion (see Ref- 
erence 22).* These are shown here explicitly, in the non-dimensional 
form. 

- for longitudinal (from Reference 22, Equation 6.1): 


mu - - mgOco.u + q^S(-(C + ^ + (C^ + 2 C^ > F S 

u 11 I a 1 

m(i - U^q) - - mgO.inO^ + 5lS(-(Cj^ + 2C ) ^ - (C^ + C )a - C |§- - - C 6 ^) [5.2(a)] 

u 11 ol al ql 0^ 

* 2'. > r * <V r fir * ir * *t> 

u \ X T T, loiT alqlo- 

u 1 a ^ E 


- for lateral (from Reference 22, Equation 6.2): 


m(v + U^r) - mg^co.O^ + ^l^Sg^ "" S 2 ^ 2 ^ ^ S 5 *A ^ 


^:;P • ^ 2^ ^ 2u7 * 

p p ji r i 


* 'a ^ «R> 


I r-I p - q.Sb (C B + C 6 + C ^+C ^+C i . C 4„) 

zz xz'^ ^1 n, n_ n 2U, n, 2U, n, A n, R 

B Tg p 1 r 1 4^ *R 

Using the definitions shown in Table 5.2, Equations [5.2] can 
be converted to the dimensional form shown below. 


- for longitudinal (from Reference 22, Equation 6.72): 

u ■ - g9cos0j^ + X^u + u + X^a + X^ 6^ 

u E 

w - U,q " - g0sin0. + Z u + Z a + Z*a + Z q + Z. 

1^ ® 1 u a a q^ E 


q*Mu+MT,u+Ma+M-,a + M*a + Mq+M- 5_, 
^uT aT a q^5_E 

u a ^ E 


*The derivatives in Reference 22 are for the stability axes system. 
See Appendix B for conversion to the Body axes used in this report. 
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Table 5.2(a) Longitudinal Dimensional Stability Derivatives 


*u 55^^ («c-i) 

q^scc^ + 2c ) 

X, 

\ ■ 5iT~ (sec*^) 

Vi 1 

-5iS(c„ - c ) 

\ ^5 (Usec-h 




«E " 


^ <ft sec"^) 


~ <se="^) 


“o ■ i ^ (>ei:"^) 

yy 


1 


^ -’1— 
“ yy 


~ (sec~^) 


q.Sc^C 
1 m* 


M* ■ ^ I 

a 21 U, ) 

yy 1 


*S> 

^0. n; (ft sec'^) 


■ ■ 2mUj^* (f‘ set'^) 


qjSCj_ : 


‘q ’ ” 2mU, 


(ft sec 


ij^Sc^C 

”q “ 21 U,'’ (sec"^) 

yy 1 


E yy 


E (sec 




*^E (ft sec 


q.Sc(C + 2C ) 
i m m, 

M " — ^ ^ **lv 

u I U, ®®c ) 

yy 1 

q,Sc(C + 2C ) 

M _ u ^1 

^ I ij (ft ^ sec 

U yy 1 


* from Reference 22, Table 6.3, page 413 



Table 5.2(b) Lateral-Directional Dimensional Stability Derivatives * 



from Reference 22, Table 6.8, page 445 
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- for le ^r«l (from Reference 22, Equation 6.141): 


+ U^r - g*co.»i + Yj8 » + Y^t + Yj^S^ + Y^^«^ 


P - f - L,e + L„p + L,c + Lj + Lj «R 


XX 


[5.3(b)] 


r - 


xz 


zz 


p - N-6 + N_ B + No + Nr + N. 6. + N. 6 


■8- Tj- -P 




Using the concept of state variable theory (see Reference 22), 
Equation [5.3] can be written in the following form: 


[R] (i(t)} - [A] {x(t)} + [B] Cu(t)} [5.4] 

where 

{x(t)} • state vector 

[R] ■ acceleration transformation matrix 

[A] • stability matrix 

[B] ■ control matrix 

(u(t)} ■ control vector. 


Equation 5.4 can be written more explicitly in the form which 
follows : 


- for longitudinal (where [R] - identity matrix) : 


• • 




Ir m 


r 1 


q 


M' M* 

q u a V 


q 




■ m 









E 

u 


0 X.' X' -g co« (e.) 

M Cl 1 


u 

+ 

X/ X . ' X’ 


S 



z + u, 






c 

a 


-3 z' 2’ s — *“77 tin <0,) co« (♦.) 

U, - Z* u a u, - 2* 1 1 

la la 


a 


2*' Zr* z; 


1 

0 


co» 0 0 0 


9 


0 0 6* 

0 


m m 

• • 


m — 


* - 





(See Table 5.3 for explicit definition of these terms.) 


[5.5(a)] 
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for lat«ral: 


i«i A 


— - 


» • 


» • 


» ■ 

p 


L’ Ll l; 0.0 

p r B 


P 


\ *r ® 

r 

• 

M' M' n; 0.0 

P r D 


r 


M * N ' N* 

\ \ % 

« 


■ In(cij) -coB(a^) YJ ^ coB(6^)coa(i^j^) 


6 


V/ If*' K 

« 

L j 


L.O eoB(p. )tan(0. ) 0.0 0.0 

^ * ■ 


♦ 

• m 


0.0 0.0 


1 

k j 


(See Table 5.3 for explicit definition of these terms.) 


To allow determination of states other than the ones contained 
In {x(t)}» the following expression can be derived: 

(y(t)} - {x(t)} + {u(t)} + t“““) * [5*6] 

where 

{y(t)} ■ computed observation vector 

[G] • observation matrix 

[H] ■ observation matrix 

{v} ■ variable bias vector. 

(See Table 5.4 for explicit definition of these terms.) 

The computed observation vector, (y(t)}, corresponds to the 
measured observation vector, shown here: 

(z(t)} - (y(t)} + {n(t)> (5.7] 

where 

{z(t)> ■ measured observation vector " {9, (j), p, q, r, A^, A^, 

^N’ ^A* ^S* ^D’ 

{n(t)} ■ measured noise vector. 

From the terms of Equations [5.4], [5.6], [5.7], the vector 
*From Reference 16 


[5.5(b)] 
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T«ble 5, 3 (a) Longltud^l, Dlncntlonal SC«tc Vector Stability Parlvatlvas 


Z_ + U 




II, - 

z* 

1 

a 


Z 

^ a 

u 

" U.-Z. 

1 

a 

M« 

Z 

a 

+ r~ 

a 


- 2 , 


ala 


-M'g aln(0.) cos ((<),) , . 

M* - — 2 L. , 0 (sec"^) 




M* Z. 
a 6 


EiC / -2. 

^ (sec ) 


"”5 U.-Z* 

E,c E,c 1 a 


pitching moment equation bias 
(aec”^) 

X' - X + X_ (sec"^) 
u u T 

u 

X’ - X (ft sec‘^) 

01 a 

X^ " ^6 

E,c E,c 


X* - longitudinal acceleration 
equation bias (ft sec ) 


II, - Z* 
1 a 


* 1 . 0 ' 


z z 

u u 

u - U^-Z. “ 


Z' ■ = 3 -i 

cx U,-Z‘ U. 
la 1 


(ft‘S 


(sec**^) 


-g 8ln(0j^)cos(<|)j^)» -g sin(0j^)ccs((j)j^) 


U,-Z* 
1 a 


7 * m EilC 

*E.c V==i 


U, 


(sec”^)^ 


E.c , -1, 

if ^ 


Z^ ■ normal acceleration equation 
bias (sec 

• -In* 

0^ ■ pitch rate equation bias (sec ) 


Note: The equation bias terms are used to allow prediction of the complete 

state which is made up of the steady state and the perturbed state. 


Note: With the approximations above, Equation [5.5(a)] is rewritten as; 
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Table 5.3(b) Lateral. Dloenaional State Vector Stability Derivatives 


K 
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Lp (sec“^) 
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‘ (sec 
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N.' - N. (sec"^) 

*^A ^A 
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^A ^1 
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r 1 


Y' 
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I 


‘■'o 


N’ 
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-1 * 

lateral acceleration equation bias (sec ) 


-1 * 

roll rate equation bias (sec ) 


- 2 . * 

rolling moment equation bias (sec ) 

•*2 ^ 

yawing moment equation bias (sec ) 


NOTE: The equation bias terms are used to allow prediction of the 

complete state which is made up of the steady state and the 
perturbed state. 


[R] 


1.0 

I 


xz 


XX 


xz 


zz 


1.0 


1.0 


for I ==0; [R] - identity matrix 

XZ 
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Tabic 5.4 MaCrlcca Used In Observation equation 


LONGITUDINAL (y(t)}^ - {q, U, a. 0. q, k^) 
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(c) - f([Aj, [Bl. [Gl. [HJ. {v}) [5.81 

(wh«r« f lndlc«t«« "« function of") Is d«fln«d at the vactor of 
unknowna. It la this vector that the MMLE method eatlmatea. 

MMLE determlnea the unknowna ({c}) by minimizing the coat function 
given by: 

T 

J - i (z(t) -y(t)}‘’'[D] (z(t) - y(t)} dt * [5.91 

(T, t: Indlcatea time) 

or approximately In the dlacrete case: 

N 

“ (i^T) ^ ^*1 “ ^1^^^°^ ^^1 ■ ^1^ * [5.101 

(where 1 Is the time Index, and N the number of time points) . 

The weighting matrix, [Dl, is used to provide emphasis on the 
various measured states; In other words, to allow greater emphasis 
on the more accurate transducers, or the transducers that are more 
Important to describe the maneuver performed. 

The value of the cost functional, J, Is minimized using the 
Newton-Raphson method. This technique Is an iterative procedure, 
utilizing an estimated value of the vector of unknowns, {c}, and the 
first and second gradients of the cost functional, J, with respect 
to the vector of unknowns, {c}. The equation 

(where L Is the Iteration number) Is used to revise estimates for 
the vector of unknoxms, {c}. The first and second gradients are 
given by; 

*From Reference 16 
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N 


(v> - 2 


[5.12] 




1-1 


1-1 


[5.13] 


The Balakrlshnan modification makes use of the fact that the 
term - y^} approaches zero with convergence and Is thus 

neglected. The expression for the second gradient becomes: 


N 


"c'^> ■ iPT - yi> - yi> 


[5.14] 


After several iterations the cost function converges near some 
small value. At this point, the parameters of Equation [5,5] have 
been modified to obtain their most likely value which results In 
the best fit of the measured states. 

5.6.3 Assumptions Used In Data Reduction 

The following inputs and modifications were made to the MMLE 
method, allowing effective use of the technique on the MING 11/03 
computer. 

Initial estimates of the derivatives In Equation [5.5] were 
obtained using the analytical methods of Reference 22. Although 
the MMLE technique does not require accurate knowledge of these 
derivatives, this procedure does speed convergence. 

The MMLE program usually uses a modified least squares method 
for the first iteration to estimate the derivatives, as an aid to 
*From Reference 16 
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•pcttd conv«rg«nctt. Thli, however, requires meesurenent of most 
of the stetes Indicated In [5.5]. The instrumentation package 
uses only a minimum of transducers, and all the states required 
for this least squares estimate are not measured. Using the least 
squares procedure would result in divergence of the first iteration. 
Therefore, the least squares estimate was nut used, which did slow 
convergence of the derivatives. 

A diagonal multiplying factor allows control over how large a 
change is made to the derivatives after each iteration. Too large 
a value of this factor causes sluggishness in the convergence, and 
too small a value will cause divergence. Further analysis Into 
this factor will Indicate Its optimum value for best convergence. 

The weighting matrix, [D], of Equation [5.9], was chosen after 
analysis of the instrumentation error magnitudes. The first run 
through the MMLE program, with measurements from this instrumentation 
package, provided a weighted error for each measurement state. As 
suggested in Reference 16, the values for the weighting matrix were 
chosen to attempt to equalize the weighted errors. After the values 
for the weighting matrix were chosen for the instrument package, 
they were then left at this for further maneuver analysis. 

5.7 Time History Plotting 

The MMLE reduction method not only produces the estimates for 
the derivatives, but also calculates the estimated time history for 
the various states. This is stored on the data disc by the MMLE 
program. The programs presented in Appendix A. 8 retrieve both the 
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predlctAd tlm* hlatorl«s and th« nwasurad tlm* hlatoriaa and plot 
them togathar on tha graphlca CRT. Thaaa grapha of tha flight 
teat manauvar ara tha vlaual Indlc ition of tha goodneaa of tha 
pradlctad alrplana darlvatlvaa. 

aaaaaa**** 

As Is evident from the many programs provided, the final 
results of a flight test maneuver are obtained only after a multi- 
step procedure. This Is primarily due to the nature of the methods 
being used In aircraft flight testing, as well as the limitations of 
computer technology being used. 
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6. KU-FRL FLIGHT TEST PROGRAM 


Two s«rlea of flight teats have been conducted using the KU-FRL 
Cessna 172. The first aeries, cord 'ted under Phase 1, Is presented 
In Reference 2. Presfmted here < . the basic concepts of the type 
of flight maneuvers required, and results of the Phase II test program. 

6.1 Flight Teat Maneuver 

Traditional flight testing methods have utilized primarily steady- 
state flight paths for data collection. This was due mostly to the 
data acquisition systems available. Unfortunately, this required rn 
highly trained and competent test pilot to obtain realistic and valu- 
able results. 

With the current transducer and acquisition system technology 
available, flight testing need no longer rely on steady-state maneu- 
vers to allow accurate state measurement. This development has re- 
sulted in the newer flight testing metnods utilizing dynamic maneuvers. 

When techniques such as the MMLE are used, the literature (Ref- 
erence 29) Indicates that the nature of the maneuver Is not critical 
to determine the aircraft characteristics. What Is Important when 
using these techniques Is to ensure that the proper aircraft modes 
have been excited. For example, a longitudinal maneuver should excite 
both the short-period and phugoid ?odes of the airplane. This reali- 
zation (i.e., non-critical flight path) leads to the possibility of 
using lesser qualified pilots and still obtaining accurate results. 

All testing done on this program has been don^.- by a pilot who had no 
previous flight test experience. 
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Tht control Inputs prsssntsd In ths trscss of Chsptsr 6.2 
sro typical of ths type of maneuver required. Several frequenclea 
are excited, which tends to Increase the validity of the results 
obtained^. Also the total energy Input is approximately symmetrical. 
In other words, the motion produced In one direction Is offset by 
the motion produced In the opposite direction a short time later. 

6.2 Results of Flight Test Program 

Presented here are results of the Phase II flight test program. 

e 

All flights were done at the conditions of Table 6.1. 


Table 6.1 Cessna 172 Flight Test Conditions 


Wing area (S) 

174 ft2 

Wing span (b) 

35.8 ft 

Inertias 


I 

XX 

1029 slug ft2 


1891 slug ft^ 

^yy 

1092 slug ft^ 

Mass (m) 

59.46 slug 

Weight 

1913 lb 

Center of Gravity (Body Station) 

41.3 Inch 

Mean chord (c) 

4.9 ft 


The plots following show the absence of noise In the measurement, 
as well as the typical maneuver required. 

Estimated by Reference 30 


176 ft/sec 
33.69 lb/ft2 
3000 ft 


Speed (Uj^) 

Dynamic Pressure (q^) 
Altitude 
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The fit of the eetlneted otmtee compared to the actual atatea 
In the longitudinal maneuvera la good (Flgurea 6. 1-6. 3) v Tha only 
state that la )ff <'onalatently la the term, which appaara to be 
affected by a phaaa ahlft. The cauae of thla phaae ahlft has not 
been determined. 

The estimated parameters have been compared with the analytical 

methods of Reference 22 and with flight teat results obtained by 

NASA Langley on a Cessna 172 (Reference 30). This correlation Is 

shown In Table 6.2 for the longitudinal maneuvera. It Is seen that 

there la good correlation between some derivatives, but not between 

others. The best correlation appears to be with the one of run 23B, 

In which the speed derivatives have been held constant for the MMLE 
* 

analysis. This would tend to be the predicted result due to the 
mismatch In the A^ term, which is a major contributor to the speed 
prediction. 

The lateral maneuvers are presented in Figures 6.4-6. 6; and the 
correlation of derivatives, in Table 6.3. The fit of the measured 
and predicted states is again reasonable. Run 11 has the best fit 
as well as the overall best fit to the parameters. Again, however, 
the predicted coefficients are not within acceptable limits. The 
cause of this is not known. 

Observing the rudder trace on Figure 6.6. what appears to be 
rudder float is evident (especially between 10 sec and 12 sec) . 

•k 

This is the same procedure performed by NASA Langley, which makes no 
attempt to predict any speed derivatives. 
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T«bl« 6.2 Coaparlton of R«sult«» Longitudinal 



( ) Aa eoaiparad wltti HAM Langlay raaulea. 

* Hroog alga. 

^ Rafaranca 30, Tabla IV, paga 28, Ka.'iiMa Llkallhood Nachod, 
avaraga at cha two niaa at Full Trta. 

^ Saa Appandla 8 for eoovaralon to tha lody aaaa ayataa uaad 
la chla raport. 

^ Rafaranca 22, Alrplaoa A, paga 190. 

'* ' haa a larga coaponant which cannot ba pradlctad. 

q a 


85 








































I 


8 

B «.e 


« 

-4M,P 



S>aoo.o 
5 «••«» 


c 

w -«a,e 
J aaa.o 
M 
\ 

O 


-O.KO 

M.O 



Ui>. 0 


• 0.0 



Figure 6.5 Flight time history; Flight 19/10/80 Run 26C; Lateral 


07 



SR 


OiO 


-ao.o 


4. O 


a. o la. a 

TZHK ZN SECONDS 


1*. Q 


Fi<?ure 6.6 Flight tiae history; Flight 19/10/80 Run 33B; Lateral 


ao. o 


88 



Table 6.3 Coaparlaon of RaaulCa, Lataral 



STICS riXDi AlRSniD 174 ft/a«'^ 

(120 rVN) 


latlaatlaa 

HaelMd 

cD-m wu 

NASA 

LANCUl'** 

ANALYTICAL*** 

Oroaa 

Uatglie 

(14) 

19U 

1848 

2140 

Cwcar of 
Oraaley 

(lady acaelaa, 
iMhaa) 

41.3 

42.3 

40.3 


nitht M«. I 23/10/80 


Itw No. 


19/10/10 


11 28C 338 

(rig. 9.4) (rig. 4.3) rig. 4.4) 


(IK) (IS) (241) 

•0.402 -<1.491 -0.331 


* (I3X) 

-0.108 0.139 


(341) (9Z) (34S) 

-0.049 -0.047 -0.034 


(2S) (3001) (211) 

-0.042 I -0.419 -0.074 -0.043 


(141) (821) 

0.037 0.008 


(421) (201) (231) 

-0.333 -0.464 -0.439 


( 101 ) ( 10001 ) 

0.009 0.128 


( ) Aa coaparad trith NASA Langlay rc-ulca. 

* Wrong atgn. 

^ Rafaranea 30, Tabla Vtt, p. 32, caaa 34. 

^ Saa Appandta B for eonvaraion to cha Body axaa ayacan uaad 
In chla raport. 

^ Rafaranea 22, Alrplana A, paga 390. 















































Thia manauvar waa parfornad by holding Cha ruddar padala fixadf 
yac a float of 2*-3* la aaan in tha ruddar. Thia aagnltuda of 
input could affact tha paraiaatara datarainad. Thia affact ia 
dua to a aacond ordar control aurfaca tarn introducad by thia 
float but not pradictad by tha MMLE nathamatical modal. 

It ia auggaatad that furthar work ba dona to avaluata whathar 
thia la tha caaa, and parhapa to Includa control aurfaca float 
into tha mathanatlcal rapraaantation. Tha affact of control aurfaca 
float can ba datamined by varying tha tanalon of tha cabla which 
movaa tha aurfaca. 


_ » • • • aaaaaaaaaa ^ m, ,m m. m, 

Poaaibla problema that could be caaponalbla for tha diffarencea 
in paramatar prediction are listed here: 

- Calibration of transducers. It is suggested that part of the 
error in parameters is due to Inaccuracies in transducer calibration. 

- Uniqueness. It has not yet been determined if the methods 
such as MMLE have a unique solution. The possibility does exist 
of more than one solution to any given maneuver. 

- Control surface float. No attempt was made in this flight 
test program tc ensure a minimum of float of the control surfaces. 

It is suggested that cable tensions be tightened to allowable maximums 
prior to flight testing. 
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7. CJESSMA FLIGHT TEST PROGRAM 


Th« v«rs«tillty of this flight test peckegc was deaonstratad 
In a spin taat progran conductad by Caasna Aircraft. In this 
program tha KU-FRL provided the data management portion of the 
Instrumentation syatam described In this report. Cessna supplied 
the InstrumancHtlon and the airplane. A block diagram of this 
Installation Is shown in Figure 7.1 

7.1 Instrumentation 

The purpose of this program was to Investigate the spin char- 
acteristics of Cessna's latest model 172 airplane. To do this, 
Cessna approached the KU-FRL as to the applicability of the Instru- 
mentation system for this type of test. After Initial evaluation 
It was decided that the measurements described In Table 7.1 would 
be required. It was apparent that the KU-FRL transducer package 
was unable to meet these needs; however, the data management 
portion of the package would be able to. 

The airplane used In this program is shown in Figure 7.2. 

The external modifications to the airplane Include a spin chute 
as well as right-hand and left-hand wing tip booms. 

The spin chute was added for safety reasons. A device for 
deploying thvs chute Is provided to the pilot, allowing him to re- 
trieve the airplane from an unrecoverable spin. Also a release 
mechanism is provided to release the chute after deployment and 
spin recovery. The pilot also wears a parachute in the event the 


91 




CESSNA TRANSDUCERS AND SIGNAL CONDITIONING 


Figure 7.1 Block diagram of Cessna spin test Installation 
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Table 7.1 Cessna Spin Test Measurement Requirements 



SPIN CHUTE 


SYMBOL 

TRANSDUCER 

RA.^E 

6 

e 

ELEVATOR POSITION 

FULL TRAVEL 

6 

a 

AILERON POSITION 

FULL TRAVEL 

6 

r 

RUDDER POSITION 

FULL TRAVEL 

P 

ROLL RATE 

±360 */sec 

q 

PITCH RATE 

±360 •/sec 

r 

YAW RATE 

±360 */sec 

A 

z 

NORMAL ACCELERATION 

-3 to -fS g 

A 

V 

LATERAL ACCELERATION 

g 


ANGLE OF ATTACK LEFT HAND 

-20 to -*-80 • 


ANGLE OF ATTACK RIGHT HAND 

-20 to +80 • 


SIDESLIP ANCLE LEFf HAND 

±45 • 

®R 

SIDESLIP ANGLE RIGHT HAND 

±45 • 

KTASj^ 

TRUE AIRSPEED LEFT HAND 

20 to 180 knots 

KTASj^ 

TRUE AIRSPEED RIGHT HAND 

20 to 180 knots 

H 

P 

PRESSURE ALTITUDE 

0 to 15000 ft 


Figure 7.2 Cessna spin test airplane 
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BLACK AND WHITE PHOTOGRAPH 

spin chute does not deploy or will not release, permitting him to 
leave the airplane in safety. 

Tlie right-hand and left-hand wing tip booms are shown in 
more detail in Figures 7.3. The booms utilize a flow direction 
and airspeed sensor (described in detail /.n Reference 31). This 
sensor allows detet^ilnation of the airspeed, angle of attack, and 
angle of sideslip (as shown in Figures 7.3). A probe is included 
on each wing tip to allow determining th-i true properties of the 
spin. Tlie axis of a spin is generally not at the center of gravity 
of the airplane. (In the C172 it appears to be ahead of the 
center of gravity.) Providing both left-hand and right-hand measure- 
ments allows determining where this spin axis is by using the dif- 
ferences between the measurements from each side. 



Figure 7.3(a) Cessna wing tip booms (supplied by NASA Langley) 
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Figure 7.3(b) Cessna wing tip booms (supplied by NASA Langley) 



Figure 7.3(c) Cessna wing tip booms (supplied by NASA Langley) 
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Inaida tha airplana rockplc tha Caasna inertial rafaranca 
tranaducara (p, q, r, A^) wara aountad on tha aenaor pallat 
aa ahovn In Flgura 7.4. Aa can ba aaan, tha KU-FRL p<;war aupply 
ayatam waa uaad in thia inatallation. Thla waa nacaaaary to provide 
power for tha computer and waa alao utilized to provide powe^- for 
aoma of the tranaducara. Flgura 7.1 ahowa the power aourcaa uaad 
for tha apaclflc dwlcea. 

Tha KU-FRL data management computer la ahown Inatalled In the 
Caaana airplane In Flgura 7.S. 

A chaaa airplane waa uaed In thla flight teat program. Thla 
waa for aafety purpoaea to provide an outside observer who could 
warn the pilot of the spin test airplane (over Che communications 
radio) of any unexpected problems. Also, a video camera was carried 
onboard the chase airplane to record the spin visually. 

7.2 Data Reduction 

Data analysis for this spin program was done by Cessna on 
their Hewlett Packard 9825 microcomputer. Data was transferred 
from the KU-FRL system to the Cessna computer, using the standard 
RS232 ports on each machine (see Appendix A. 9 for Hewlett Packard 
9825 programs). After transfer, the data was plotted on Cessna's 
computer using the program in Appendix A. 9. Figures 7.6 present 
the traces of several of the spins. It can be seen that the data 
recorded produce results capable of analysis. 
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Figure 7.4 Cessna spin test instrumentation installation 



Figure 7.5 Cessna spin test computer installation 
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No. 10 SrU TUNKS » DIUCtlOM 




7.3 lUtulf of Spin Progr>» 


Tht results of the spin progrsm, from ths aspect of this report, 
show the edeptebillty of the KU-FRL-designed date nanageaMiit s/stea. 
The versatility specifically designed into this portion of the system 
allows virtually any 16>channel inat rumen tat ion combination to pro- 
vide the meaaurepMnts. Also shorn is the feasibility of using dif- 
ferent data raduction computers by the use of the standard RS232 port 
for data transfer. 

Mo real problems were encountered by Cessna personnel in using 
the data maiiagement system, even though none of them had any extensive 
microcomputer experience. 
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8. CONCLUSIOWS 


Th« t««e ■ystcfli d«tlgn«d «id «v«luac«d under this program 

has mot tha objactlvas outllnad in Chapter 2. The system 

• Is easy to Install, 

• is virtually self contained « 

• is simple in operation, 

• requires no complex flight maneuvers, 

■ is applicable to general aviation airplanes, 

• is capable of longitudinal and lateral stability analysis, and 
« is low in cost. 

This system has shovm that the technology used is capable of the tasks 
to be performed. 

In the data reduction method all the derivatives contained in 
Equation [*^.5] can be determined. The method also allows determining 
any combination of these derivatives. It must be noted that these 
are the state vector dimensional derivatives which can be converted 
to the normally accepted stability derivatives (as per Reference 22) 
using Tables 5.2 and 5.3. 

Areas have been discovered where further work is required. 

A comprehensive list is Included in Chapter 9. 
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9. RECOMMENDATIONS FOR PHASE III 

Four ar««s hmvo boon suggootod throughout this rsport for in** 
provsmsnt of the KU-FRL Instrumentation system. These are summarized 
here. 


9.1 Equipment 

- Equipment Is required for accurate transducer calibration. 

A pendulum arrangement as per Reference 32 Is suggested as an excel- 
lent means of calibrating the transducers. 

- Size reduction of equipment la suggested. To allow easier 
placement In aircraft, the size of the system could be reduced sig- 
nificantly, especially If the number of packages Is increased to 
form more efficient space utilization. 

9.2 Calibration 

- All transducers should be calibrated as a system. Using the 
actual data acquisition package for transducer calibration Is sug- 
gested as a means to reduce calibration errors. This should be done 
In conjunction with the calibration pendulum of 9.1 above. 

9 . 3 Data Reduction 

- Refinements are required to the current MMLE "BONES" program 
to simplify Its use and add to the versatility. 

- Further study Is suggested to allow Performance analysis 
(l.e., Drag Polar) of the test airplane. Methods similar to those 
of References 4-11, 13, and 18 seem to provide promising solutions. 
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- Some of th« f«atur«« of th« lataat varalon of MMLC (Rafaranca 
33) should ba addad* Spacifically* tha Cramar-Rao bounds addition, 
and Che correction for center of gravity offsets can ba added directly 
into the MMLE program. 

- The addition of the acceleration transformation matrix ([R] of 
Equation [5.4]) to the MMLE program should be explored. 

- Determine the validity of the prediction of a and $ by 
comparing with measured values. 

9.4 Effect of Control Surface Floa t 

- The effect of control surface cable tension should be evalu- 
ated to determine the Influence this has on the parameters predicted. 

9.5 Proof-of-Test Capability 

- Tests are recoiranended In other general aviation airplanes 
to demonstrate the system's adaptability. Recommended are tests 
on a high performance, single-engine retractable, and on a light- 
to-medium, twin-engine airplane. 

- The tests suggested above would also aid in providing further 
insight into the possible "Uniqueness" problem. This should be a 
definite area of research, to validate the ^!MLE (or similar) con- 
cepts. 
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tronics Conference, Kansas City.) 

Development of a Simple, Self-Contained 1981 

Flight Teat Data Acquisition System. 

(Paper presented at SAE Business Air- 
craft Meeting, Wichita, Kansas.) 
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APPENDIX A 


PROGRAMS 

This appendix Includes descriptions, flow cherts, and listings 
of the computer programs required by this flight test system* 

A.l Data Acquisition (AIM-65) 

A. 2 Data Transfer (AIM-65) 

A. 3 Data Receive (MING 11/03) 

A. 4 Engineering Conversion (MING 11/03) 

A. 5 Quick Look Plots (MING 11/03) 

A. 6 Detailed Engineering Conversion (MING 11/03) 

A. 7 MMLE BONES Routines (MING 11/03) 

.1) MMLE Set-Up 
.2) Main MMLE Programs 
.3) MMLE Output Format 
A. 8 Time History Plotting (MING 11/03) 


Cessna Programs 

.1) 

Data 

Acquisition 

.2) 

Data 

Read back 

.3) 

Data 

Receive 

.4) 

Data 

Plotting 
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A.l) DATA ACQUISITIOW PROGRAM 


PRECEDING PAGE BLANK NOT FILMED 


Datcrlptlon: This progrcn, which run* on Ch« AIM 65 » collocta 

and aavaa tha maaaurad atata tlma hlatorlaa. Tha information ia 
collactad and atorad on tha caaaatta tapa in onc-aacond raal-tima 
blocks. Tha data for aach channal la codad as two binary alght-blt 
words totalling slxtaan bits. Tha first word holds tha alght most 
significant bits. Tha sacond word holda tha four moot significant 
bits and tha four laaot significant bits. This glvas a redundancy 
check of tha highest order bits. 

Flowchart : 



* 20 u sac 

baewaao 

sueeassiva 

chaaaals 
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PROGRAM LISTING 


;OATA ACqVISITION 


t 


mforr-o 

LOA 10 

STA UTIL 

BLKCMT-2 

STA KDDM2 

LDA #>C34E 

BUrCNT-4 

STA KDRA2 

STA UTICH 

lUIF-S 

LOA #9C0 

CLI 

OBUr-7 

STA UACR 

LOA 10 

OIT-9 

START 

STA BLXCNT 

iun-$2oo 

LDA #9U 

STA BUCMT41 

IUr2-$300 

STA FORO 

MCI 

n>DKA2-$A481 

LDA fREU 

JSR SWAP 

KDDBB2-$A483 

JSR COMD 

JSR WRITE 

HDRA2-$A480 

LDA fRIU 

REC2 

KDIIB2-8A482 

JSR ooeto 

JSR GKEY 

DBR-89008 

MAIN 

CMF fRECK 

WDC-99009 

JSR GREY 

BNE RECX 

COR-89O0A 

CMF fLOAOK 

IDA BUFCNT 

KDRO- 99008 

BNE MAIN 

CMF #220 

CSR-9900C 

LOA #REU 

BNE REC2 

ESR-9900D 

JSR COMO 

INC BLKCNT 

ISR-8900E 1 

LDA ISLE 

BNE RECl 

MDR1-9900F , 

JSR COMD 

INC BLKCNT-t-1 

V(RT-$C1 

MAIN2 

JMF RECl 

vnH-9C2 

JSR CKEY 

RECX 

OA-$C3 

CMF IRECK 

LDA BUFCNT 

SLB-$C9 

BEQ RECORD 

CNF #220 

REU-9CA 

CMF ICLOSEK 

BNE RECX 

NROY-910 

BNE MAIN2 

LOA #940 

rPT-$04 

CLOSE 

STA UIER 

DA-920 

LOA iwni 

JSR SWAP 

DBRE-940 

JSR COMD 

INC BLKCNT 

CCE-980 

LDA fUTM 

BNE RECXl 

UDRB-9A000 

JSR COMD 

INC BLKCNT41 

UACK-9A00B 

LOX #12 

RECXl 

UIER-9A00E 

CLOSEl 

JSR WRITE 

UT1L-9A004 

LOA #ERA 

JSR ENDREC 

UnCH-9A005 

JSR COMD 

JSR SWAP 

UT2L-9A008 

DEX 

LDA #$FF 

UT2H-9A009 

BNE CLOSEl 

STA BLKCNT 

UIFk-9A00D 

JMF START 

STA BLKCNT41 

B1T5-920 

RECORD 

JSR WRITE 

LOADK-9EF 

LDA #0 

INC RNCNT 

RECK-9BF 

STA IBUF 

BNE RECX2 

CL0SEK-9DF 

STA OBUF 

INC RNCNT41 

TIMElH-927 

LOA #>BUF1 

R2CX2 

TIMElL-910 

STA IBUF41 

JMF MAIN2 


LOA #>BUF2 


*-90400 

STA OBUF<fl 


LDA #$92 

JSR ENDREC 


STA MORO 

LOA #>ttrr 


LOA 11 

STA 9A40S 


STA RNCNT 

LDA #<INT 


LDA #0 

STA 9A404 


STA RNCNT+1 

LDA #SC0 


LDA #$FF 

STA UIER 


STA KDDRA2 

LOA #<9C34E 
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PROGRAM LISTING (continued) 


CKcr 

LOA KDM2 
PNA 

UtA ITIMEIL 
STA UT2L 
LOA fTIHElH 
8TA UT21I 

cmi 

LOA UIFR 
AND lom 
■IQ CKEYl 
PLA 

CMP KORB2 
•NE CKCY 
RTS 


COMO 

PHA 

LOA ESR 
CONDI 
LOA CSR 
AND fNROY 
BNE COMDl 
LOA CSR 
AND IPPT 
BNE COMOl 
PLA 

STA COR 
COMD2 
LOA ISR 
AND ICCE 
BEQ COMD2 
RTS 


ENUREC 
LOT *0 
ENDRl 
LDX 110 
LDA $8000 
JSR WAIT 
ENDR2 
LOA $9001 
JSR WAIT 
DEX 

BNE EN0R2 
LOX #5 
ENDR3 
LDA $8001 
JSR WAIT 
LDA $8002 
STA (IBUF),Y 
INY 

LDA $8003 
STA (IBUF),Y 
INY 
DEX 

BNE ENDR3 
CPY #220 
BNE ENDRl 
RTS 


WAIT 

JSR WAITX 

UAITX 

RTS 


WRITE 
LDA ESR 
LDA #224 
STA WDC 
LOA #WRT 
STA COR 
LOA RNCNT 
JSR WUORO 
LOA RNCirm 
JSR UWORD 
LOA BUCNT 
JSR WUORO 
LOA BLKCNT-t-1 
JSR WUORO 
LDY #0 
URITEl 
LDA {OBUF),Y 
JSR WUORO 
INY 

CPY #220 
BNE WRITEl 
URITE2 
JMP C0MD2 


SWAP 

LOA OBUF-fl 
PHA 

LOA IBUF-l-l 
STA 0BUF41 
PLA 

STA IBUF+l 
LDA #0 
STA BUFCNT 


WUORO 
PHA 
UUOROl 
LOA ISR 
AND #DBRE 
BEQ UUOROl 
PLA 

STA DBR 
RTS 


INT 

PHA 

LOA UDRB 
BPL INTEX 

t;a 

PHA 

LOT BUFCNT 
LOA #11 
STA CNT 
LOA $8000 
JSR WAIT 
ILOOP 
LOA $8002 
STA (IBUF).Y 
LOA $8001 
INY 

LDA $8003 
STA (1BUF),Y 
INY 

DEC CNT 
BNE ILOOP 
STY BUFCNT 
PLA 
TAT 
INTEX 
LDA UTIL 
PLA 
RTI 
END 
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A.2) DATA TRANSFER 


PREaa>INQ PAGE BUNK NOT HLMED 


Description: This program allows the AIM 65 to read tha 

information stored in the DATA ACQUISITION program. The information 
is pasaed to the MINC 11/03 computer using the RS 232 port. 

Flowchart : 



DONE 
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PROGRAM LISTING 


tMTA ucovnv 


RNCNT-0 

STA UTIL 

■UCNT-2 

VRUN-4 

LOA #0 

VBU-S 

STA UTiai 

CNT-A 

LOA #111 

iUPl«l200 

EOR #>rp 

DBR- 19008 

STA SCR 

HOC* >9009 

t 

MAIN 

CDR-S900A 

MSR0-S900B 

JSR CCOM 

CSR-I900C 

CMP #LOADC 

ESR-$900D 

SEQ MAXN2 

ISR»S900E 

JSR INVAL 

MDRl«$900r 

JMP MAIN 

SLP-SC8 

MAIN2 

RDL-$C4 

LOA #$12 

REW-SCA 

STA MDRO 

NRDY-910 

LOA #REW 

TDRE-S02 

JSR COMOA 

CR-SOD 

LOA #SLP 

SCR- 99006 

JSR COMOA 

SDR- $9007 

MAIK3 

L0A0C-$4C 

JSR CCOM 

REA0C-SS2 

CMP IREAOC 

CLOSEC-$43 

BEQ READ 

INAU-9E993 

CMP CLOSEC 

NUMA-$EA46 

BEO CLOSE 

REA0M-$E93C 

JSR i;<VAL 

OUTALL-9E9BC 

OUTPUT-SE97A 

JMP MAIH3 

UTIL-SA004 


UTICN-$A00S 

CLOSE 

UACR-SAOOB 

LDA #REH 
JSR COMDA 

*-$300 

LOA #REU 

CCE 

JSR COMDA 

.BYTE 9B0 

JMP MAIN 

DA 


.BYTE $20 

READ 

MO 

LDY #MRUN-M0 

.BYTE CR.'TAPE ERROR ',$A0 

JSR MESS 

MHUN 

JSR CCNT 

.BYTE CR, 'WHICH RUN NUMBER*. $BF 

CMP #0 

MBLK 

BEQ CLOSE 

.BYTE CR.'HOW MANY BLOCKS ',$BF 

STA VRUN 
READl 

MEND 

JSR RBLK 

.BYTE CR.'LAST BLOw. XIS RU',$CI 

' BCS CLOSE 

MINV 

LOA BLKCNT 

.BYTE CR,' INVALID C0»iAN',$C4 

ORA BLKCNT-Pl 

MERRl 

BIIE REA03 

.BYTE CR.'FILE MARK F0UN*,$C4 

U'Y #MRNCNT-MO 

MR)>;CNT 

.BYTE CR.'RUN NL’MBER'.SAO 

JSR MESS 

MERROR 

LOA RNCNT 

.BYTE CR.'DATA ERROR', $BF 

JSR NUMA 

TEMPO 

READS 

•BYTE $00 

LDA RNCNT 

C»1P VRUN 

*-$400 

BCC READl 

RESETS 

BEQ REA02 

LDA mi 

LDA #REW 

.STA MDRO 

JSR COMDA 

LDA #$C0 

LDA ISLP 

STA UACR 

JSR COMDA 

LDA l$68 ;$68-300BAUD 

JMP READl 

{$34-600 

READ2 

;$IA-1200 

LDY #MBLK-M0 

; $0D-24O0 

JSR MESS 


JSR CCNT 
CKP 10 
•CQ CLOIE 
STA VILK 


SCNDI 


JMP 

SENDIl 

LDA 

RNCNT 

JSR SEND 
LOA RMCNT4>1 

JSR 

SEND 

LOA 

BLKCNT 

JSR 

SEND 

LOA 

BLKCNT41 

JSR 

SEND 

SEND81 

LOX 

#0 

CNVT 

LOA 

BUFl.X 

AND #$F0 

STA 

TEMPO 

INX 


LOA 

BUFl.X 

AND #>P0 

CMP 

TBffO 

SEQ CNVTl 

JSR 

FIX 

OIVTl 

DEX 


LOA 

BUFl.X 

SEC 


■U)R 

BUFl.X 

LSR 

BUFl.X 

AND 

#$3 

CLC 


ROL 

A 

ROL 

A 

ROL 

A 

ROL 

A 

STA 

TEMPO 

INX 


LDA 

BUFl.X 

AND 

#$0F 

ORA 

TEMPO 

ORA 

1940 

STA 

BUFl.X 

INX 


CPX 

#220 

BNE 

CNVT 

LDY 

#0 

SENDB2 

LOA 

BUFl.Y 

JSR 

SEND 


INY 

CPY #220 
BNE SENDB2 
LOA BLKCirr 
CMP #«FF 
BNE SENDB3 

CMP BLKCNT+1 
BEQ END 
SENDB3 
DEC VBLK 
BNE SENDBS 
SENDB4 

LDY #MBU-MO 


{CHANCE TO NOP$ 

:T0 TKANSHIT COUNTI 


{COMPARE HIGH BITS 

{IF BAD FIX HERE 

{CONVERT TO PRINTABLE 
(CHARACTERS 
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JSR 

HESS 

JSR 

GCNT 

CMP 

#0 

BEQ 

CLOSEl 

STA 

VBU 

SENDRS 

JSR 

RBLK 

•CS 

wLOSEl 

JMP 

SENOB 

CLOSEl 

JHP 

CLOSE 

END 


LOY 

#NEN0-N0 

JSR NESS 

JMP 

MAIN) 

HESS 

LOA 

M0,Y 

PHA 


AND 

#$7F 

JSR OUTPUT 

INY 


PU 


BPL 

MESS 

RTS 


GCOM 

JSR 

READM 

JSR 

OUTALL 

RTS 



INVAL 

LOY «^MINV-M0 
JSR HESS 
RTS 

GCNT 
LOA 10 
STA CUT 
GCNTl 
JSR INALL 
JSR DPACK 
see GCHTl 
LOA CNT 
RTS 


OPACK 
CMP #'0' 
see RSPAC 
CMP l$?A 
BCS RSPAC 
AND #$0F 
PHA 

LOA CNT 
ASL A 
ASL A 
CLC 

AOC CNT 
ASL A 
STA CNT 
PLA 
CLC 

AOC CNT 

STA CNT 

CLC 

RTS 

RSPAC 

SEC 

RTS 


SEND 
PHA 
SENDl 
LOA SCR 
AND «TDRE 
•NB SEMDl 
PU 

lOR #irp :CHC roR known char 

STA SOR 

RTS 


RBLK 

LOA CSR 
ANO INROY 
INE RBLK 
LOA #22A 
STA woe 
LOA *ROL 
STA COR 
JSR RWORO 
BCS RBLK2 
STA RNCNT 
JSR RUOKO 
BCS RBLK2 
STA RNCNTH 
JSR RWORO 
BCS RBLK2 
STA BUCNT 
JSR RWORO 
BCS RBLK2 
STA BLKCNT+l 
LOY 10 
RBLKl 
JSR RWORO 
BCS RBLK2 
STA BUFl.Y 
INY 

CPY #220 
BNB RBLKl 
JNP C0M0A2- , 
RBLK2 

JKP C0KDA4 


COMDA 

PHA 

LOA ESR 
COHDAl 
LOA CSR 
AND INROY 
BNE COHOAl 
PU 

STA COR 
COHDAL 
LOA ISR 
AND CCE 
BEQ C0MDA2 
C0M0A4 
LOA CSR 
PHA 

AND 12 

BEQ COMDAS 

LDY #MERRi-HO 

JSR MESS 

PU 

SEC 

RTS 

COMDAS 

PU 


AND #M1 
BNE C0NDA3 
CLC 

RTS 

'C0MDA3 
LOY IMO-NO 
PHA 

JSR HISS 
PU 

JSR MJHA 
LOA ESR 
JSR NUMA 

CLC 

RTS 

RWORO 
LDA ISR 
BIT CCE 
BNE RW0R02 
BIT DA 
BEQ RWORO 
LDA DBR 

CLC 

RTS 

RW0RD2 

SEC 

RTS 


FIX 

LOY #MERROR-MO 
JSR MESS 
RTS 


; INITIAL SET UP 

*-S700 
IDA #$C0 
STA UACR 
LOA #$68 
STA UTIL 
LDA #$00 
STA UTICH 
JMP RESETS 
END 
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A. 3) DATA RECEIVE 


preceding page blank not filmed 


Description: This program accepts raw data from the 
AIM 65. This raw, formatted data is collected in files 
to be used in the ENGINEERING CONVERSION routine. 

Program listing 


0001 

c 

c 

c 

c 

c 

c 

c 

c 

c 


0002 



C 

0003 


0004 


0005 


0004 


0007 

7 

0001 


0009 


0010 


0011 


0012 

998 

0013 


0014 


ooia 


0014 


0017 


0011 


0019 



C 

0020 


0021 

10 

0022 


0023 


0024 


0025 

5 

0024 


0028 



C 

0029 

100 

0030 


0031 


0032 


0033 


0034 

120 

0035 


0034 

no 

0037 


0038 


0039 


0040 

169 

0041 

135 

0042 



c 


PROORSn AININ 

AIN TO HINC fROORAN 
URirrCN IV NARK A NOItCR 


Thlft intuit data fpoa *ha AIN - 45 throush SLU-1 

aa eharactara (22 at a tiaa ) to fill a 400 X 22 chafaetar 
arraw • Uhaii full thla arraw la outautad to tha uaar 
aaaciflad fila. 

OINCNtlON lAODR ( 4 > • lOATAi 400>22 ) • ICHAR< 22 > 


TYPE It 'THIS PROORAN RCAOt 22 
TYPE ar^TlNE AND PtACEl THESE 
TYPE *f‘ • 

TYPE 7 

rORNAT {' WHAT IS THE NANE OF 
CALL ASSION (lif-lfNEW) 

TYPE tt'lST I MILL ATTACH SLU 
ICRR"NTATCH<2) 

TYPE »»S»IERR 
FORNAT <aCRR • '>12) 

TYPE If 'NEXT I WILL SET IT UP 

lADDR(l) > *50010 

IADDR(2> 0 

1A0DR(3> ■ 0 

IADDR(4) ■ 0 

lERR - NTSET(2flADDR<l>) 

TYPE SaStlERR 


CHARACTER WORDS FRON SLUl 
WORDS ON A FILE' 


YOUR OUTPUT FILE T'7) 
1 ' 


FOR READINO' 


TYPE a.'NOW YOU HAUE 2 CHOICES . STOP OR READ IN DATA' 
TYPE a.'l > READ IN DATA ' 

TYPE *f'2 ■ STOP' 

TYPE tf'SHHICH T' 

ACCEPT SiIFCT 
FORNAT (ID 

IF (IFCT .OT. 2) OOTO 10 
GOTO (100»200>pIFCT 


TYPE »f' READ FUNCTION' 

DO no I>1>400 

DO 120 J«lf22 

lERR - NTIN(2.ICHAR(J)f 1) 

XDATAd* J)-ICHAR( J> 

CONTINUE 

TYPE a. 'COUNT - ',I 
CONTINUE 

TYPE tf'WAIT FOR DISK OUTPUT' 

DO 135 K>1>600 

WRITE (in49> (IDATA<K.L>rL-n22) 

FORNAT (22A1) 

CONTINUE 
OOTO 10 


C STOP 
C 


(>043 200 

TYPE 

a. 'STOPPING' 

0044 

TYPE 

a»'DETATCH INPUT PORT 

0045 

lERR 

■ NTDTCH<2) 

0044 

TYPE 

998f lERR 

0047 

STOP 


0040 

END 



AT A' 
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PRECEDING PAGE BLANK NOT FILMED 


0001 

0002 

0003 

0004 

0005 
0004 

0007 

0008 
OOOf 
0010 

0011 

0012 


0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 
0021 
0022 

0023 

0024 

0025 

0026 


0027 

0028 

0029 

0030 

0032 

0033 

0034 

0036 

0037 

0038 

0040 

0041 

0042 

0044 

0045 

0046 

0047 


A. 4) ENGINEERING CONVERSION 


Dftscrlptlon: Tbit proctsa praaantly uaaa two programa 
which run aapcriitaly . The firat routine, AIMCNV, converta 
the raw coded AIM 65 data to voltagea. The AIMCNV program 
makea uae of a macro aaaembly language routine ( CONVRT) 
which per forma the bit manipulations neceaaary to turn the 
AIM 65 coded data into a form uaeable by the MING Fortran 
programa. The aecond routine ,EGRCNV, converta the voltages 
into engineering units. 


Program listing: 


c 

C aXH TO VOLT C0NVER810N 7R00RAH 
C IT MUST 8£ LINKED TO CONVRT TO WORK 
C 

PROORAH AINCNV 

DIHENSION IN(22) >V0LT< 1 1 ) t IV0LT< 1 1 ) 

TYPE AIH FORMAT TO VOLT C0NVER8I0N ' 

TYPE a*' INSERT A DATA DISKETTE' 

5 TYPE «t' Vh«t i« the full nsn* of tho INPUT filoT' 

TYPE * 

CALL A88X8N(l»f-l>RDO) 

10 TYPE «t' Uhat i« tho full nooo of tho OUTPUT filo?' 
TYPE t 

CALL ASSIGN (2f>-liNEU> 

C 

20 TYPE How oonw SECONDS (10 tino oto) do wou wont?' 

ACCEPT tfNUH 
C 

C convort this ooss 
C 

NUH>NUH«10 
DO 100 I-lrNUM 

READ (If 109>END>300) ( IN( J ) f J«1 t 22 ) 

109 F0RHAT(22A1) 

DO 150 K-lfll 

KH-(K«2)-1 

KL>K»2 

CALL CONVRT(IN(KH)f IN(KL)fIVOLT(K) > 
V0LT(K)>IV0LT(K>/16 

150 V0LT(K)-(V0LT(K)«.002435632861)f .06330029324 

WRITE (2fll9> (V0LT(J>f J-lfll) 

119 FORMAT (11F12.9) 

TYPE »f' COUNT- '.I 
100 CONTINUE 

C 

C Hopo It tho toriot of auottiont 
C 

TYPE *f' Mopo in this filot' 

ACCEPT 209fN0RY 

209 FORMAT (Al) 

IF (NORY .EQ. 'Y') GOTO 20 

TYPE *f' Anothop OUTPUT filo fpoo thit INPUT flloT' 

ACCEPT 209fN0RY 

IF (NORY .EQ. 'Y') GOTO 210 

TYPE *f' Anothop INPUT filo in thit OUTPUT filoT' 

ACCEPT 209fN0RY 

IF (NORY .EQ. 'Y') GOTO 310 

TYPE »f' Apo wou donoT' 

ACCEPT 209fN0RY 

IF (NORY .EQ. 'N') GOTO 5 

GOTO 500 

210 CALL CL0SE(2) 

GOTO 10 

300 TYPE «f' END OF INPUT FILE ERROR' 
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310 


004t 

004? 

0050 

0051 
0053 

0053 

0054 

0055 


500 


CALL CLOtKl) 

Tm Vhti i% tht n*Mi INOUr fll* ntMt' 
TYPK t 

CALL AMION<li»-l»RDO> 

•OTO 30 

TY?t I*' OOODBYC' 

•YO? 

fNP 



•TITLE 

CONVRT 




•BBTTL 

CONVERSION ROUT ME 



•QLOBL 

CONVRT 



CONOR Ti 

CLC 





CLR 

R1 




BIC 

•177700iB3<RSl 

IREEL OFF 1ST 3 BITS 



BIC 

•177700»I4<RS> 

lOF THE AROS 



NOUB 

t3(RS)iRl 

IRUT HIOH 4 BITS INTO 

R1 


A8N 

•4iRl 

(SHIFT LEFT 4 



ADD 

i4(R5>fRl 

(ADD IN LOU 4 BITS 



NOV 

R1»B4(R3> 

(RLACE IN RESULT ARO 



ROL 

B4(R5> 




ROL 

B4(RS) 

(ROTATE LEFT 4 



ROL 

t4(R5) 




ROL 

t4(RS) 



TSTi 

CNR 

B4(RS)»4100000 

(SEE IF R08ITIVE OR NEGATIVE 


BNI 

NEQA 



ROSt 

CON 

B4(R5) 

(R08IT1VE NUNBEN8 ARE 

I'S CONRLINENTED 


ADI 

t4(R5)#«l 

(THEN ADD 1 



RTS 

RC 

(DONE 


NEOAt 

NEO 

■4(RS) 

(NS<:iAriVE 4'8 ARE 2'B 

COHRLINiNTED 


RTS 

RC 




.INO 


0001 PROORAN CORCNV 

C 

C Volt to En«ln»«rir>t Unlit cunvvrtlon troaran 
C Iw Nark Ai Nottar Dac> 1?80 
C 


0003 

0003 

0004 

0005 
0004 

0007 

0008 

0009 

0010 


CONNON U0LT8<Kfll)feN0(t0iU) 

10 TYRE Ooltt to Enilnaarlna Unlit convartion' 

TYRE INSERT DATA DISK * 

TYRE »»' Uhat It iha full naaa of tour INRt'T fila?' 
TYRE « 

CALL A8SI0N(l#f-l>RD0) 

30 TYRE «f' Uhai it iha full naaa of wour OUTRUT fllaT' 
TYRE • 

CALL A88I0N(2m>1>NEU> 

C 

C inaui flrti block <T>RdfRt> tconvari t ouiaui 
C 


0011 


READ(1>109> <<VOLTS< J>K) •K>l>U>tJ-l»10) 

0012 

109 

F0RHAT(11F12.9> 

0013 


CALL T0TR8D(T>R8>RD) 

0014 

119 

F0RNAT(3E12.4> 

OOlS 


READ <l»10?)((V0LT8<JfK>tK-lrll)Fj>lfl0) 

0014 


CALL VTOEOR 

0017 

100 

READilf 109fCND>300fERR-1000)((V0LTS( J*K)»Kolfll>f Jalf 10) 

0018 


WRITE (2»129)( <ENQ( JfK)fK«l>ll)t Jal>10) 

0019 

129 

FORNATt 11E12.4) 

0020 


L. ‘ VTOEOR 

0021 


GOTO 100 


C 

C LOOR EXITED BY EOF IN READ 


C 


0022 300 

WRITE 

<2*119> 

TfRS.RD 

0023 

CALL 

TOTRSOiT 

• P8*PD> 

0024 

WRITE 

(3»119) 

TrRStPD 

002S 

CALL 

CLO'SE(l) 


0024 

CALL 

CL08E(2) 
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0027 

002t 

002f 

0030 

0032 


0033 

0034 

0035 


0001 

0002 

0003 

0004 

0005 
0004 
0007 


0008 

OOOf 

0010 

0012 

0013 

0014 

0015 
0014 
0018 

0019 

0020 
0021 
0022 
0024 
0024 

0027 

0028 


0029 

0030 

0031 


0032 

0033 
0035 

0034 

0037 

0038 

0040 

0041 

0042 

0043 

0044 
0043 


TY8C 8t* ANOTHER ritCT ' 
ACCEPT 777fN0RY 
777 PORNATCAl) 

IF< NORY .EO. 'YM OOTO 10 
QOTO 500 
C 

C ERROR 
C 

1000 TYPE ERROR XN READ' 

500 8T0P 

END 


C 

C Subroutin* T0TP8D 

C Thl« tubrouiln* cnnvvrtt tht data to T < tasaaratura) 

C PS (Btatle araatura) 8 PD Idwnaal'e araatura) 

C 

SU8R0UTXNE T0TP8D(TiP8fPD> 

COHPON V0LrS<10tll>«EN0<10fll) 

IX>2 

JJ-1 

T>0 

P8>0 

PD>0 

C 

C Looa to fill TrPSfPD 
C 

11 T«T4M0LT8(n»JJ> 

Jj-jjfl 

IF (JJ .LE. 11) OOTO 21 

JJ>JJ-11 

Xl-Ilfl 

21 PD-PDYVOLTSdlf JJ) 

JJaJJfl 

IF (JJ .LE. 11) OOTO 31 

JJ*JJ-11 

II-II41 

31 PS>P84V0LTS(IIf JJ) 

JJ«JJ93 

IF (JJ .LE. 11) OOTO 11 
If (II .EO. 10) OOTO 101 
JJ«JJ-11 
11-1141 
OOTO 11 
C 

C Avaraaa TiPDfPS 
C 

101 T-T/20 

PD-PO/20 

P8-P8/20 

C 

C Convart to andlnaarlna unlia 
C 

T-( ( ( (280«T)-487)-273. 16)4(9/5) >432 
IF (PD .LE. 4.524) OOTO 206 
PD-64S8.47-(1436.6iaPD) 

OOTO 214 

204 PD-264B.5S-(550.4aPD) 

216 IF (PS .LE. 2.303) OOTO 226 

PS- ( 12919. 894aP8)- 19754. 521 
GOTO 234 

224 P8-(9451 .794aPS)-11747.49 

234 CONTINUE 

RETURN 
END 


c 

c 

c 

c 

c 

c 

0001 

0002 

0003 12 

0004 

0005 
0004 
0007 
0000 

0009 

0010 
0011 
0012 

0013 

0014 

0015 102 
0014 
0017 


tubpoutln* OTOCOR 

This tubroutlnt eonvvrtt tht <S»i» from volt* to onotnoopina 
unit** It than poopOop* thoo a* rollou*i 
that* •*tAztAKfOaltaE»PhitPfAwrPfOaItaA«daltaR 

tUlROUTINC VTOCOR 

COHRON VOLrt(lO»ll)>ENO<10>ll) 

00 102 XA-lilO 

CN0<IA»UXV0LTt<IAt9>«. 211 1)4.00738 

CNO<IA>2)«><CVOLTt(IAr4)4.0Ol)/.0fft8)/S7.3 

EN0<IAf3)a< <VOLTS<IA»S>4.002>/1.001>-1.0 

ENO<lA>4)«<VOLT8<IAt3)-.0O2)/2.499 

EN0(IA»9)a<V0LT8<XA»10)«4.2442.4)/S7.3 

ENO<XAf4)><OOLT8<XA>2>«.3117)4.01342 

CNO< XAf 7)>-( (MOLTS (XAi8>4.003S)/.09933>/97.3 

EN0(XA>8)>V0LTS<XAf4)/10.027 

CN0(XA*9>— ((V0LT8(XA»7)-.011)/.10134)/97.3 

ENG(XA»10)-(-<V0LT8(XA»ll)t4.73091>4.99940)/97.3 

ENO(XA> 11)«( (M0LT8( XA> 1 > «7 . 92849 ) 48 . 20487 ) /97 . 3 

CONTINUE 

RETURN 

END 


1 


0001 

0002 

0003 

0004 

0005 

0004 

0007 

OOOB 

0009 

0010 
0011 
0012 

0013 

0014 

0015 
0014 

0017 

0018 

0019 

0020 
0021 


0022 

0023 

0024 

0025 
0024 
0027 


0028 

0029 

0030 

0031 


0032 


0033 

0034 

0035 
0034 


A. 5) QUICK LOOK PLOTS 

Description: The QUICK LOOK PLOT program Is used as an 

aid In choosing appropriate flight data for further analysis. 
The routine collects the engineering units data and plots It 
on a graphics CRT terminal. 

Program listing: 


c 

c 

c 

c 

c 

c 

c 

c 


c 


1 

c 

c 

5 

25 


24 

4 

7 

C 

C 

C 


99 


21 

C 

C 

c 

223 


22 

24 

C 

C 

c 

c 


« 

« NAIN PROORAH FOR RUOTTINO NHI.E DATA 

« 

t use FOLtOUINO SUBROUTINE! TO LINK : 
« XNlT.FL0T55>F0UR»0RIDt0RAFH. 


C0HN0N/8TATU8/ISTAT(14) 

0XNEN8I0N XARRAY(512) . HH( U > 1 ) t DATA( 240 » 1 1 ) 
COHMON /F:VE/ KlfK2.K3>K4.N 
COHNON/FIOEA/ DATA 
BYTE YE8fN0>ANS«NANE<15) 

DATA YE8.no /'Y'.'MV 
DATA XSTAT/1440/ 

CONTINUE 


CALL XNXT 
TYPE 5 

FORNATi' TYPE IN NANE OF FILE UXTH MEASURED DATA'/) 
format: 14A1 ) 

ACCEPT 25. (NAME<X>.I«1.14> 

QPEN<UNXT>2. NAME-NAME. TYPE- 'OLD' .ACCESS- 'SEQUENTIAL 
1 READONLY. FORM-'UNFORMATTED' ) 

DO 24 1-1.14 

name: I) -' ' 

CONTINUE 
TYPE 4 

format:' type in number of time points'/) 

REA0:S.7> N 

fornat:i3> 


READ IN DATA FROM FILE 
DO 21 I-l.N 

read:2) :hh:k.i>.k-i.ii) 

CONTINUE 
DO 21 J-1.11 
data:i .j)-hh: j.i) 
CONTINUE 


CONTINUE 
CALL INXT 
F0RMAT:11E12.4) 

format:ai) 


CALI 8RID<500.4S) 


. 
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C 

c 

c 


CONTINUE 
CALL four:i) 
CALL four: 2) 
CALL four: 3) 
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0037 

003t 

003f 

0040 

0042 

0043 

0044 

0045 
0044 


0040 

0049 

0050 

0051 


0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 


0001 


0002 

0003 

0004 

0005 
0004 

0007 

0008 


0009 

0011 

0013 

0014 

0015 
0014 

0017 

0018 

0019 

0020 
0021 
0022 

0023 

0024 

0025 
0024 

0027 

0028 

0029 

0030 

0031 

0032 

0033 

0034 


rm 51 

51 rORNATl' DO YOU WANT TO TAKE ANOTHER 1.00K AT THE DATAfCY/N)') 
REA0<5»24> AN8 

IF<AN8.EQ.YE8) 00 TO 223 
TYRE 52 

52 FORNATf' DO YOU WANT TO LOOK AT ANOTHER DATA FILET<Y/N)'> 
READ(5»24) ANI 

CALL CL08E(2> 
ir(AN8.E0.YE8> 00 TO 1 
C 

c 

c 

CALL RL0T55( 2>512> 142f 4432444 tlSTAT) 

CALL RL0T55(0>>li0»X8TAT) 

RETURN 

END 


SU8R0UTINE TNIT 
C0NN0N/8TATU8/I8TAr(14) 

DATA ISTAT/UtO/ 

CALL RL0T55(13»72>>I8TAT) 
CALL PL0T55(13»74.tISTAT) 
CALL PLQT53<2rI4Sl2**ISTAr> 
RETURN 
END 


SU8R0UTINE FOUR<NZ) 

C 

C tttttttttttttttt********* 

C « THIS SUDRQUrtNE WILL DISPLAY FOUR VARIABLES ON 
C « THE CRT EVERY TINE IT IS CALLED FRON THE MAIN 
C « PROORAN, 

C 

C0MN0N/8TATU8/ISTAT ( 14) 

INTEGER IARRAY(512> 

DINEN8I0N DATA(240ftl) 

COHHON /FIVE/ K1 rK2>K3tK4t N 
INTEGER GAINiil) 

CONHON/FIVEA/ DATA 

DATA GA1N/I50r90.39>157f 225 f 75 >90 f 157 >90*225 f 225/ 

C 

C 

IF<NZ.EQ.2.) GO TO 52 
IF(NZ.EQ.3.> GO TO S3 
C 

K1 * 3 
K2 - 2 
K3 > 1 
K4 ■ 5 
GO TO 54 

52 CONTINUE 
K1 ■ 9 
K2 - 11 
K3 ■> 4 
K4 ■ 5 
60 TO 54 

53 CONTINUE 
K1 - 8 
K2 ■ 7 
K3 ■ 4 
K4 - 10 

54 CONTINUE 
C 

DO 41 K>lrN 

IARRAY<2«K) • DATA <K>K1 )*. 5729S«GA1N< K1 ) 445 
IARRAY(2«K>1) •DATA(K>K2)».S7295«6AIN(K2)490 
41 CONTINUE 

CALL 0RAPH(2«N>1ARRAY) 


126 


0033 

0034 
0037 
0031 
003f 


0040 

0041 

0042 
0044 

0044 
0048 
0048 
0050 

0031 

0032 

0033 

0034 
0053 
0034 

0037 

0038 

0039 

0040 
0061 
0062 

0063 

0064 

0045 
0066 

0067 

0068 

0069 

0070 

0071 

0072 

0073 

0074 
0073 

0076 

0077 

0078 

0079 

0080 
0081 
0082 
0083 


0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 
0011 
0012 


0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 
0011 


C 

DO 42 K"1>N 

I4RRAr(2«K> • DATA<K»K3>«.37293*Q4IN(K3)>133 
IARRAY(2*K-1 > - DATA(KtK4)i.3729S*0AIN<K4>tl80 
42 CONTINUE 

CALL 0RARN(2«NrXARRAr> 

C 

C 

C 

CALL rL0T33(9f23tl.ISTAT> 

CALL RL0T55(12f QUICK LOOK DATA RLOTS *««'»ISTAT) 

C 

IFlNZ.CQ.l.) 00 TO 31 
XF(NZ.C0.2.) 00 TO 32 
ir<NZ.EQ.3.> 00 TO 33 

31 CONTINUE 

CALL RL0T55(9>5Ot4»X8TAT) 

CALL RL0T53(12»> 'ELEVATOR ROSN. 20 DEO.'flSTAT) 

CALL PL0T3S(9rS0f8>I8TAT) 

CALL PL0T35(12>f 'PITCH ATTITUDE* 30 DE0.'*I8TAT) 

CALL PL0TS5(9tS0*13U8TAT) 

CALL PL0T55<12i»'PITCH RATE* 50 DEO/SEC .'* I8TAT ) 

CALL PL0T35(9»50tl7»I8TAT) 

CALL PL0T55(12f p'NORHAL ACCEL.* 2 0.'*I3TAT) 

C 

00 TO 34 

32 CONTINUE 

CALL PL0T55(9*S0*4*ISTAT) 

CALL PL0T55< 12* * 'ELEVATOR POSN.* 20 DEG.'*ISTAT) 

CALL PL0T55(9*!>0*8*ISTAT) 

CALL PLOT53(12**'LONOITUDINAL ACCEL.* .3 0'*I8TAT) 

CALL PL0T55(9*50tl3*ISTAT> 

CALL PL0TS3( 12** 'RUDDER POSN.* 20 DEO. '*1STAT> 

CALL PL0T55(9*50*17*ISTAT> 

CALL PL0T33(12** 'TAM RATE* 50 DEO/SEC •'* ISTAT > 

00 TO 34 

33 CONTINUE 

CALL PL0T55(9*S0*4*XSTAT) 

CALL PL0T55( 12** 'AILERON POSN.* 20 DEO.'* ISTAT) 

CALL PL0T55<9*50.8*ISTAT) 

CALL PL0T5S( 12* * 'BANK ANGLE* 60 DEO. '*ISTAT> 

CALL PL0T55<9*S0*13*ISTAT) 

CALL PL0T55<12**'R0LL RATE* 50 DEG/SEC .'* ISTAT > 

CALL PL0T55(9*50.17*ISTAT) 

CALL PL0T53< 12** 'LATERAL ACCEL.* .5 0.'*ISTAT) 

34 CONTINUE 

CALL PL0rS5(9*0*20*ISTAT) 

DO 1 J>1*S000 
DANNY> C0S(4S> 

1 CONTINUE- 

RETURN 
END 


SUBROUTINE ORID( IDX* XDY > 
C0NM0N/8TATU8/ISTAT(16) 

CALL PL0TSS<2*14S12**ISTAT) 
CALL PL0TSS(9*0*0*ISTAT> 

CALL PL0T33<10***IS7AT> 

CALL PL0T55<2*14-32-l-64* * ISTAT) 
DO 3 I-1*512*IDX 

3 CALL PL0T55<5*I-1*1*ISTAT) 

DO 4 I-1*236*IDY 

4 CALL PL0TS5(4*l*I-lf ISTAT) 

RETURN 

END 


SUBROUTINE ORAPHTN* lARRAY) 

COHNON/STATUS/ISTAT( 16) 

DIMENSION XARRAY(S12) 

NUNBER-ISTAT(8)/8 
CALL PL0TSS(7*O*O*ISTAT) 

CALL PL0TSS<8*512*0*ISTAT) 

CALL PL0TS3< 2 * U (NUMBER'*-! ) «2 * (NUMBER-)'! ) «10 * ISTAT ) 
CALL PL0T55(3*-N*IARRAY*ISTAT) 

CALL PL0T55(1*1-NUMB£R**ISTAT) 

CALL PLOT33(9*10*1*XSTAT) 

END 
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A.6) DETAILED ENGINEERING CONVERSION 

Description: The CRINST program performs Che detailed 

corrections for instrument offsets from Che body axes. Biases 
on the accelerometers are also removed in the corrections. 

Program listing: 


0001 


0002 

0003 

0004 

0005 

0004 

0007 

0005 

0009 

0010 
0011 
0012 
0013 


0014 

0015 
0014 

0017 

0018 

0019 

0020 


0021 

0022 

0023 

0024 

0025 
0024 

0027 

0028 

0029 

0030 

0031 

0032 

0033 

0034 

0035 
0034 

0037 

0038 


PR08RAN CRINST 

C.... RROORAH TO HODIFY THC RAN CNOXNCERINO 
C.... DATA FOR INSTRUMENT CORRECTIONS 
C.... THAI MEASURED FROM lODY TO INST AXES 
C.... XSARi YSAR> AND ZSAR FROM SODY TO INST AXES 
BYTE NAME<15> 

DIMENSION FI(ll)>FIMl(ll)>FIFl(ll>fDATA(400f 11) 

DATA DORfO /57.2957Sf 32. 174/ 

DATA THAI Z-5.494/ 

DATA XSAR>YSARiZiAR /40. 052*41 . 179»41 .430/ 

C.... TRANSDUCER F08TI0N8 RECALCULATED ON 4-FEB-81 

4 FORMAT (14A> 

5 F0RNAT(I10tSF12.4) 

C.... SET LAST BYTE OF CHARACTER STRING TO NULL 
NANE(lS)-0 

C.... ENTER THC FILE NAME FOR THE DATA TO BE CORRECTED 

10 FORHATi' ENTER THE FILE NAME FOR THE RAW ENOINEERINO DATA'*/) 
TYPE 10 

ACCEPT 4*(NAME<I)*I>1*14) 

OPEN < UNI T>1 * NAME-NAME * TYPE- ' OLD ' * ACCESS- ' SEQUENTIAL ' * 

. READ0NLY*F0RN-'F0RHATTED'*REC0R0SIZE-132) 

C.... ENTER THE FILE NAME FOR THE DATA TO BE 8AUED ON 
C.... CLEAR OUT OLD FILENAME 
DO' 22 1-1*14 
NAME(I)-' ' 

22 CONTINUE 

20 FORMAT (' ENTER THE FILE NAME TO HOLD THE CONVERTED DATA'*/) 
TYPE 20 

ACCEPT 4*(NAME(I)*I-1*14) 

OPEN ( UN I T-2 * NAME-NAME * T YPE- ' NEW ' * ACCESS- ' SEQUENTIAL ' * 

. FORM-'UNFORMATTEO'* 

. BUFFERCOUNT-2) 

C.... READ DATA FROM HP DATA FILE 

28 FORMAT(' ENTER THE NUMBER OF TIME POINTS TO BE CORRECTED'*/) 
TYPE 28 

29 FORMAT (I 10) 

ACCEPT 29*IEND 

30 F0RNAT(11E12.4) 

DO 31 1-1*IEND 

READ ( 1 * 30 ) ( DAT A ( I * J ) * J-1 * 1 1 ) 

31 CONTINUE 
CLOSE(UNIT-l) 

C.... TRANSFER FIRST TWO DATA POINTS 
DO 35 1-1*11 
FZN1(I)-DATA(1*I) 

35 CONTINUE 

C.... CORRECT FOR SION ERRORS IN CALIBRATIONS 
FINK 4)— FIMK 4) 

FINK 4)— FIMK 4) 

FINK 8)— FINK 8) 

FINK 9)— FIMK 9) 

FINKID— FIMKll) 

C... CORRECT FOR GYRO MISALINEMENT 
C08THI-C0S(THAI/DGR) 


0039 

0040 

0041 

0042 

0043 

0044 

0045 


8INTHI-SIN(THA1/DGR) 

FIMKD-FINKD-THAI/DGR 

PI- FINK7)«C08THI4FIMK9)«8INTH1 

RI— FINK7)«SINTHI4FIMK9)«C08THI 

FIMK7)-PI 

FINK9)-RI 

C.... WRITE THESE VALUES TO THE OUTPUT FILE 
URITE(2XFIHK I) *1-1*11) 
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i 


C.... TMNtrCR NEXT OAT* POINT 
0044 ‘ 00 34 

0047 PKX) ■DATA<2»1) 

0040 34 CONTINUE 

C..*. CORRECT FOR IION ERROR! IN CALIORATIONI 
004f FI< 4) •-FK 4) 

0050 FI( 4) — FI( 4> 

0051 FI( •) ■-FK t) 

0052 FI( f) —FK t) 

0053 FKll) —FKll) 

C.... ADD CORRECTION FOR OYRO NltALINENENT 

0054 FI(1)«FIU)«THAI/D0R 

0055 PI- FI<7>*C00THHFI(9)«IINTHI 

0054 RI»FI(7)ttINTHI4FI(9)«C00TNZ 

0057 FK7)-PI 

0051 FI(V)-RI 

C.... PRINT OUT AVERAGE VALUEI FOR THIl AND PHIl 
0059 THAl-<<FINl<l)9FI<l))/2.)tD0R 

0040 PHIl-<(FINl<4)4FI<4)>/2. )«00R 

0041 40 FORNAT(' THAI (IN DEO) - '>F12.4»' PHIl (IN DEO) ■ '»F12.4> 

0042 TYPE 40>THA1»PHI1 

0043 KOUNT-2 

0044 COtTHA«COt(TMAl/DOR) 

0045 IINTHA-IIN( THAI/DOR) 

0044 COtPHI-COt(PHIl/DOR) 

0047 IINPHI-IIN(PHI1/D0R) 

C.... START CORRECTION LOOP 
0041 50 C'^NTINUE 

0049 kOUNT-KOUNT41 

C.... IF KOUNT It GREATER THAN lEND GO TO 1000 
0070 1F(K0UNT.QT.IEND>00 TO 1000 

TRANFER NEXT TINE POINT 

0072 00 37 X-l>lt 

0073 FIP1(I)-DATA(K0UNT>I) 

0074 37 CONTINUE 

C.... CORRECT FOR SION ERRORS IN CALIBRATION 

0075 FIPK 4>— FIPK 4) 

0074 FIPK 4>— FIPK 4> 

0077 FIPK 8>— FIPK B> 

0078 FIPK 9>— FIPK 9) 

0079 FIPKIU— FIPKll) 

C.... CORRECT FOR NISALINEHENT ANGLE 

0080 FIPK1)-FIPK1)>THAI/DQR 

OOSl PI> FIPK7)tC08THK>FIPK9)*8INTHI 

0082 RI»-FIPK7)*8INTMKFXPlf9)tC0STHI 

0083 FIPK7)-PI 

0084 FIPK9)-RI 

C.... GET AN BACK TO ORIGINAL SIGNAL 

0085 FK3) ■FI(3)41.00 

C.... CONFUTE PDOTt QOOTt AND ROOT 
0084 PDLT -(FIPK7)-FINK7))/0.2 

0087 QOOT -(FIPl(2)-FINK2))/0.2 

0088 ROOT ■(FIPK9)-FINK9) )/0.2 

C.... CORRECT ACCELERATIONS FOR OFFSET 

0089 P -FI (7) 

0090 G >FI(2) 

0091 R •FI(9> 

0092 AXPl ■FI(4)«C0STHI 

0093 AXP2 •FK3)*8ZNTHI 

0094 AXP3 •(R«R4a«Q)«XBAR/G 

0095 AXP4 ■(P«Q-RD0T)«YBAR/0 

0094 AXP5 -(P«R40D0T)«ZBAR/0 

0097 AXFIX -AXP1-AXP24AXP3-AXP4-AXP5 

C098 AX -AXFIX 

0099 AY -FK8) -(P«QfROOT)«XBAR/Gf (P«P4R»R)»YBAR/G 

, >(Q»R-PD0T)«ZBAR/0 

0100 AZ — FK4)«81NTHX-FK3>«C08THI-(P«R-QDOT>»XBAR/0 

. -( Q«R<fPDOT ) «YBAR/6-f ( P«PtQ«Q > XZBAR/G 

0101 FK3) — AZ 

0102 FK4> - AX 

0103 FK8) « AY 
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0104 

0105 
0104 


c • • • • 


CORRECT FOR ACCCLCRONCTCR RIAI 
FX<3> •FI(3)-C0ITHA«C0IPHI 
FX<4> >FX<4)-tXNTHA 
FX<t> ■FX<t)4C0«THA«SXNFNX 
C.... WRXTC VAtUCt ON OUTPUT 
C.*.. TYPE OUT KOUNTER 


0107 


TYPE StKOUNT 

OlOt 

c< • • • 

URtTE<7.MFX(X>»X-lill) 
tUCKET IRXOADE VALUES THRU TXHE 

OlOf 


DO 100 X"l«ll 

0110 


FXNim-FX<X> 

0111 


FX<I) >FXP1(X> 

0112 

100 

CONTI, iUE 

0113 


00 TO SO 

0114 

1000 

CONTINUE 


c • « • • 

TRANSFER LAST DATA POINT 

Oils 

c • • • « 

WRXTC<2)(FX<X>>I«1>11) 
CLOSE DATA FILE 

0114 


STOP 

0117 


END 
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0001 


0002 

0003 

0004 

0005 
0004 

0007 

0008 

0009 

0010 
0011 
0012 

0013 

0014 

0015 
0014 

0017 

0018 


0019 

0020 
0021 
0022 

0023 

0024 

0025 
0024 

0027 

0028 

0029 

0030 

0031 

0032 

0034 

0034 
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A.7) MMLE BONES ROUTINES 

This appsndlx describes the MMLE programs. Tiie first 
program required is the one that sets up the input matrices, 
as well as defining for the MMLE program which parameters 
it is to estimate. The MMLE programs, as well as their 
output format is also presented. 

A. 7.1) MMLE SETUP 

Description: The setup program is an interactive 

program which sets up the input data for the MMLE BONES 
routine. Non-dimensional derivatives, geometric, and inertia 
data are input and used to form the initial estimate to 
the MMLE program. 

Program listing: 


peOORAH lETUR 

C.... THIS RROORAH SETS UP THE DATA USED IN SONES NNLE. 

C.... DIHEN8I0NAL DERIVATIVES ARE 8UILT UP FROH NON-DIHENSIONAL 
C.... INPUT DATA AND AIRPLANE OEOHETRIC DATA. 

C.... DEFAULT VALUES (IF THEY EXIST) ARE SHOWN AFTER EACH CUESTION. 
DIMENSION A(Si4>»S(5t4)»AA(S»4>fBI(St4)fAP(ef4)»iP<8i3> 
DIHEN8I0N ZER0(4)>DIA8(4)fDl<7f7) 

DOUiLE PRECISION CASEtTEHP 
BYTE BANNER (4» 80) 

DATA VALUE »IVALUEfAA>BBrAPf BP/0. tO >20*0. »30«0.f32«0. >24*0./ 
DATA Bl>BIA8*ZER0/49«0.r4«V.»4«O./ 

DATA CA8E.TEHP/' '/ 

DATA BANNER/320«' V 
C.... SET DEFAULT VALUES 


NN 

-200 

ITR 

-10 

HZ 

-7 

NAPR 

-0 

HH 

-0. 10 

EPS 

-0.0 

TIME 

-0.0 

ALPHA 

-0.0 

XLA 

-1.0 


C.... UNIT 1 WILL BE THE FILE NUMBER OF THE FILE FOR 
C.... THE DATA DISK WHICH IS ASSUMED ON DYi: 

C.... OPEN UNIT 1 

2 F0RMAT<80A1) 

9 F0RMAT(' ENTER A BANNER OF UP TO FOUR LINES.') 

TYPE 9 

1 format:' enter line: 'fii) 

DO 3 I>1>4 
TYPE 1*1 

ACCEPT 2*(BANNER(I*J)*J>1*80) 

3 CONTINUE 
10 CONTINUE 

30 format:' ENTER ’LONG* OR *LATR* FOR THE TYPE OF CASE'*/* 

. ' TO BE SET UP.') 

TYPE 30 
40 F0RHAT(1A8) 

ACCEPT 40*CASE 

IF ( CASE. EQ. 'LONG' )QPEN (UNIT-1 *.\'AME> ' DYI : HMLELO. DAT '* TYPE-' NEW ' * 
. RECORDSIZE-94*lNlTIALSlZE-SO.riSPOSE-'SAVC' ) 

IF (CASE. EG. 'LATR')OPEN (UNIT-1 . NAME- ' OYl : MMLELD. OAT '. TYPE- 'NEW' * 
. REC0RDS1ZE-94*INITI4LS1ZE-50*DI8PG?E-'SAVE') 

IF((CASe.NE. 'LONG'). AND. (CASE. NE.'LATR'>)00 TO 10 
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oosa 


003f 

0040 

0041 

0042 

0043 

0049 

0044 

0047 

0040 

0049 

0091 

0092 

0093 

0094 
0099 
0097 
0099 

0099 

0040 

0041 

0042 

0044 

0049 

0044 

0047 

0040 

0070 

0071 

0072 

0073 

0074 
0074 
0077 
0070 
0079 
0000 
0002 

0003 

0004 
0009 
0004 
0000 
0009 

0090 

0091 

0092 
0094 
0099 
0094 
0097 
0090 

0099 

0100 

0101 

0102 

0103 

0104 
0109 
0X04 
Ui07 
0100 

0109 

0110 


C.... CRROO TOM ir OI090NCC 10 NOT *tONO* OR *L*TR* 

OMIC 0*T* FOR IITHCR LONOITUOINRL OR LRTERAfOlRICTXONAL CAtl 
90 rORNAT<' INTER TNE NUNOER OF OATA FOlNTt TO OE FROCEOOED. ' f /> 

' (BEFAULT It 200)') 

TTFE 90 

40 F0RNAT(2F10.0) 

41 F0RNAT(1I9) 

ACCEFT 41*IVALUE 
XFdOALUEtOT.OtNNalVALUE 

IVALUE>0 

70 FORMAT (' ENTER THE NUMBER OF ITERATXONO TO BE FERFORMEB. ' */* 

. ' <BEFAULT 10 10)') 

TTFE 70 

ACCEFT 4 I 1 IVALUE 

XF<X9ALUE.0T.0)XTR>iIVAI.UE 

XVALUE«0 

00 FORMAT(' ENTER THE NUMBER OF OBOERVATXONO. ' 

. ' (BEFAULT XO 7>') 

TTFE 00 

ACCEFT 41»XUAI.UE 

XF(XVALUE.OT.O)MZ«XVALUE 

XVALUE>0 

90 FORMAT(' ENTER THE CONTROL NUMBER FOR THE AFRORX OFTXON.'f/* 

. ' (BEFAULT XO 01 WHICH XO NO AFRORX VALUED') 

TTFE 90 

ACCEFT 41fXVALUE 
NAFR >0 

IF ( IVALUE . NE . 0 >HAFR>XVALUE 
XVALUE-0 

100 FORMAT(' ENTER THE DELTA TIME INCREMENT .'»/ > 

. ' (DEFAULT XO OdO)') 

TTFE 100 
ACCEFT 40* VALUE 
XF(VALUE.OT.O. IHH-VALUE 
VALUE >0. 

110 FORHAT(' ENTER THE VALUE FOR EFS.'*/*' (DEFAULT It 0.0)') 

TTFE 110 
ACCEFT 40*VALUE 
IF(VALUE.OT.O. >CFS-VALUE 
VALUE >0. 

120 FORHAT(' ENTER THE VALUE FOR TIME.'*/*' (DEFAULT XO 0.0)') 

TTFE 120 
ACCEFT 40* VALUE 
XF(VALUE.OT.O. )TIME>VALUE 
VALUE >0. 

130 FORMAT(' ENTER THE VALUE FOR ALPHA.'*/*' (DEFAULT IS 0.0)') 
TTFE 130 
ACCEFT 40*VALUE 
I F ( VALUE . OT . 0 . > ALFHA- VALUE 
VALUE -0. 

140 FORMAT!' ENTER THE VALUE FOR XLA.'*/*' (DEFAULT IS 1.0)') 

TYPE 140 

ACCEFT 60*VALUE 

IF ( VALUE . OT . 0 . > XLA«VALUE 

DO 141 I>1*4 

HRITE(1.2>(BANNER(I*J)*J>1*80) 

141 CONTINUE 

UR ITE ( 1 * ISO ) Nt^ * I TR * MZ * MAFR 
URITE(1*140>HH*EP8* TIME .ALPHA* XLA 
190 F0RHAT(7I10> 

140 F0RHAT(0F10.4> 

C.... ENTER THE HASS AND GEOMETRIC DATA 

170 FORMAT!' ENTER THE AIRPLANE WEIGHT. (IN LBS)') 

TYPE 170 

ACCEFT 4O*UEI0HT 
AMS8 >UEX0HT/32.174 

ISO FORMAT!' ENTER THE AIRPLANE WING AREA. (IN FT««2)'> 

TYPE 180 
ACCEPT 40*8 

190 FORMAT!' ENTER THE AIRPLANE CBAR. (IN FT)') 

TYPE 190 
ACCEPT 40* CBAR 
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<1N fT>') 


0111 

0112 

0113 

0114 

0115 
OIU 

0117 

Oil! 

0120 

0121 

0122 

0123 

0124 

0125 
0124 

0127 

0121 

0120 

0130 

0131 

0132 

0133 

0134 

0135 
0134 
0137 
013t 

013V 

0140 

0141 

0142 

0143 
0145 

0144 

0147 
0140 
OMV 

0150 

0151 

0152 

0153 

0154 

0155 
0154 
0157 
0150 
015V 

0140 

0141 

0162 

0164 

0166 

0167 

0148 

0169 

0170 

0171 

0172 

0173 

0174 

0175 

0176 

0177 
0170 


1V5 rOiMATr CNTIR TNC UXNO tPAN. (IN fTI') 

rrvt 1 V 5 

ACCEVT 40»0VAN 

200 VOAHATi' ENTER TNE ALTITUDE OF THE FLIONT/RUN. 

TYRE 200 
ACCERT 40 *H 

C.... COHRUTE ATHOIRHERIC CONOITIONI FRON ARRROXIHATE RELATIONE 
TA >51t.7*H40.0035i 

IF<TA.LT JV0.)TA«3V0. 

RA >2114. 33»(l.-0.000004t7E4«N>*»5. 2532 

RNO >RA/U714.54«TA> 

AUEL >4V.02«tQRT(TA) 

C.... ENTER THE ETEADY-fTATC FLIGHT CONOITIONI 

210 FORNAT(' ENTER THE ITEADY ITATE VELOCITY. (IN FT/IEC>'> 

TYRE 210 
ACCERT 40»U1 

CLl >2.«UEXOHT/(RHOtUl«Ul«l> 

C.... AtSUNE L/D OF 10. 

CDl >CL1/10. 

C.... ABIUHE CXT1>C01 
CXTl >CD1 

220 FORHAT(' ENTER THE ITEADY ITATE THETA. (IN DCG)'f/f 
. ' (DEFAULT II 0.0) '> 

DOR >57.2V57B 
TYRE 220 
ACCERT 40tTHA 
THA >THA/OOR 

230 FORNAT(' ENTER THE STEADY ITATE . '^NK ANGLE. (IN DEO)'»/i 
. ' (DEFAULT IB 0.0)') 

TYRE 230 

ACCERT 40rRHX 
RNI >RHI/OOR 

240 FORHATI' ENTER TNE STEADY ITATE ANGLE OF ATTACK. (IN DEO)'f/. 
. ' (DEFAULT IB THETA)') 

TYRE 240 

VALUE >0. 

ALR >THA«OOR 
ACCERT 40*VALUE 
IFXVALUE.NE.O. >ALR«VALUE 
ALR xALR/OOR 
8 INALR>BIN(ALR) 

COBALR-COKALR) 

IINTHA>IIN(THA) 

COBTHA>COI(THA) 

SINRHI> 8 IN(RNI) 

COIRHI-COS(RHI) 

TANTHA>BINTHA/COITHA 
C.... ENTER THE INERTIAL DATA 

260 FORHAT(' ENTER lYYB. (IN SLU0*FT««2 > ' ) 

TYPE 240 
ACCEPT 40fAIY 

270 FORNAT(' ENTER IXXB. (IN ILU04FTM2 ) ' ) 

TYRE 270 
ACCERT 40fAXX 

280 FORNAT(' ENTER IZZB. (IN SLU0«FT*«2 ) ' > 

TYRE 280 
ACCEPT 6 O 1 AIZ 
C.... SPLIT FOR CASES 

IF(CAIE.EQ. 'LONG' )00 TO 300 
IF(CASE.EQ. 'LmTR' )00 TO 500 
STOP 

C.... LONGITUDINAL CASE 
300 CONTINUE 

310 FORHAT(' ENTER CDU. 0 OR 1 . ') 

TYRE 310 

311 FORHAT(' ( 1 IF THIS IS A VARIABLE! 0 OTHFT^UISE )') 

TYRE 311 

ACCERT 60rCDU>AA(2.2) 

320 FORHAK' ENTER CXTU.') 

TYRE 320 
ACCERT 60.CXTU 

330 FORNAT(' ENTER COA> 0 OR 1.') 

TYRE 330 

ACCERT 60»CDAtAA(2f3) 
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Ol7f 

340 

rOKHATt' CNTfR CDOIi 0 OR 1.') 

otto 


TYRf 340 

otit 


ACCtRT 40*C08f»M<2»l> 

Ott2 

150 

90RHAT(' CNTIR CtU> 0 OR 1.') 

Otl3 


TYRE 350 

Olt4 


ACCERT 40»CLU»AA<3»2> 

otts 

340 

RORHATl' ENTER CLAi 0 OR 1.') 

OtM 


TYRE 340 

0tf7 


ACCERT 40»CLA»AA<3»3> 

OIM 

370 

FORNATl' ENTER CLDEi 0 OR 1.') 

Oltf 


TYRE 370 

OtfO 


ACCERT 40«CLDEi88i1il> 

0191 

380 

FORNATt' ENTER CHAD.') 

0192 


TYRE 380 

0193 


ACCERT 40 t CHAD 

0194 

390 

FORHATt' ENTER CNO* 0 OR 1.') 

019S 


TYRE 390 

0194 


ACCERT 40>CNQ*AA<lf 1) 

0197 

400 

FORHATl' ENTER CNUt 0 OR l.'> 

0190 


TYRE 400 

0199 


ACCERT 40rCHUfAA(l>2> 

0200 

410 

FORHATt' ENTER CNTU.') 

0201 


TYRE 410 

0202 


ACCERT 40 f CNTU 

0203 

420 

FORHAT(' ENTER CNA» 0 OR l.'> 

0204 


TYRE 420 

020S 


ACCERT 40>CNA*AA(lf3> 

0204 

430 

FORNAT< ' ENTER CHTA. ' ) 

0207 


TYRE 430 

0201 


ACCERT 40»CNTA 

0209 

440 

FORHAT(' ENTER CNDEt 0 OR 1.') 

0210 


TYRE 440 

0211 

C 4 9 9 • 

ACCERT 40tCHDE*8l(l>l> 

DEFINE D1MCN8I0NAL DERIUATIVES 

0212 


01 - RH09U19U1/2.0 

0213 


XU B Q1«8X(CXTU42.9CXT1-CDU-2.«CD1>/<AN88«U1> 

0214 


XA a.0l«t«(C0A-CLl)/AHS8 

021S 


XDE b.qusxcDDE/AHSS 

0214 


ZU B-Ql«|»(CtU42.*CLl>/(AH88»Ul«Ul) 

0217 


ZA b.Q1«8«(CI.A4CD1>/(AH88«U1> 

0218 


ZDE ■-Q1«8«CLDE/(AH88«U1> 

0219 


AND • Ql«S«C8AR«CiAR«(CNAD4CHa)/(2.«AIY*Ul> 

0220 


ANU • Q1«8«CIAR«(CHU4CHTU)/(AIY«U1) 

0221 


ANA - 01«8«CIAR«(CNA4CNTA)/AIY 

0222 

C 9 • • 9 

ANDE B 01989C8AR9CNDE/AIY 
DEFINE A NATRIX ELEHENTS 

0223 


A(ltl>-ANQ 

0224 


A<lt2)>AHU 

0223 


A(lt3)BANA 

0224 


A(1>4)b0.0 

0227 


A(2f n-o.o 

0228 


A(2>2)bXU 

0229 


A<2f3)BXA 

0230 


A(2»4)b-C08THA«32.174 

0231 


A(3»1)b1.0 

0232 


A<3f2)BZU 

0233 


A(3f3)BZA 

0234 


A( 3 » 4 ) >-8INTHA»C08RHI932 . 174/Ul 

0235 


A(4f 1 )bC08RHI 

0234 


A(4f2)B0.0 

0237 


A(4f3)B0.0 

0238 

C 9 • • • 

A(4t4>B0»0 

DEFINE B NATRIX ELEHENTS 

0239 


B(1>1)bAHDE 

0240 


B(l>2)-0.0 

0241 


B( 1 i3)b0.0 

0242 


B(2f 1 )rXDE 

0243 


B(2f 2)-0.0 

0244 


B(2t3>>0.0 

0243 


8<3>1)*ZDE 

0244 


B(3t2)B0.0 

0247 


B(3»3)«0.0 
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0241 

024f 

02S0 


I 

I 


02SS 

0234 

0287 

0231 

0239 

0240 

0241 


0242 

0243 

0244 
0243 
0246 
0267 
0248 

0269 

0270 

0271 

0272 

0273 

0274 
0273 

0276 

0277 

0278 

0279 

0280 
0281 
0282 

0283 

0284 
0283 
0286 

0287 

0288 

0289 

0290 

0291 

0292 

0293 

0294 
0293 

0296 

0297 

0298 

0299 

0300 


• (4iU«0.0 
•<4t2)«0.0 
•<4»3)"0«0 

C.... ALL ELCHCNT8 09 THE AA MATRIX ARC DCF2MCD 
C.... DEFINE ADDITIONAL CLCMCNT8 OF THE il MATRIX 
ID<lt3>>1.0 
•D(2»J)al.O 
Dt<3i3)>1.0 
8i<4»3)«lfC 

C.... DEFINE AP MATRIX <AltUK<CD ORDER OF THE 0I8ERVATI0N VECTOR III 
C.... Op ALPHAf THCTA> ODOV t AX-AXDIAI* AND AN-ANIIAI) 

DO 430 I-liS 
DO 440 J«l>4 
AP<IiJ)-1.0 
440 CONTINUE 
480 CONTINUE 

AP<4»2)>1. 0/32. 174 
AP<7r3>—Ul/32.174 

C.... DEFINE DP MATRIX (ASSUMED ORDER OF THE CONTROL VECTOR 18( 
C.... DCi 8IA8) 

DO 470 I>lf8 
DO 4S0 J-lf3 
RP(I»J>-1.0 
480 CONTINUE 
470 CONTINUE 

DO 490 I>lf3 
8P<4»I)-1. 0/32. 174 
BP(7>I)«-U1/32.174 
490 CONTINUE 

C.... SKIP LATERAL DIRECTIONAL INPUT CASE 
80 TO 700 

C.... LATERAL DIRECTIONAL CASE 
300 CONTINUE 

310 FORMAT!' ENTER CLPr 0 OR 1.') 

TYPE 810 

311 FORMAT!' ! 1 IF THIS VARIESi 0 OTHERWISE )') 

TYPE 311 

ACCEPT 60>CLPrAA!lfl) 

320 FORMAT!' ENTER CLRt 0 OR 1.') 

TYPF 320 

ACCEPT 60fCLR>AA!l>2) 

330 FORMAT!' ENTER CLDt 0 OR 1.') 

TYPE 330 

ACCEPT 60fCLD>AA(l>3) 

340 FORMAT!' ENTER CLDA> 0 OR 1.') 

TYPE 340 

ACCEPT 60>CLDA>BB!lrl) 

380 FORMAT!' ENTER CLDRr 0 OR 1.') 

TYPE 530 

ACCEPT 60tCLDRfBB!lf2> 

560 FORMAT!' ENTER CNP> 0 OR 1.') 

TYPE 360 

ACCEPT 60pCNPf AA(2il> 

370 FORMAT!' ENTER CNRt 0 OR 1.') 

TYPE 370 

ACCEPT 60iCNR.AA;;2.‘2) 

380 FORMAT!' ENTER CNBt 0 OR 1.') 

TYPE 380 

ACCEPT 60pCNBpAA!2i3) 

390 FORMAT!' ENTER CNOAp 0 OR 1.') 

TYPE 390 


0301 ACCEPT 60pCNDA>BB!2»1) 

0302 600 FORMAT!' ENTER CNDR» 0 OR 1.') 

0303 TYPE 600 

0304 ACCEPT 60»CNDRpBB!2f 2) 

0305 610 FORMAT!' ENTER CYBi 0 OR 1.') 

0306 TYPE 610 

0307 ACCEPT 60fCYBtAA!3»3) 

0308 620 FORMAT!' ENTER CYOAt 0 OR 1.') 

0309 TYPE 620 

0310 ACCEPT 60pCYDAfBB(3f 1) 

0311 630 FORMAT!' ENTER CYDRt 0 OR 1.') 

0312 TYPE 630 

0313 ACCEPT 60»CYDRf BB!3»2) 
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0314 

0315 
03U 
0317 
0310 
031f 

0320 

0321 

0322 

0323 

0324 

0325 
0324 
0327 

0323 

0329 

0330 

0331 

0332 

0333 

0334 

0335 
0334 
0337 
0333 

0339 

0340 

0341 

0342 

0343 

0344 
0349 
0344 
0347 
0343 

0349 

0350 

0351 

0352 

0353 

0354 
0359 


0354 

0357 

0358 

0359 


0340 

0361 

0342 

0343 
0364 
0345 


0344 

0347 

0348 
0347 

0370 

0371 

0372 

0373 

0374 

0375 
0374 
0377 


C.... OfMNC DINCNflONAL DlftlVATIOIf 
01 > RH0«Ul«Ul/2.0 

ILI* • 0im8MNfffAN«CL9/(2.*AlX«Ul) 
iLR • Ql«8ffPAN8IPAN«CLR/<2.«AIX»Ul) 

111 • ai8f88PAN«CLI/4!X 

INP > QlS«8tPAN«tP4NICNP/(2.94IZtUl) 

•NR ■ QimtPAN»fPAN«CNR/(2.«AlZ9Ul) 

•Nl • ai«888PAN*CNt/AlZ 
'«• ■ Q19MCYi/(AN'*5«Ul) 

vOA ■ Ql«t8fPANACl.ii,4/AXX 
ilDR • Ql«ft8PAN«Cl.DR/AXX 
•NDA ■ Ql«t«8PANtCNDA/AXZ 
•NDR • Ql«f»tPAN«CNOR/AIZ 
YDA ■ Qlff*tPAN«CYDA/(AHII«Ul> 

YDR > Q1M*SPAN«CY0R/(AH8I«U1) 

C.... DCPXNC A NATftXX ELCHCNTI 
A<lf l)>iLP 
A<lr2>-il.R 
A(1»3)«SI.I 
A<li4)>0.0 
A<2f 1>->NP 
A(2»2)-0NR 
A<2f3>*tNI 
A<2i4>«0.0 
A(3>l)>8XNAt.P 
A<3f2)«-C08ALP 
A(3>3)>Yi 

A(3i4>-32.174«C08rHA«C0SPHI/Ul 

A<4tl)>l,0 

A<4t2)>C08PHX»TANTHA 

A<4i3)>0.0 

A<4>4)-0.0 

C.... DEFINE 8 HATRXX ELEMENTS 
•(ltl)«iLDA 
•(lt2)>»LDR 
Klt3)>0.0 
l<2>l)>8NDA 
•(2i2)alNDR 
•(2t3)>0.0 
D(3f l>-rOA 
•(3>2>«V0R 
•(3»3)>0.0 
•(4«1)>0.0 
B(4f2)>0.0 
I(4f3)>0.0 

C.. ■ ALL ELEMENTS OF THE AA MATRIX ARE DEFINED 

C...« DEFINE ADDITIONAL ELEMENTS OF THE 8B MATRIX 
•B( lf3)-l .0 
•D(2i3)-1.0 
BB(3t3)>1.0 
BB(4>3)>1.0 

C.... DEFINE AP MATRIX (ASSUMED ORDER OF THE OBSERVATION VECTOR 181 

C.... P> R> BETAf PHIr PDOTr RDOT> AND AY-AYBIA8) 

DO 450 I-lf4 
DO 440 J«l»4 
APdfJl-l.O 
440 CONTINUE 
450 CONTINUE 

AP(7.3)-UI/32.174 

C.. . DEFINE AP MATRIX (ASSUMED ORDER OF THE CONTROL VECTOR IS! 

C>... DA> DRt AND BIAS) 

DO 470 I-l»4 
DO 480 J-l>3 
BPdf Jl-l.O 
680 CONTINUE 
470 CONTINUE 

DO 690 I«lf3 
BP(7d)«Ul/32.174 
490 CONTINUE 
700 CONTINUE 
MU -3 

NZ -4 

MY ■? 

C.... ECHO DATA BACK 
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0371 


037t 


0310 

0312 

0304 


0315 

038A 


0387 


0388 

0389 

0390 

0391 
07'.»2 

0393 

0394 
0393 

0396 

0397 
0393 

0399 

0400 

0401 

0402 

0403 

0404 

0405 

0406 

0407 

0408 

0409 


0410 

0411 


710 rORNAr<* AimANC ZNAUT 0ATA'«/i 

!>/* 

!•/» 
l»/t 
l»/» 
l>/t 
l»/» 
l»/f 
l»/» 

!>/• 
l»/» 

)i/i 
)»/f 
)t/i 
li/> 

!»//> 

’»AIX» 

ItCOlfCXTl 

C.... SPLIT FOR CASES 

IF(CASE.CQ. 'LONG' >00 TO 790 
tFCCASE.CO.'LATR'lOO TO 890 

STOP 

C.... LONOITUOINAL CASE 

C..«. URITE OUTPUT TO DISPLAY* OET UEZOHTINO FACTORS* AND FINISH FILE 
790 CONTINUE 


0 

WINi AREA 

(IN FT8S2) 

• 'fF12.4t 


WCIOHT 

(IN LtS) 

■ '*F12.4* 

4 

WINt SPAN 

(IN FT) 

■ '*F12.4» 

4 

CIAR 

(IN FT) 

■ '*ri2.4* 

4 

AIRSPEED 

(IN FT/SEC) 

■ '»F12*4. 


DENSITY 

(IN SLU0/FT«»3) 

■ '*F12.4» 


ALPHA 1 

(IN RAD) 

■ '*F12.4* 


THETAl 

(IN RAD) 

■ '*F12.4* 


PHIl 

(IN RAD) 

• '.F12.4* 


lYY 

(IN SLU8tFTt«3) 

> '*F12.4* 

4 

I XX 

(IN SLU0«FT«t2> 

• '.F12.4* 

4 

ZZZ 

(IN SLUOtFTSSZ) 

• '*F12.4* 

$ 

CLl 


• '*F12.4* 

4 

CDl 


• '*F12.4* 

4 

CXTl 


■ '*F12.4* 

TYPE 710.S.WCIOHT 

*SPAN*CtAR*Ul*RHO*ALP*THA*PKI*AIY* 

> 

AlZ.CLl* 

CDl.CXTl 



760 FORHAT(' 


LONOITUOINAL DERIVATIVES'*/* 

4 

9 

CDU 

' *FS.4* ' 

XU 

' *FS.4*/* 

. ' 

CXTU 

' >F8.4*/* 




CDA 

' *F8.4* ' 

XA 

'*FS.4*/* 

. ' 

CODE 

'*F8.4»' 

XDE 

' *F8*4*/* 

. ' 

CLU 

'*F8.4*' 

ZU 

'*FS.4*/* 


CLA 

' *FS.4*' 

ZA 

'*FS.4*/* 

. ' 

CLDE 

' *F8.4* ' 

ZDE 

'*F8.4*/* 

. ' 

CHU 

' *F8.4* ' 

NU 

'*F8.4*/* 


CHAD 

' *F8.4*/* 



. ' 

CHQ 

'*F8.4*' 

HQ 

'*F8.4*/* 

. ' 

CHA 

'.F8.4*' 

HA 

'*F8.4*/* 

» 

• 

CHDE 

' *F8.4* ' 

HDE 

'*F8.4*/'//) 

TYPE 760 

*CDU*XU*CXTU*CDA* 

XA.CDDE*XDE*CLU*ZU 

*CLA*ZA* 

6 


CLDE*ZDE*CNU*AHU 

*CHAD*CHQ*AHQ*CHA* 

AHA* CHDE *AH0E 


C • • • • 
770 


780 

790 

800 

810 

820 

830 


C • • • t 

c • • • • 

c « • • • 

890 

860 


OET THE UEIOHTINO HATRIX DIAGONAL VALUES 
FORHAT<' ENTER THl 
TYPE 770 

ACCEPT 60*D1(1*1) 

FORHAT(' ENTER THE WEIGHTING FACTOR FOR VELOCITY.') 

TYPE 780 

ACCEPT 60*01(2*2) 

FORMAT(' ENTER THi 
TYPE 790 

ACCEPT 60*01(3*3) 

FORNAT(' ENTER THi 
TYPE 800 

ACCEPT 60*01(4*4) 

FORHAT(' ENTER THi 
TYPE 810 

ACCEPT 60»D1(S*S) 

FORHAT(' ENTER THi 
TYPE 820 

ACCEPT 60*01(6*6) 

FORNAT(' ENTER THI 
TYPE 830 

ACCEPT 60*01(7*7) 

SKIP PAST LATERAL DIRECTIONAL CASE 
GO TO 990 

LATERAL DIRECTIONAL CASE 

WRITE OUTPUT TO DISPLAY* OET WEIGHTING FACTORS. AND FINISH FILE 
CONTINUE 

FORHAT(' LATERAL DIRECTIONAL DERIVATIVES'*/* 


WEIGHTING 

FACTOR 

FOR 

0' > 

weighting 

FACTOR 

FOR 

VELOCITY 

WEIGHTING 

FACTOR 

FOR 

ALPHA.'} 

WEIGHTING 

FACTOR 

FOR 

THETA. ' > 

WEIGHTING 

FACTOR 

FOR 

QDOT. ' > 

WEIGHTING 

FACTOR 

FOR 

AX. ' > 

WEIGHTING 

FACTOR 

FOP 

AN. ' ) 


' CYB 

'*F8.4*' 

YB 

' *F8.4*/. 

' CYDA 

' *F8.4* ' 

YDA 

' *F8.4*/* 

' CYDR 

'*F8.4*' 

YOR 

' *F8.4*/* 

' CLB 

' *F8.4* ' 

LB 

' .F8.4,/, 

' CLP 

' *F8.4* ' 

LP 

' *F8.4*/* 

' CLR 

' *F8.4* ' 

LR 

' *F8.4*/* 

' CLOA 

' *F8.4* ' 

LDA 

' *F8.4*/. 
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0412 


0413 

0414 

0419 
0414 
0417 
0411 
041» 

0420 

0421 

0422 

0423 

0424 
0429 
0424 
0427 
0429 

0429 

0430 

0431 

0432 

0433 

0434 

0439 

0434 

0437 

0438 

0439 

0440 

0441 

0442 

0443 

0444 
0449 

0446 

0447 
044B 
0449 

0490 

0491 

0492 

0493 

0494 
0499 

0496 

0497 

0498 

0499 

0460 

0461 

0462 

0463 

0464 
0469 

0466 

0467 

0468 


* 

CLDR 

'iF8.4»' 

LDR 

'*F8.4*/* 

t 

CNB 

'»F8.4.' 

NB 

' •F8.4*/* 

t 

CNR 

'.F8.4*' 

NR 

'*F8.4*/* 

§ 

CNR 

•tF8.4f ' 

NR 

' *F8.4»/* 

* 

CNDA 

'»F8.4.' 

NDA 

' iF8.4*/f 

» 

CNDR 

' »F8.4t' 

NDR 

' *F8.4t////) 


TY9E t«OiCyityi»CYDA>rDA»CYORiYD8>CLitBt8»CLPi8l.PfCLRtlLR» 
CLOA » 8L0A » CLDR • 8I.0R i CNI > tNR t CNR > t*NR I CNR > INR • CNDA f 
•NOA»CNOR»BNDR 


C.... OCT THE UCIONTINO MATRIX DIAGONAL 0ALUC8 
870 FORNATl' ENTER TNC UCIONTINO FACTOR FOR F.') 
TYRE §70 


880 

890 

900 

910 

920 

930 


ACCCRT 60»DK1«1> 
FORMAT!' ENTER THE 
TYRE 880 

ACCCRT 60tDl(2f2) 
FORMAT!' ENTER THE 
TYRE 890 

ACCERT 60fDl!3>3> 
FORMAT!' ENTER THE 
TYPE 900 

ACCERT 60f0t!4f4) 
FORMAT!' ENTER THE 
TYPE 910 

ACCCRT 60fDl!9»9> 
FORMAT!' ENTER THE 
TYRE 920 

ACCERT 60tDl!6i6) 
FORMAT!' ENTER THE 
TYPE 930 


UCIONTINO FACTOR FOR R.') 
UEIOHTINO FACTOR FOR BETA.') 
UCIONTINO FACTOR FOR PHI.') 
UEIOHTINO FACTOR FOR PDOT.') 
UEIOHTINO FACTOR FOR ROOT.') 
UEIOHTINO FACTOR FOR AY.') 


ACCEPT 60>D1!7>7) 

930 CONTINUE 

C.... URITE MATRICES TO FILE 
960 F0RMAT!2I10) 

1160 F0RMAT!7F12.6> 
URlTE!lr960)4>4 


DO 970 I-l>4 

URITE!lrll60)!A!It J) >J-1>4) 
970 CONTINUE 

URITE!l>960)4fMU 
DO 980 I-lf4 

URITE!lrll60)!B(It J)» J-lfMU) 
980 CONTINUE 

URITE!lr960)4f4 
DO 1000 I>li4 

URITE!l>1160>(AA!lt0)f J-l>4) 
1000 CONTINUE 

URITE! 1 t960)4>MU 
DO 1010 I-lf4 

URI TE ! 1 • 1 1 60 > ! BB ! I f J ) f J- 1 1 MU > 
1010 CONTINUE 

URITE!lf960)7>4 
DO 1020 I>1«7 

UR1TE!1>1160)(AP(I>J)*J-1*4) 


1020 CONTINUE 

URITE!lf960)7tMU 
DO 1030 1-1*7 


URITE! 1 *1160)! DR! I* J> >J-1 *MU) 
1030 CONTINUE 

URITE(l»960>7f7 
DO 1040 I-li7 

URITE! 1 > 1160)!D1! I* J> f J-1 r7> 
1040 CONTINUE 

URITE!l*1160>!ZER0!I>*I-lr4> 
URITE! 1* 1160) !DIAS(nf 1-1*4) 
STOP 


END 
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A. 7. 2) MAIN MMLE PROGRAMS 

Description: The mein program of the MMLE BONES routines 

acts as a controller In calling the subroutines as needed. 
Initially, It reads the Input data for the starting conditions 
of the test case. If all states are measured a least squares 
process Is used to compute the Initial estimate of the derivatives. 
If the states are not completely measured this feature must 
be skipped over or errors In the solution of updates to the 
coefficients will result. 

Flowchart: 


START 



STOP 
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Program listing 


0001 


0002 


0003 


0004 

0005 

0006 

0007 

0008 

0009 

0010 


0011 

0012 

0013 

0014 


0015 

0016 

0017 

0018 


0019 


C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

0 

c 

c 

c 

c 


c 


c 

c 

c 

3511 

3512 


ROQRM MAIN 


:S: 


■0NC8 - FRL 


MAIN PROORAN OF THE HAXINUH LIKELIHOOD CtTIHATOR 
TECHNXQUEf <NHLE). THIS PROORAN IS DERIVED FRON 
THE *SONES* PROORAH THAT HAS ORIOINALLY DEVELOPED 
lY NASA. THE FOLLOUINO SUSROUTXNES ARE REQUIRED 
FOR THE OPERATION OF THIS PROORAH 1 
OIRLf EAT> CRAHERt SPITl* REDUCE# HULT# OUTPUT 
ADO# HAKE# ZOT> LOAD# LOADl# SPIT# SOLVE# AND 
DIAOIN. 

THE OUTPUT OF THE PROORAH IS TUO FILES THAT 
CONTAINS THE HATRICES CAI AND C13 FOR EACH 
ITERATION AND THE ESTINATED TIHE RESPONSES 
RESPECTIVELY. SEE THE SPECIFIC INSTRUCTIONS 
OF THIS PROORAN FOR FURTHER INFORHATION 
CONCERNING THE INPUT AND OPERATION OF THIS HHLE 
PROORAH. 

THIS HODIFIED ‘SONES* PROORAH HAS HRITTEN DYt 
ALEX KOTSABASIS 

DATE 22<-N0V>80 


NEUTON-RAPHSON HETHOD FOR OBTAINING STABILITY DERIVATIVES 
LONOl ALPHA# Q# V# THETA# AN# QDOT# AX 
L-D 1 P> R# BETA# PHI# PDOT# RDOT# AY 

COHHON HAX#MA»nAH#HAT.Z#U#02#El#APHX>DUH#PHIl#Dl#A*B#AA#BB# 

2 BJI#XJI#SUH#FB#XT1# ZERO* DS4#0D4#E#XTX#CCC# BIAS# 

3 IZE#XBIA8#IC#XLA#APR#NAPR#XT4#JKKH#XTS#AP#BP 
COHHON HH#ENN#HX#HUHX#RB#IA# JK#NB#EPS#P0B#I#HUMX1# 

2 TRACE#K#IJ#TINE#RDB#MXP1#JKH>LH#TT#JKHH1> 

3 ALPHA# NNH1#HZ#J#KJ#BD#LL#L#NN#HU#HUX#EN#PP#AAA#KH#FAC> 

4 ITR#N1#XT2#XT3#HZH1 
COHHON/HATAB/ALX # BLX #ERX 
COHHON/CON8T/ KABC#KCDF 
COHHON/TRNSFR/ XL(300#7) 

DIMENSION ALX(20#10> #BLX(20#10)#ERX<10#10> 

DIMENSION XT5(25) #APR<25) 

DIMENSION AP(8#4># BP{8#3># XT4(4) 

DIMENSION Z(7#3) #U(3#3) #D2(7)#DD4(S#4) #BIAS<5> #APHI(5#4> # 

2 XT1(7>#PHI1(5#4>#D1<8#7)#A<S#4)#B<5#4)#AA<5#4># 

3 BB(S#4)#BJI(2S#4> #XJI<2S#7)#8UH(25#2S> #PB(2S) # 

4 DUM(25#4> #XT2(7) #ZER0(S)#DS4(S#4) #XT3(7> 

BYTE INAHE(1S)#CONNT(00#4) 

C0MM0N/AN8UER/ KBUOO 

BYTE L00#DIR#AN8 
DATA LOO#DIR/'L'#'D'/ 


TYPE 3511 

F0RMAT<///////////#20X#'« * * THE MMLE PROGRAM * * *'//) 
TYPE 3512 

FORHATdOX# 'THINGS YOU HAVE TO KNOW TO RUN THIS PROORAMl'#/ 


1 #10X#'1.IS IT A LOGITUDINAL OR A LATERAL DIRACTIONAL RUN?'# 

2 /#10X# '2.NAME OF FILE UITH INITIAL CONDITIONS.'/ 

3 #10X#'3.NAME OF FILE UITH MEASURED DATA.'/ 

4 #10X#'4.H0U TO DESIGNATE OUTPUT FILES.'////) 

C 

TYPE 3513 
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0020 3813 FORMAT</////flOX>' INDICATE TYNE OF RUNt'*/ 

1 ilOXt'IF LONGITUDINAL TYPE 

2 tlOXf'IF LATERAL-DIRECTIONAL TYRE *0* 

3 flOX» 'DELECT RUNi') 

C 


0021 


RCA0(8f3814) AN8 

0022 

3814 

C 

C 

FORMAT (Al) 

0023 

8 

CONTINUE 

0024 

S3 

CONTINUE 

0028 


NI • 28 

0026 


XXXX - 1.0 

0027 


NAX > 8 

0028 


NA B 4 

0029 


HZ ■ 7 


C 

C 


C ATTACH DATA FILE CONTAINING HATRICES AND INT. CONSTANTS. 


0030 


TYPE 119 

0031 

119 

FORMAT!//* 10X» 'ENTER DATA FILE NAME WITH INITIAL CONDITIONS*'* 
1 /*10X*'AND MATRICES A* D* AA* AP* ETC.') 

0032 

128 

FORMAT (14A1) 

0033 


ACCEPT 128* ( INAME ( 1 ADC) * lADC-l * 14 ) 

0034 


INAME(15)>0 

0038 


0PEN<UNITb2*NAME«INAHE*TYPE>'0LD' *ACCES3b'8EQUENTIAL'* 
1 READ0NLY*F0RM>'F0RMATTED'*REC0RD3IZE>132) 


C CLEAR OUT OLD FILE NAME 

0036 


DO 127* IADC«1*14 

0037 


INAME! IABOb' ' 

0038 

127 

CONTINUE 

0039 

1700 

FORMAT!10X*10E12.4) 

0040 

3777 

F0RNAT!12X*7I10) 

0041 

3700 

F0RMAT!12X*8F10.4) 

0043 

700 

F0RMAT!8F10.4) 

0044 

777 

F0RHAT!7I10) 

0043 

1010 

FORMAT! 10E12. 4) 

0046 

1011 

F0RMAT!I3*9E12.4) 

0047 

1012 

F0RMAT!12E12.4) 

0048 


FACT -1.0 

0049 


BLANC - 0.0 


C READ C0HNENT8 FROM INPUT FILE <4 LINES OF 80 CHARACTERS) 

0050 


READ! 2* 1301) ! !COHNT ! 1 * J) * I>1 *80 ) * J>1 * 4 ) 

0081 

1301 

F0RHAT!80A1) 


C READ 8TATMENTS 

0032 


READ!2*777) NN* ITR*HZ*MAPR 

0053 


READ ! 2 * 700 > HH * EPS * T I ME * ALPHA * XLA 


C LOAD MATRICES 

0054 


ENN - NN 

0035 


KCDF > NN 

0056 


KABC > ITR 

0057 


CALL L0A0!4*A*B*AA*BB) 

0058 


MAX • 8 

0059 


MA ■ 7 

0060 


CALL L0A0!3*AP*BP*D1*01) 


C 

C 

c 


c 


0061 


MAX - 5 

0062 


MA » 4 

0063 


NNHl - NN-1 

0064 


MU ■ B!HAX*2)4’ .01 

0065 


MX > A!MAX*2)6 .01 


C 

READ IN ZEROS' AND 3IA8E8. 

0066 

C 

READ !2*700) ! ZERO! I ) * I-l * MX ) * ! BI AS! lA) * IA-1 *MX) 


C 

c 

C 

0067 CL08E<UNIT-2) 
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OOAt 

004f 

0070 

0071 

0072 

0073 

0074 

0075 
0074 
0077 

007t 

007» 

OOtO 

OOfl 


0012 

0013 

0014 
0085 
0084 

0087 

0088 
0087 
0090 


0091 


0092 


0093 

0094 

0095 
0094 

0097 

0098 

0099 

0100 
0101 

0^02 

0103 

0104 

0105 
0104 
0107 


0108 

0109 


0110 

0111 


0112 




c 

NXPl • HX 4 1 
nZHl ■ NZ > 1 
NUX • HU 4 NX 
NUHX > NU*HX 
NUNXl HUNX 4 1 
YY • 0.0 
XX • 1.0 

C ADO 8XA8C8 AND ZCK08' 

DO 4f lalfHX 
XT4<X> ■ 0.0 
XT3(X) > 0.0 

C XX ■ XX 4 ZCRO(X) 4 iXA8(X) 

XX » XX 
DO 48 J-1»HU 

YY > YY 4 AA<X>J) 4 88<XfJ> 

48 XX > XX 4 AA<XfJ> 4 88(X>J) 

C 

C ........ ...... 

C 

YY ■ YY 4 AA(X.NX) 

49 XX ■ XX 4 AA(XrNX) 

C 

JKNN ■ YY 4 .01 
JKN ■ XX 4 .01 
JKHNl ■ JKH - 1 
8UH<NXfl> > JKH 
8UH<NX>2) > JKH 
HAX • NX 
HA - NX 

C ...... ..... 

c 

C INXTIALIZE HATRXCE8 TO ZERO 
C 

CALL Z0T(8UH) 

C SELECT APRXORX OPTXON THRU HAPR. 

C 

IF<HAPR) 174>178»177 
C 

C READ IN APRXORX HATRXX 
C 

177 DO 261 XB-lfJKH 
DO 643 lA-l.JKH 
8UH<XBf XA)-0.0 

663 CONTINUE 

DO 261 XA-lfJKH 
261 8UH(XBf IA)-SUH(XB>XA) 

APRUB) - SUH<XB>XB) 

80 TO 178 
176 CONTINUE 
C 

DO 664 XA«1»JKHH1 
APR(IA)>0.0 

664 CONTINUE 

DO 263 XAaltJKH 
263 APR<IA) > APRdAXPACT 

178 CONTINUE 
C 

C ENTER NAHE OP DATA PILE UXTH HEA6URED PLIGHT TEST DATA 
C 

TYPE 139 

139 PORHAT(/7flOXt 'ENTER PILE NAHE CONTAINING THE HEASUREO DATA') 
C 

C ATTACH 8TATHLNT POR PILE CONTAINING HEASURED DATA 
C 

ACCEPT 128. ( INAHE< lABC) . IABC-1 . 14 ) 

0PEN(UN1T>4.NAHE>INAHE.TYPE-'0LD' .ACCES8>'SEQUENT1AL' . 

1 PORH-'UNPORHATTED'. READONLY) 

C REUIND TAPE 
REWIND 4 
C 
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c 

c 

C HEAD IN DATA AND MINT OUT INITIAt. CONDITIONS 
C SIASCS AND ZERO'S 
C 

0U3 1303 FORHAT<'- - -- -- -- -- -- -- -- 

1 ' - 

OlH 1303 F0RNAT<24X>' INITIAL CONDITIONS 

Oils 1304 F0RNAT(10X»'NUHiER OF DATA MINTS I '»I3i 

1 SXt'NAXINUN NUNIER OF ITERATIONS i 'fI3> 

2 /rlOXt'OATA SANFLINO INTERVAL t '>F10.4r 

3 5X» 'FIRST DATA POINT AT TINE t '«F10.4» 

4 /flOXt 'DIAGONAL HULTIPLTINO FACTOR t 'iFlOtA* 

5 SX» 'NUNIER OF STATES I '»I3*/> 


0114 

1328 

FORMAT dOXf' DIAGONAL ELEMENTS OF THE WEIGHTING MATRIX Dll'f/ 
1 8X>7F13.3*/) 

0117 

1324 

FORMAT</UOXi 'ESTIMATES OF THE CA3 AND CI3 MATRICES' p/) 

0118 

Q 

1308 

FORMATUOXp'INITIAL input matrices CA3 AND CI3.'p/p 

1 IOXp'A star («> FOLLOWING THE VALUE OF A MATRIX' p/p 

2 IOXp 'ELEMENT INDICATES THAT THE RESPECTIVE DERIVATIVE 

3 IOXp 'IS NOT ESTIMATED IT THE NMLE METHOD.'/) 

OUT 

1304 

FORriATC/plOXp'STAilLITY MATRIX CA3') 

0120 

1307 

FORNAT(/p IOX p 'CONTROL MATRIX CI3') 

0121 

1309 

F0RMAT<10Xp 'ITERATION' pI3p ' WAS COMPLETED'/) 


C 

PRINT OUT INPUT DATA 

0122 

w 

PRINT 1302 

0123 


DO 1407 J>1p4 

0124 


PRINT 1421p (COMNT(IpJ)pI-IpBO) 

0125 

1421 

F0RMAT(10Xp80A1) 

0124 

1407 

CONTINUE 

0127 


PRINT 1302 

0128 


PRINT 1303 

0129 


PRINT 1304p NNpITRpHHpTIMEpXLApMZ 

0130 


PRINT 1302 

0131 


PRINT 703 

0132 

703 

FORMAT </p IOX p 'ZEROS AND BIASES') 

0133 


PRINT 1700p (ZER0<I)pI>1pMX) 

0134 


PRINT 1700p (DIAS(IA)pIA-IpMX) 

013S 


PRINT 1328p(D1(IBCDpIBCD)pIBCD>1p7) 

0136 


PRINT 1302 


C SET MAX AND MA TO CA3 AND CB3 DIMENSIONS 

0137 


MAX>8 

)138 


HA -4 

0139 


PRINT 1308 

0140 


PRINT 1306 

0141 


CALL SPITKApAApI) 

0142 


PRINT 1307 

0143 


CALL SPITKBpBBpI) 


C 

c 

c 

c 

c 

C STARTING ITERATION LOOP 


0144 

C 

TT • TIME - HH 

0148 


DO 1 LM>1pNN 

0146 


TT ■ TT f HH 

0147 

c 

U(MUp1) ' 1.0 

0148 

« 

CONTINUE 


c 


0149 

c 

DO 272 lA >1 pNI 

0180 


XT8(IA) > 0.0 

0181 

272 

PDUA) > 0.0 

0182 


IZE «1 

0183 


DO 276 IA»1pMX 

0134 


IF< ZERO<I A) )277p 276.277 

0188 

277 

IZE - IZE 4 1 

0186 

276 

CONTINUE 
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01S7 


OlSt 

OlSf 

OlAO 

01*1 

0142 

0U3 

0144 

0145 
0144 
0147 
014G 


0149 

0171 

0172 

0173 

0174 

0175 
0174 


0177 

0178 

0179 

0180 
0181 


01B2 

0183 

0184 

0185 
0184 

0187 

0188 

0189 

0190 


0191 

0192 

0193 

0194 

0195 


0194 

0197 

0193 

0199 

0200 
0201 


C 

C - 

c 

C HAXN LOOP FOR 'INT* NUHIIR OF ITERATIONS 
C 

DO 12 Lt ■ IfITR 
C 

C RCUIND TAFC FOR EVERY ITERATION STEF 
C 

REMIND 4 

C 

TYFE 345» LU 

345 FORMAT (//' NUNSER OF ITERATION PRESENTLY CONFUTED l'rI3i//> 
C 

MAX ■ 5 
NA • 4 

DO 31 JK aliS 
DO 31 IK>1>4 
KT ■ XK44«(JK-1> 

ALX(KTtLL) - A(JK*IK> 

DLX<KT>LL> > i(JKtlK) 

31 CONTINUE 

C 

C 

C CALL SPECIAL MATRIX OUTPUT ROUTINE 

C 

ir<LL.EQ.l} 00 TO 1308 
PRINT 1324 
PRINT 1304 
CALL SP1T1(A>AA>LL) 

PRINT 1307 

CALL SPITKBfDDtLL) 

1308 CONTINUE 

C 

C 

MAX ■ 5 
MAT > 5 
MAM ■ 4 

CALL EAT (AfHHfPHIlrAPHIfD54>DD4) 

U<3rl) • 1.0 
C 
C 

U(3>2) - 1.0 
U<3r3) - 1.0 
XJKNItl) ■ JKN 
XJI(NI»2) - MX 
DJKNIil) • JKN 
BJI(NI>2) > MX 
8UM(NIfl) > JKN 
SUM(NI»2) - JKN 
NA • NI 
C 

C INITIALIZE AND READ DATA FROM TAPE 
C 

DO 778 IJK-lrJKM 
DO 778 JKL-lfIJK 
778 SUM<IJK>JKL> « 0.0 
MAX > NI 
CALL ZOT(XJI) 

C 

C READ IN THE FIRST TUO SERIES OF MEASURED DATA 
C FROM THE DATA TAPE. 

C 

C 


XTK3) 

■ 

0.0 

XT2<3> 

a 

0.0 

XT1(5> 

a 

0.0 

XT2{5) 

a 

0.0 

XTl (4> 

a 

0.0 

XT2(4> 

a 

0.0 
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D 




0292 

0203 


0203 

0204 


0207 

0200 

020f 

0210 

0211 

0212 

0213 

0214 

0213 

0214 

0217 

0218 


C 

KlUOO - 1 

c 

IFlANt.tQ.LOO) 00 TO 2012 
C 

c 

READ <4) DXY>DXYtDXY»DXY>DXY»Xri<4)tXTl<l>>XTl(7>fXTl(2)» 
1 U(l»l>iU(2»l) 

READ (4) DXY»DXYtDXY»DXYrDXY»XT2(4>>XT2<l)>XT2(7)iXT2(2)» 
1 U<lf2)»U<2»2$ 

C 

c 

00 TO 2013 
2012 CONTINUE 

XT1(2) - O.C 
XT2(2) ■ 0.0 
XT1(3> > 0.0 
XT2<3> ■ 0.0 
XTU3) • 0.0 
XT2(3> ■ 0.0 
U(2.1> ■ 0.0 
U(2.2) - 0.0 
C 

KlUOO -0.0 

c 

READ (4) XTl(4)tXTl<l)rXTl<7>.XTl(4>>U<l.l)fDXY.DXYfDXY> 

1 DXYfOXY.DXY 

C 


0219 


READ (4) XT2(4)tXT2(l)tXT2<7>.XT2<4>fU<1.2).DXY>DXYtDXY. 
1 DXYfDXYtDXY 


C 

0220 2013 
C 

C 

0221 
0222 

0223 31 

0224 

0223 44 

0226 

0227 67 

0228 
0229 
C JO 
0^31 

0232 46 

0233 


CONTINUE 


IC ■ 0.0 
DO 31 I-l.HX 
XJI(JKHil> - XT2<1> 

IF(LL-l) 44.63.64 
DO 64 lA-l.NX 
IF<ZER0<1A>>67.64.67 
1C - 1C f 1 

XT3<IA> - XT3<1A) ^ FKJKH-IZE -f IC) 
XTKIA) > XTl(IA) 9 XT3(1A> 
XJKJKH.IA) - XJKJKN.IA) * XT3(IA) 
XT2(1A) > XJKJKH.IA) 

CONTINUE 
IC - 0.0 


C 

C ADD BIASES 
C 


0234 


DO 166 lA-l.HX 

0233 


IF(IIAS(IA))167.166.167 

0236 

167 

IC - IC f 1 

0237 


XT4(IA) > XT4(IA) 9 FB<JKMH4>IC) 

0238 


XTKIA) - XTKIA) - XT4UA) 

0239 


XT21IA) - XT2CIA) - XT4(IA) 

0240 


XJKJKH.IA) • XT2(IA) 

0241 

166 

CONTINUE 

0242 

65 

P 

CONTINUE 


W 

c 

HAIN HHLE LOOP 

0243 

U 

DO 260 lA-l.JKHH 

0244 


XTS<IA) • XT5(IA)4PB(IA) 

0243 

260 

CONTINUE 

0246 


DO 13 lA-l.HZ 

0247 


D2(1A) - 0.0 

0248 


Z(IA.l) - XTKIA) 

0249 


Z<IA.2) > XT2(IA) 

0250 

13 

CONTINUE 

0251 


IC >0.0 


C 

ZERO SPLIT 

0252 


DO 62 I>1.HX 


i " » 

I 
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02S3 


IF<ZCft0<I)>43«42f43 

0294 

43 

1C - IC91 

0299 


XJKJKN-tZC 9 ICiX) > l.C 

0294 

42 

CONTINUE 

0297 


CALC 0I9C 

0299 


NAX ■ NI 

0299 


NA - Nl 


C 1 

c 

9UT9UT or ITERATION LOOP 

0240 

c 

DO 329 IA-1*JKH 

0241 

329 

9UN<IA»IA> ■ 9UN<1A*IA)«XLA 


C 

CALL iriT<9UN) 

0242 


9UN(NI»1) ■ JKN-1 

0243 


8UH<NIf2> > JKN-1 

0244 


PRINT 1309tLL 

0249 


PRINT 1302 

0244 


Xr<LL-ITR) 249>248*248 

0247 

244 

CALL CRAHERinUtNXfHZiNI) 


C 

CALL 8PIT<9UH) 

0240 


CALL OUTPUT 

0249 


PRINT 1302 

0270 

c 

e 

STOP 

0271 

c 

249 

CONTINUE 

0272 


CALL BOLVEtSUHrPI) 

0273 


Nt ■ 8UN(NIrl) 4 0.01 

0274 


IJ > 0.0 

0279 


DO 18 I-l.HX 

0274 


00 21 J-lrHU 

0277 


ir(89(l> J>) 22 f 21 r 22 

0271 

22 

IJ » IJ 91 

0279 


KlfJ) > 9<1>J) 9 P»<IJ) 

0290 

21 

CONTINUE 

0291 


00 19 J-1>NX 

0292 


IF<AA(I»J)>19rl8f 19 

0293 

19 

IJ > IJ 9 1 

0284 


A(I.J) ■ A(I.J) 9 PI(IJ> 

0299 

18 

CONTINUE 

0284 

12 

CONTINUE 

0287 


60 TO 83 

0298 


RETURN 

0289 


END 
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Subrout ln« GIRL 

Doscrlptlon: Subroutine GIRL perfomo the parameter 

identification. 

Important varlablea; 

SUM Contalna the aecond gradient In lower triangular and 
diagonal locatlona and off-diagonal a priori weighting 
In upper triangular. Diagonal a priori welghtlnga are 
•tored In APR. The flrat gradient appears as an extra 
column In SUM (the JKM column) 

xji -y^) 

PHIl e^*^ 

APHI e^^ dt 

Z|U measured values of observations and controls 

XT1,XT2 computed values for observations 

XT3 variable Initial conditions on states 

XT4 variable bias on the observations other than states 

XT5 difference between estimated coefficients and the a priori 
values 

PB solution vector for the change In the estimates of the 
coefficients 

MX number of states 

MZ number of observations 

Subroutine listing; 


OOOt tUIROUTINC OISL 

r. 

0002 COdHON HSX>HA»MAHpMAT iZ>Ui02tEl >APHI pDUNiPHIt >D1 pApIp AA fiOf 

1 IJXpXJI p8UNpPepXripZEROpOS4pOD4pCpXTXpCCCp|IA8p 

2 IZEpZiXA8pICpXLA>APRpNAPRpXT4pJKMMpXT5pAPpiiP 

0003 COMMON HHp ENNpMX pMUMXpRIp IAp JK pNt p EPSpPOI p I pMUMX I t 

1 TRACE pKp Up TINE pRDSpMXPIp JKM pLMpTTpJKMMIp 

2 ALPHA pNNMl pMZp JpKJpIDpLLpLpNNpNUpMUXpENpPPp AAA p 

3 KMpFACpITRpNIpXT2pXT3pNZM1 

0004 DIMENSION AP(Sp4> pBP(8p3> pXT4(4) 

0005 DIMk^NSION XTS(25> pAPR<29) 

0006 DIMENSION Z( 7 p 3) pUi 3 p3 ) pD2< 7 ) pD04( S p 4 > pIIAS ( S) p APHI ( 3 p 4 ) p 

1 XTK7) pPHn(Sp4>pDl(8p7> pA(Sp4> pB<5p4>pAA(5p4)p 

2 ID(Sp4) pBJK2Sp4) pXJI(2Sp7)p8UM(25p 28) ppB<23> p 

3 OUN<2Sp4)pXT2<7)pZERO<S)pOS4<5p4>pXT3(7) 

C 


0007 


C0MN0N/AN8UER/ KDUOO 

0008 


COMMON/EQOATA/ ANPT 

000? 

C 

DIMENSION XL(300p7> 

0010 


COMMON/TRNSPR/XL 

0011 


COMMON/MATAB/ ALX»BLXpERX 

0012 

c . 

DIMENSION ALX<20p10>pBLX(20p10)pERX(10p10> 

0013 

c 

777 

F0RNAT(7I10> 

0014 

606 

F0RMAT(10Xp10E12.4> 

OOlS 

1011 

F0RMA7(I3p9E12.4) 

0016 

1012 

F0RMAT(12E12.4) 

0017 


ANPT-FL0AT(NNH1)41 .0 
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oott 

0019 

0020 
0021 
0022 

0023 

0024 

0025 


002* 

0020 

0029 

0030 

0031 

0032 

0033 


0034 


003S 

003* 

0037 

003S 


0039 


0040 

0041 

0042 

0043 

0044 

0045 

0046 

0047 
0046 

0049 

0050 

0051 

0052 

0053 

0054 

0055 
005* 

0057 

0058 

0059 
00*0 
0061 
0062 

0063 

0064 

0065 

0066 

0067 

0068 
00c9 

0070 

0071 

0072 

0073 


C TZHC LOOP 

TT • TINC 9 HH 
DO 41 I«2»NNN1 
TT ■ TT 9 NH 
DO 28 JK«ltJKN 
DO 21 JaHXPlfHZ 
28 XJKJKfJ) • 0.0 
DO 170 lA • UHX 

170 XJKJKNrIA) > XT2<XA) 

C 

C READ HEA8URCD DATA FROH DATA FILE 
C 

C 

C 

IFIKDUOO.EO.O.O) 00 TO 3071 
C 

REA0(4) DXr.DXY.DXV.DXV.DXYtZ(4t3) tZ<1.3)fZ(7.3)fZ(2«3) t 
1 U(li3)»U<2i3> 

Z<3f3) • 0.0 
Z<St3) > 0.0 
Z(*>3> > 0.0 
C 

00 TO 3011 
3071 CONTINUE 

C 

C LONOITUDINAL 
C 

RCAD<4> Z(4»3>>Z(lf3)>Z<7t3)fZ(*i3>»U<1.3>>DXY»DXY»DXY»DXY» 
1 DXV.DXY 

C 

c 

Z(2t3) ■ 0.0 
Z(St3) ■ 0.0 
Z(3f3) > 0.0 
U<2r3) - 0.0 
C 

C 

3011 CONTINUE 

C 

C 

DO 171 lA-l.NX 

171 Z(IA<3) > Z(IAi3>-XT4(IA> 

MAX • N1 

HA • 4 

CALL ZOT(BJI) 

JK - 0 
DO 44 J«1.HX 
DO 43 K«1>HU 

BJKJKHtJ) - 8JK JKH. J)9B<J>K)«<U(K>3>9U(Kr2) )«0.S 
IF(BB<J*K>)4Si43>45 
45 JK • JK 91 

XJZ(JKfJ9HX) - U<K»2)»BP< J4NX>K) 

BJKJK.J) >0.S«(U<Kf2)9U(K»l) ) 

43 CONTINUE 

DO 44 K ■ liHX 

IF(AA( JfK) >4*>44>46 
4* JK ■ JK 9 1 

IF(LL-1)4»4>4 
2 CONTINUE 

BJI(JKfJ)> 0.5«(Z(K.2)9Z(K.i)) 

XJI(JK>J9HX) - Z<K.2)«AP< J9HX.K) 

00 TO 44 
4 CONTINUE 

BJKJKfJ) ■ (XT2(K) 9 XT1(K))*0.5 
XJI(JKfJ9HX) - XT2(K)«AP( J9HX.K> 

44 CONTINUE 
NAX > NI 
HA > 4 
HAH - 4 
HAT « 5 

XJKNIrZ) ■ H). 

CALL HULTCXJI.PHIl.XJIfDUM) 

CALL HULT(BJI>APHI.OUH.DUH) 

CALL ADD(1.0iDUH>1.0>XJI>XJI) 
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0074 

0075 
0074 
0077 
0071 
007f 
OOtO 
0001 
0002 
00S3 

0084 

0085 
0084 

0087 

0088 

0089 

0090 

0091 

0092 

0093 

0094 
^095 
V. 596 
0097 
0090 

0099 

0100 
0101 
0102 

0103 

0104 

0105 
OlOA 
0107 


0108 

0110 

0111 

0112 

0113 

0114 

0115 

0116 
0117 
01 • 8 
0U9 
0120 
0121 
0122 

0123 

0124 

0125 

0126 

0127 

0128 


0129 

0130 

0131 

0132 

0133 

0134 

0135 

0136 

0137 

0138 

0139 

0140 

0141 


XJI(NX»2> - 
I8I4S • 0.0 
00 162 lA • l.MX 
XF<BlA8<IA))163»162tl63 
163 I8IA8 > IiIA8 6 1 
00 175 X0-1»HZ 

175 XJ1<JKHH4X8IA8>IB> - 0.0 
XJK JKHMfI8ZAS>IA) « 1.0 
162 CONTINUE 

JKNHl •• JKN >1 
DO 7 JK > IrJKHN 
DO 7 L ■ NXPlrNZ 

DO 7 K * IrHX 

XJKJKfL) • XJKJKrL) t A(L-NX »K > «X JX ( JK r K) «AP(L rK) 

7 CONTINUE 

DO 9 L -HXPlrNZ 
XJKJKM.L) « 0.0 
DO 8 K>liNU 

XJI(JKN>L) > XJI ( JKNiL)4B(L-MXiK)WU<K>3>»EiP(L>K/ 

8 CONTINUE 

DO 9 X>1»NX 

XJI(JKNrL) - XJK JKN»L)4A<L-NX>K)«XJX< JKHrK)«AP(LrK> 

9 CONTINUE 

00 3 J>1»NZ 
XTKJ) ■ XT2<J) 

XT2(J) - XJKJKNfJ) 

XJKJKH.J) > Z( J»3)>XT2( J> 

3 CONTINUE 

DO 27 K>1>HZ 

D2<K) > D2(K> 4 XJl<JKNrK>««2 
27 CONTINUE 

NAX • NI 
NA - 7 

81 CONTINUE 

C PRINTS OUT TIME HISTORIES 

C TYPE 606» <XT2(IA) »IA«lf7) »TT 

ir(LL.LT.ITR) 00 TO 80 
DO 1013 IK«lt7 
XL(XrlK) - XT2(XK> 
lOi;^ CONTINUE 

80 CONTINUE 

DO 91 Jal.JKN 

DO 91 I4>J>JKM 

DO 92 K-lfHZ 

92 SUH<I4.J) - SUM(I4>J)4XJX(14>K)«01(K.K)«XJX< J»K> 

91 CONTINUE 

DO 69 IA>1 fMZ 
Z(1Af1) - Z<IAf.2> 

69 Z(1Af 2) - Z(IAf3) 

U(lfl) - U<1>2) 

U(2rl> a U<2 f2> 

Ua>2) > U(lr3) 

U<2f 2) a U<2>3> 

41 CONTINUE 

C 

PRINT 607» SUH(JKHfJKM) 

607 F0RMAT</.10X» 'WEIGHTED ERROR SUM « '»F12.4> 

C TYPE 606t SUM(JKMfJKM) 

MAX a 8 

HA a 7 

C CALL SPIT(Dl) 

PRINT 608 

608 F0RHAT(/f10Xf 'WEIGHTED ERRORS!'./) 

PRINT 606><D2(IA)«Dl(IArlA) .lAal.HZ) 

TYPE 606F<D2(IA)»IAai.HZ> 

DO 2101 lAal.NZ 
ERX(lAfLL) a D2(IA) 

2101 CONTINUE 

DO 888 UK a l.JKM 
888 SUM(IJK.JKM) > SUN(JKMfIJK) 

IF(MAPR) 180.181.180 
180 DO 182 IBai.JKH 
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0142 

0143 

0144 

0143 

0144 

0147 

0148 
ei4» 
0150 
0131 


0001 


0002 

0003 


0004 

0005 
0004 

0007 

0008 


000? 

0010 

0011 

0012 

0013 

0014 

0015 
0014 

0017 

0018 

0019 

0020 
0021 
0022 

0023 

0024 

0025 
0024 

0027 

0028 


lUHdtfJKM) ••XT3(Ii)9A^R(I8)4fUN<Ii> JKN) 
lUHdtflt) « 8UMdldB>4APR<ll) 

!iNl - Il-l 
00 182 lA >1>I8H1 

182 SUHl.TltlA) - BUHdIdA) » 8UH<IAdO) 

181 CONTINUE 

531 FORNAT(///' END OF ITERATION '//) 

TYPE 331 

RETURN 

END 


Subroutine EAT 

Description: Subroutine EAT computes e^^'and dt 

using the Taylor series expansion to ten terms. These are 
returned as PHI], and APHIl respectively. 


subroutine listing: 

SUBROUTINE EAT < A» T >PHI » APHl > A2 » A3) 

C 

C - - 

C THIS SUBROUTINE CONPUTES THE TRANSITION HATRIX 
C AND IT'S INTEGRAL USING A TAYLOR SERIES EXPANSION 
C TO 10 TERNS. 

C 

C A > STABILITY MATRIX 

C T > DELTA TIME INCREMENT 

C PHI • TRANSITION MATRIX 

C APHI > INTEGRAL OF THE TRANSITION MATRIX 
C A2 - DUMMY MATRIX 

C A3 « DUMMY MATRIX 

C - - 

C 

COMMON MAXfMAXliMIXltMIX 

DIMENSION A(l>»PHI(l>fA2(l>rAPHIC .<>A3<1> 

C 

C CALLS MULTIPLICATION AND ADDITION SUBROUTINES 

C 

C 

HAX2 •> MAX»2 
II - A(MAX) 

JJ • A<MAX2) 

PHI(MAX) - A(MAX) 

PHKMAX2) - A(HAX) 

C 

C INITIALIZE TO ZERO AND CREATE NEW MATRICES 
C 

CALL ZOT<PHI> 

CALL HAKE(APHI>PHI) 

CALL HAKE(A3fPHI) 

MI • -MAX 
DO 1 I-ldl 
MI a MI4MAX 
Mil a HI+I 
PHKNI + I) a 1,0 

1 CONTINUE 

CALL MAKE(A2iPHI> 

0 a 1.0 
DO 2 I-ldO 
BB a I 
6 a G»T/BB 

CALL A00<l.tAPHIfG>A2>APHI) 

CALL MULT(A»A2>A2>A3) 

CALL A0D(l.>PHI»G»A2fPKI) 

2 CONTINUE 

BO 10 laldl 
DO 10 Jal,I 
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002V 

0030 

0031 

0032 

0033 

0034 

0035 
0034 


0037 

0038 


0001 


0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 
0011 
0012 
0013 


0001 


0002 

0003 

0004 
0006 
0007 

0009 

0010 
0012 

0013 

0014 


10 

C 

C 


J1 • <I-1)*HAX4J 
IJ - (J-DVHAX4I 
TEHA ■ AHKIJ) 
AHl(lJ) > ANKJI) 
PMKJI) ■ TEMP 
TEMP - APHKIJ) 
APHKIJ) > APHIUX) 
APNKJ!) > TEHP 
CAUL IPIT(PHI) 

CALL 8PIT(APHX) 
RETURN 
END 


Subroutine ZOT 


Description: Subroutine ZOT .Intializes the elements of 
a matrix to zero. 


Subroutine listing: 


SUBROUTINE ZOT(X) 

C 

e 

C THIS SUBROUTINE SETS ALL ELEMENTS OF A MATRIX 
C TO ZERO. 

C 

C X : MATRIX TO BE ZEROED 

C 

COMMON MAXiMAXlfNIXltMIX 
DIMENSION X(l> 

C 

MAX2 <* MAX » 2 
IlMl > X(MAX) -1.0 
JJMl ■ X(MAX2) -1.0 
LEND - JJMlVMAXfl 
00 1 L>1.LEND.MAX 
KENO - L IIMl 
DO 1 K«L>KEND 
1 X(K> > 0.0 

RETURN 
END 


Subroutine LOAD 


Description: Subroutine LOAD loads matrices from the 

input file. 


Subroutine listing; 

SUBROUTINE LOAD (NfA.B.C.D) 

C 

g . 

C THIS SUBROUTINE LOADS MATRICES AfBtC AND D FROM 
C AN INPUT FILE. THE VARIABLE N SPECIFIES THE NUMBER 

C OF MATRICES TO BE LOADED. 

C 

C 

REAL A(l)iB<l).C(l)fD(l) 

C 

CAl.L LOADKA) 

if<n.lt.2) return 

CALL LOADKB) 
if(n.lt.3) return 
CALL LOADKC) 

IF(N.LT.4) RETURN 
CALL LOADl(D) 

RETURN 

END 
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0002 

0003 

0004 
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0004 

0007 

0008 

0009 

0010 
0011 
0012 

0013 

0014 


0001 


0002 

0003 

0004 

0005 
0004 

0007 

0008 

0009 

0010 

oon 

0012 

0013 

0014 

0015 

0016 


Subroutine LOADl 

Description: Subroutine LOADl actually loads the matrix 

from the Input file. 


Subroutine listing: 

SUIROUTINE LOADKA) 

C ROUTINE CALLED DY LOAD LOADS NATRIX A FROM PILE 
C 

CONHON MAX 
REAL A<1) 

READ<2>100) llfJJ 
100 P0RNAT<ex»I2*I10) 

KE >(JJ>1)«NAX 
DO 10 I-1»II 
KEND - I4KE 

10 READ(2>1001) (A(K> rK«IfKENDfNAX> 

A(NAX> > II 
A<HAXX2> > JJ 
1001 P0RNAT<8P12.6) 

RETURN 

END 


Subroutine ADD 

Description: Subroutine ADD adds scalar multiples of 

two matrices, Z»g X + h Y. 

Su^'routine listing: 


SUBROUTINE ADD <0>X»HfY»Z> 

C 

c ------ 

C THIS SUBROUTINE ADDS SCALAR MULTIPLES OP TUO 
C MATRICES AS POLLOUSt 
C 

C CZ3 - G«CX3 i HXCY3 UITH t 0 •• 1.0 

C ( NO CHECKING IS HADE FOR MATRIX COMPATIBILITY ) 

C 

C 

COMMON MAXfMAXl>MIXl>MIX 
DIMENSION X(l)fY<l)>Z<l) 

C 

MAX2 - MAX X 2 
II « X(MAX) 

JJ « X(MAX2> 

JEND » ( JJ-DXMAXfl 
IlMl » II-l 
DO S3 J>l>JENDfMAX 
KEND « J-t-IIMt 
DO S3 KoJfKEND 
S3 Z(K) > X'(K)4HXY<K) 

Z(MAX> - X(MAX) 

Z(MAX2) > X(MAX2) 

RETURN 

END 
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1 


i 

i 


Subroutine MAKE 

Description: Subroutine MAKE moves a copy of the matrix 

Y into X. 


Subroutine listing: 


0001 tUIROUTXNE HAKE(X»V> 

C 

C - - - 

C THIS SUBROUTINE GENERATES A HATRIX X THAT 18 
C A COPY OF HATRIX V. 

C 

C X : NEW HATRIXt COPY OF Y 

C Y : MATRIX TO BE COPIED 

C - - - 

C 


0002 

COMMON MAXtHAXlfMIXlfNXX 

0003 

DIMENSION X(1)>Y(1) 

Q 

0004 

HAX2 

> MAX»2 

0005 

IINl 

■ Y(MAX) -1. 

0006 

JJMl 

■ Y(MAX2> -1. 

0007 

LEND 

■ JJMIBNAX 41 

0008 

DO 1 

L«1»LEND>MAX 

0009 

KEND 

- L4IIH1 

0010 

DO 1 

K>L»KEND 

0011 1 

X(K> 

• Y(K> 

0012 

X(HAX> - Y(MAX) 

0013 

X(MAX 

12) ■ Y<MAX2> 

0014 

RETURN 

0015 

END 



Subroutine MULT 

Description: Subroutine MULT computes the matrix product 
,C ■ A B. The matrix C can not be the same as matrix A or B. 


Subroutine listing: 


0001 


0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 
0011 
0012 

0013 

0014 

0015 

0016 

0017 

0018 
0019 


SUBROUTINE MULT (AfBrCfDi 
C MULTIPLIES A AND B AND PUTS THE PRODUCT 
C IN C AND D (USING SUB. MAKE) 

C 

COMMON MAXiMAXl>MXXl>MIX 
DIMENSION A(l)rB(l)>C(l>fD(l) 
C.. 

MAX2 > MAX»2 
MIX2 » MIX«2 
II > A(MAX) 

JJ - A(MAX2) 

KK > B(MIX2> 

JE - (JJ-1)»HAX 
KE ■ (KK-1)*MAX 
DO 20 1-1»II 
KEND - KEfI 
JCND - JEfl 
L ■ I 

DO 20 K-IfKEND>MAX 
D(K) > 0.0 
JB - L 

C INITIALIZATION LOOP 

DO 10 J-l7JENDfMAX 
0(K) > A( J)BB( JB)4D(K) 
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C020 

10 

JB ■ Ji ♦ 1 

0021 

20 

L > t 6 NIX 



D(NAX) • A(NAX> 

0023 


D<HAX2> - i(MIX2) 


C 

COPT D INTO C 

0024 


CACL MAKE<C»D) 

0025 


RETURN 

0026 


END 


Subroutine SPIT 

Description: Subroutine SPIT prints out a marix. 


Subroutine listing: 


0001 


SUBROUTINE 8PIT<X> 




W 

c 

c 

SUBROUTINE USED FOR THE PRINTOUT 

OF 

MATRICES 

0002 

c 

COMMON MAXpNAXliHIXltMlX 



0003 


DIMENSION X(l) 



0004 

100 

F0RNAT(13Xp' DIMENSION 'p8XiI3f' 

BY 

' *13) 

0005 

101 

F0RMAT<12Xf 10E12.4) 



0006 

W 

MAX2 - MAX»2 




0007 

II • X<MAX> 

0008 

JJ s X<MAX2) 

0009 

PRINT lOOf IIpJJ 

0010 

KE - <JJ-1)«MAX 

0011 

DO 1 I-1>II 

0012 

KEND - I4^KE 

0013 1 

PRINT lOlf (X(K> pK>I>KENDpMAX> 

0014 

RETURN 

0015 

END 


Subroutine SPITl 

Description: Subroutine SPITl prints out the A and B 

matrices with 's to show which of the parameters have 

been allowed to vary. 


Subroutine listing: 


0001 


0002 

0003 

0004 

0005 

0006 


ooov 

0008 

0009 

0010 
0011 
0012 


c 

c 

c 

c 

c 


100 

101 

c 

c 


SUBROUTINE SPITl ( X t XX > K1 ) 

SUBROUTINE USED FOR THE PRINTOUT OF MATRICES 


COMMON MAXrMAXlf MIXltMIX 
DIMENSION X(l>fXX(l> 

BYTE CHAR i 4 ) 

format: lOXf 'DIMENSION 'tI3f' BY'fI3) 
F0RMAT(10 XiS(PE12.4pA1)) 


HAX2 > MAX«2 
II > X(MAX) 

JJ - X(MAX2) 
PRINT 100> IIpJJ 
KE » <JJ-1)*MAX 
DO 1 I*lfII 
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0013 

0014 
0013 
OOU 
0011 

0019 

0020 
0021 


0001 


0002 

0003 

0004 
0003 
0006 

0007 

0008 

0009 

0010 
0011 


0012 

0013 

0014 


0013 

0016 

0017 

0018 

0019 

0020 

0021 

0022 


KCND - 14KC 

DO 2 KaltKENDtHAX 

CHA(I(<K-1>/HAX91)>' ' 

IF(XX<K).EQ.O. )CHAK<<K>I>/8AX^l>-'«' 

2 CONTINUE 

1 PRINT 101»((X(K)»CHAR<<K>I>/HAX41>)fK*IfKEN0fNAX) 

RETURN 
END 


Subroutine SOLVE 

Description: Subroutine SOLVE solves the system of linear 

equations, A x ■ b where A is symmetrical. Only the lower 
triangular and diagonal elements of A are used. The b vector 
Is assumed to be stored In the N+1 column of A, where N Is the 
dimension of the system. 


Subroutine listing: 

SUBROUTINE SOLVE(ArX) 

c 

C SOLVES SYSTEM AX • B UHERE A SYMMETRIC MATRIX 
C AND B A MATRIX IN N91 COLUMN OF A 
C 

REAL A<25>23) >X<23> 

CALL REDUCE(A) 

N - A<25»1) 

NMl > N-1 
NPl > N41 

C MULTIPLY MATRICES! (L> « (S> . . . 

DO 70 I-2!N 
X(I) - AdfNPl) 

IMl - I-l 
DO 70 J-lflHl 

70 X(I) - X(I)fA(l!j)«A(j!NPl> 

C MULTIPLY BY iDI) 

C 

A(l!NPl) « A<l!NPl>/A(l.l> 

DO 80 I-2fN 

80 AdtNPl) > Xd)/Ad!l) 

C MULTIPLY BY (LS) TO FORM < L« > « ( DI > » ( L ) « ( B > 

C 

DO 90 I-l.NMl 
Xd) - AdfNPl) 

IPl ■ I+l 
DO 90 J*XPlfN 

90 Xd) > Xd)4A( Jf I)«A( JfNPl) 

X(N) - A(NfNPl) 

C 

RETURN 

END 
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Subroutine DIAGIN 

Description: Subroutine DIAGIN obtains the diagonal 

elements of the inverse of a symmetric matrix. 

Subroutine listing: 


SUiaOUTINC DIAOXNCA) 

c 

C FIND DIAGONAL CLEHENTI OF A INVERSE FOR 8YHNETRIC A 
C 

REAL A<28f2S> 

CALL REDUCE (A) 

N • A(2Stl) 

NNl « N-1 
DO 90 I-lfNNl 
A(I»I) ■ 1.0/A(ItI> 

IPl • 141 
DO 90 J>IPlrN 

90 A<IfI) -A<I>I)4A<JtI>««2/A(Jf J) 

A(N>N> ■ 1.0/A(N>N) 

RETURN 

END 


Subroutine REDUCE 

Description: Subroutine REDUCE factors a symmetric 

matrix A by Cholesky's matrix decomposition. 


Subroutine listing: 

SUBROUTINE REDUCE: A) 

C REDUCES STMHETRIC HATRIX A STORED IN LONER TRIANGULAR LOCATIONS 
C TO THE FORN ( LI ) S ( D) «<L1« > UHERE L IS A LONER TRINOULAR MATRIX 
C NITH UNITY DIAGONAL TERMS* D IS A DIGONAL MATRIX* I DENOTES 
C INVERSE AND « TRANSPOSE. 

C 

C - 

C 

REAL A (25 *25) 

N ■ A(2S*1> 

NMl * N-1 
C MAIN DO LOOP 

DO 20 K-1*NH1 

KPl ■ K+1 

KMl - K-1 

AKKI > 1.0/A(K*K> 

C 

DO 20 I>KP1*N 
AKKIK « A(I*K>«AKKI 
DO 10 J-I*N 

10 A(J>I) * A(J*I)-AKK1K»A(J*K> 

A(I*K) «-AKKIK 

IF(KM1 .EQ.0.0) GO TO 20 

DO IS J«1*KM1 

IS A(ItJ) > A(I*J) -AKKIKSA(K*j;i 

20 CONTINUE 

C 

C (L) IS NON STORED IN LONER TRIANGULAR PART OF (A) 

C EXCEPT OF MATRIX DIAGONAL* NHICH CONTAINS 0 
RETURN 
END 
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Subroutine CRAMER 

Description: Subroutine CRAMER computes tUe confidence 

levels of the estimated derivatives. 


Subroutine listing: 

SUiaOUTINE CRAHER(HU>HX>HZ>N1) 

C 

C 

c 

C THIS SUBROUTINE CONFUTES THE CRAHER-RAO BOUNDS 

C ALSO KNOWN AS THE CONFIDENCE LEVELS OF THE 

C E8TINATED DERIVATIVES. 

C NU * NUHBER OF CONTROL INPUTS 

C NX - NUNBER OF STATES 

C NZ - NUNBER OF OBSERVATIONS 

C NI « NAX. NUNBER OF UNKNOWNS (2S> 

C 

c ... 

C 

CONNON NAX>NA>NAN>NAT»ZtU>D2rElf APHIf DUNiPHIl»Dl»A>BfAAtB6f 

1 BJI>XJI»SUNtPBfXTlfZER0»DS4fDD4>E»XlXrCCCrL:AS» 

2 IZEfIBIASiICf XLA»APRfNAPRtXT4» JKNN>XT5»APf BP 
CONNON /EQDATA/ ANPT 

DINENSION AC<5»4> f BC<5f 4) 

OINEHSION XTS<2S>>APR(2S) 

DINENSION AP<e»4>>BP<8*3)fXT4(4> 

DINENSION Z(7>3) fU(3f 3) rD2<7)f DD4<5>4> >BIA8(5) rAPHI(Sr4)f 

1 XTl(7>fPHIl(5r4)>Dl(e»7)fA(5«4)>B(Sf4)tAA(5>4)« 

2 BB(5f 4)f 8JI(25>4)>XJI(25t7)f 8UN(2St2S)f PB(2S> r 

3 DUN(25>4) >XT2(7) fZER0(5) •C54(5f 4) f XT3<7) 

C 

C NORNALLY THE APRIORI CONTRIBUTION TO HESSIAN IS SUBTRACTED 
C FOR THIS CONFUTATION BUT THIS ROUTINE ASSUNES NO APRIORI 
C OPTIONS ARE BEING USED AND HENSE THERE ARE NO CONTRIBUTIONS 
C 

AC(Sfl) > NX 
AC(S>2) > NX 
BC(Sfl) > NX 
BC(Sr2) ■ NU 
JKNNl > SUN(NIfl) -f 1.01 
C 

C STORE WEIGHTED ERROR SUN IN ERRSUN 
ERRSUN ■ SUN( JKNNlfJKNND/ANPT 
C 

C OBTAIN DIAGONAL ELENENTS OF INVERSE 

C 

CALL DIAGIN<SUN) 

C 

C CONFUTE CRANER-RAO BOUNDS 

C 

WTS >0.0 
DO 1 I»lfNZ 

IFdiKItl) .NE.O.O) WTS > WTS 1.0 
1 CONTINUE 

COE. F • ERRSUN/WTS 
PRINT lOf ERRSUNtCOEFFf WTS 

10 FORNAT(' ERRSUN « '*F12.4t' COEFF ■ '»F12.4.' WTS ■ '.F12.4) 

L » 0 

DO 2 I«1.NX 
DO 3 J-ltNU 
BC(I>J) * 0.0 
1F(BB(1>J).NE.1.) GO TO 3 
L ■ L+1 

BC(IfJ) > SORT(ABS(SUN<LfL> )«COEFF> 

3 CONTINUE 
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0011 

0012 

0013 


0014 

0015 

0016 

0017 

0018 

0019 

0020 
0021 
0022 

0023 

0024 

0025 

0026 

0027 

0028 

0029 

0030 


DO 4 J-1fHX 
ACUfJ) - 0 

IF(AA<I»J).NC.l.) GO TO 4 
L ■ L41 

4C(I>J) ■ t0RT<AB8(8UH<L»L>)*C0eFF) 
4 CONTINUE 

2 CONTINUE 

NAX • 5 
FRINT 4 

4 FORHAT(' AC MATRIX') 

CALL 8RIT<AC> 

FRINT 7 

7 FORMAT <' DC MATRIX') 

CALL SPIT (DC) 

RETURN 

END 


Subroutine OUTPUT 

Description: Subroutine OUTPUT provides the output of 
time histories and matrices to user defined files for later 
plotting. 


Subroutine listing: 

8UDR0UTINE OUTPUT 

C - ------------- 

C THIS 8UDR0UTINE UlLL PROVIDE MMLE RESULTS 

C IN A FILE TO BE SPECIFIED BY THE USER. 

C THE FILE WILL CONTAIN INFORMATIOM ABOUT 

C THE MATRICES A AND D AT EACH ITERATION. 

C - 

C 

COMMON/MATAB/ ALXvBLX 
COMMON/TRNSFR/ XL 
COMMON/CONST/ ITRtNN 
BYTE INAHE(IS) 

DIMENSION ALX(20f 10)»BLX(20rl0) 

DIMENSION XL(300>7) 

C 

C 

TYPE 10 

10 FORMAT</' ENTER FILE NAME FOR OUTPUT OF MMLE MATRICES'/* '»' ) 
ACCEPT ll*(INAME<IAB)rIAB>l*14) 

11 F0RMAT(14A1) 

INAME(15)-0 

0PEN(UN1T»2*NAME-INAME*TYPE«'NEU' fACCESS-'SEQUENTIAL' * 

1 FORN-'FORMATTED' *BUFFERC0UNT«2) 

C . 

35 FORMAT( 10E12.4) 

40 FORMAT(' $«« MATRIX A «««'/) 

50 FORMAT(' ««» MATRIX B *»$ '/) 

60 FORMAT(' ITERATION ... '*I2>/) 

70 FORMAT(//.' ESTIMATED TIME RESPONSES'//) 

DO 130 1>1*ITR 
URITE(2f60) I 
URITE(2t40) 

00 30 J"ltl4>4 

URITE(2t35) (ALX(J-14K>l)>K>li4) 

30 CONTINUE 

HRITE(2*50) 

DO 131 J>1»14*4 

HRITE(2f35) (BLX( J-14^K* I ) >K-lt3) 

131 CONTINUE 

130 CONTINUE 

C 

CALL CLOSE (2) 
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0031 

C 

TY*E to 

0032 

10 

rO«M4T</' TY« IN FILE THAT MILL CONTAIN'. 

1 ' CAtT INERATION TINE KEtFONtEI. ' /. 't' ) 

0033 


ACCCFT ll.<INANE(IAl>.IAt-1.14> 

0034 

c 

OFEN ( UNI T>3 . NANE> INANE . T YFE- ' NEU ' . ACCEtl> ' lEQUENT I AL ' > 
1 FORN- ' UNFORNATTED ' » 1UFFERC0UNT>2 > 

0033 

c 

DO 330 N>2.NN 

0034 


URITE(3) (XL<N.I).I>1.7) 

0037 

330 

CONTINUE 

0031 


WRITE(3) 

003* 


RETURN 

0040 


END 
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A. 7. 3) tflLE OUTPUT FORMAT 

Following Ic oxnaple and dascrlptlon of tha MMLE oucpuc. 
Longitudinal: 


KU rRL lONII HHLC miULTI 
CCttNA 172 tONtITUaiNAL CAff 
3000. rr ALT. AT 17A. 7M AlAOACfO 
PLIOHT If/tO/tO AUN 23 


INITIAL CONBITIONI 

NUHBKR or OATA POINTO t 240 NAXINUN NUHBIA or ITINATlONf I f 

DATA lAMPLINO INTCAVAL t O.tOOO FIAOT DATA POINT AT TINK I 0.0000 

OIAOONAL NULTIPLTINO PACTOA I l.OOOO NUNOtA OP OTATIt t 7 


ZCAOt AND DIAOCI 

0.00001400 O.OOOOC400 0.00001400 0.0000C400 

O.0000C40O O.0O0OC4OO 0.OOOOC4OO O.OOOOC400 
DIAOONAL CLCNCNTS OP TNC NIIOHTINO NATAIX DU 
100.000 0.000 0.000 700.000 0.000 40.000 S.OOO 


INITIAL INPUT HATAICCa CA} AND CDl. 

A ITAA <») FOLLOUINO TNC VALUC OF A HATAIX 
ELCnCNT INOICATCI THAT THE ACIPECTIVC DCAIVATIVE 
If NOT CITINATCD »V THE MNLC HCrHOD. 


OTAIILITY HATI 
DIHENIION 4 
-0.SIA3E401 
O.OOOOE4004 
0.1000E4014 
0.f»PSC4004 


IX CAl 
DT 4 

O.OOOOE4004 

0.0000c400> 

O.OOOOE4004 

O.OOOOE4004 


•0.2130E402 

0.1f34C402 

-0.299SE401 

O.0OOOE4OO4 


O.OOOOE4004 

-4.32UC402k 

-0.4SA0E-024 

0.0000C4O04 


1 


CONTAOL HATAIX CD] 
DIHENIION 4 DT 3 
-0.3043E402 O.OOOOE400k 
-0.S91SE401 O.OOOOE4004 
-0.2409E400 O.OOOOE400* 
O.O0OOE40O* O.OO00E4OO* 


O.OOOOE400 

0.000OC4O0 

0.0000E400 

O.OOOOE400 


2 


UEIONTED ERAOA SUN • 2197421.2300 


UCIOHTED CAAOABt 

3 

0.14I1E404 O.OOOOC400 0.00OOC400 0.2193E407 0.0000C400 0.1147C402 0.290SC404 

ITERATION 1 UAS COHPLETED 


■ 



Mi 

M* 

M* 

<1 

u 

a 

0 

z + u. 

X* 

u 

X' 

a 

q 1 


f • 

u, - z- 
1 a 

z 

u 

z 

a 

COB (Oj) 

• 

0 

0 


*(q,U,a,e,A^,Ajj) 


-f coa (Oj^) 

^ " *■ 2: 

1 o 

0 



0 0 i' 

o 
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Lattral 


KU r»L •ONfl NNLI RIIUI.T* 

CCtIN* 173 LATCML-tmCCTtONAL CMC 
3000. fT. ACT. AT 17A. 7P0 AIRIPCCO 
FLIONT If/tO/iO RUN 4S 


INITIAL CONIITIONR 

NUNIIR or DATA POINT! t 140 HAXIHUH NUNOCR OP ITIRATIONt i f 

DATA lANPLINO INTCRVAL I 0.1000 PIRIT DATA POINT AT TINC I 0.0000 

DlAflONAL NULTIPLTINO FACTOR t 1.0000 NUNOCR OF ITATCi i 7 


ZCROI AND DIAOC! 

0.0O00C40O 0.0000C40O O.0000C4O0 0.0000C400 

0.00O0C40O 0.0000C400 O.OOOOC400 0.00O0C40O 

DIAGONAL CLCHCNTt OF THC UCIOHTINO NATRIX D1 i 

40.500 40.000 0.000 150.000 0.000 0.000 100.000 


initial input NATRICCC CA3 AND CD). 

A iTAR (•) FOLLOMINO THC VALUC OP A NATRIX 
CLCNCNT INDICATI! THAT THC RCIPCCTIVC OCRIVATIVC 
II NOT crtinatcd it the NHLC HCTHOD. 


ITADILITY NAT 
DINCNSION 4 
-0.»74fl+01 
-0.33I4C400 
0.3943C*01« 
O.lOOOCTOlt 


IX CA) 

IT 4 

0.ir*^E401 

-0.I117C401 

-0.9f93C400« 

0.3744C-01I 


-0.1I1SC403 

0.7214Ev01 

-0.1737C400 

O.OOOOC400* 


0.0000C40O> 

O.OOOOE400* 

0.1I27C400* 

0.00O0C40O4 


I 


CONTROL NATRIX Cl) 
DIHENIION 4 IT 3 

0.3430C402 0.299Ivl.^01 

-0.SII2E401 -0.739IC401 

O.OOOOC400 0.37S0C4O1 
O.OOOOC4004 0.000OE40O4 


0.0000C400 

0.0000E400 

O.0OOOC4OO 

O.OOOOE400 


2 


UEIOHTED ERROR lUN ■ 52015.3211 

UEIOHTED ERRORIi 

0.23I7C403 0.I472E403 O.OOOOEfOO 0.123fEi05 0.0000C400 O.OOOOE400 0.3IS2E405 

ITERATION 1 HAI COMPLETED 


1 

• 

K 



« 

0.0 

2 





% 

M' 

r 

% 

0.0 




K 


aln(aj^) 

-coa(a,) 


^ C0«(ej^)C0R(4j) 



r 

n 


1.0 

coR(4j^)Can(»j^) 

0.0 

0.0 

m 


0.0 

0.0 

♦o. 
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A.8) TIME HISTORY PLOTTING 


PRECEDING PAGE BI.ANK NOT RLMED 


Dcecriptlon: Th« PL0T03 program collacta tha time hlatories 

of both the actual inflight measurementa» and the predicted 
states of the MMLE BONES program. These are then plotted on 
the graphics CRT terminal. This process allows the user to 
observe the fit to the measured states. For hard copy plots 
the time histories are transfered to the KU~FRL’s Hewlett 
Packard computer. 

Program listing: 


coot 

0002 

0003 

0004 

0005 
0004 

0007 

0008 

0009 

0010 


0011 

0012 


00 3 
0014 


0015 

0016 


C t 

c « 
c t 
c « 
c s 
c t 
c » 
c « 
c « 
c s 
c 


c 

c 

c 

1 


c 

c 

c 

301 


THIS PROORAN CAN SC U8CD TO PLOT « 
HCA8URCD AND E8TIHATED DATA. * 
DATE 18-N0V-80 * 
THE REQUIRED lUSROUTINES ARES « 
PL0T5SiINIT AND GRAPH. LINK AS « 
FOLLOWS t PL0T«HAIH.PL0T55.ir’T. * 
GRAPH (CR>. * 


C0HH0N/STATU8/ISTAT<16> 

DIHCN8I0N IARRAY<512) >HH<9f n fRDATA(240.9> »CDATA(240f9) 
DIliivNSION 0AIN(9) »GAIN1(9) »0AIN2(9) 

CONHON/FIUEA/ DATAfOAIN 

lYTE YC8fN0rAN81.L00»DIR>AN82»NAH£(15) 

DATA LOGfDIR /'L'»'D'/ 

DATA YESfNO /'Y'»'M'/ 

DATA I8TAT/1640/ 

DATA 0AINl/114.6f 114.6f 143.2>47.7>S7.3>57.3» 100 . » 1 43 . 2 r 143. 2/ 
DATA GAIN2/114,6f 2.5>143.2.95.S>S7.3f 100. f 25. f 143.2>0./ 


TYPE 1 

FORHAK/////' MELLOOOOOO. . './///. 

I' THIS IS A PLOTTING PROGRAH USED TO PLOT'/. 

2' HEASURED VERSUS ESTIMATED TIME HISTORIES'//. 

3' instructions: '/. 

4' I, TYPE IN TYPE OF MANOEUVRE.'/. 

5' 2. TYPE IN NUMBER OF TIME POINTS N (0-250>.'/. 

6' 3. INDICATE NAME OF FILE CONTAINING MEASURED DATA.'/. 

7' 4. INDICATE NAME OF FILE CONTAINING ESTIMATED DATA.'/. 

8' 5. WHEN DATA IS PLOTTED HIT 'CR* TO CONTINUE.'/. 

9' 6. GOOD LUCK ALEX'//) 


TYPE 301 

format:' indicate type of manoeuvre.*/ 


1 .' *L* FOR LONGITUDINAL'/ 

2 .' *D* FOR LATERAL DIRECTIONAL'/) 
RCAD<S.302) ANSI 

302 FORHAT(Al) 

C 

C 

C 

C READ IN NUMBER OF TIME POINTS 
C 
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0017 

10 

r0ftHAT<I3) 

OOlt 

20 

rOKNATi' TYPE !N NUH8CR OF TIME POINT'/) 

OOlf 


TYPE 20 

0020 


READ(S»10> N 


W 

c 

f* 

READ DATA PILE WITH HEAIURCD DATA 

0021 

C 

TYPE 30 

0022 

30 

PORMAT(' TYPE IN NAME OF DATA FILE UITH MEASURED DAlA'/) 

00/3 

c 

NAME(IS) ■ 0 

00', 4 

31 

FORMAT( 14A1 > 

0025 


ACCEPT 31f (NAME(I)>I«1>14) 

0026 


0PEN(UNIT-2>NAME-NAME*TYPE*'QLLI' fACCESS-'BEQUENTIAL' r 
1 READONLY>FORM*'UNFORMATTED' > 

0027 


DO 32 I>lfl4 

0028 


NAME < I) ■' ' 

0029 

32 

C 

CONTINUE 

0030 

C 

40 

FORMAT< 12E12.4) 

0031 


DO 50 I>1>N 

0032 

C 

IF(AN61.EQ.L00) GO TO 310 


c 

c 


0034 


REA0<2) OHY>OMY»OMYrDMYf DMY»HH<4 p 1 ) »HH(1>1 ) »HH(7*1) r 
1 HH(2>l)rHH<8f U>HH<9f 1) 

0035 


HH(3>1) « 0.0 

0036 


HH(S*1) ■ 0.0 

0037 


HH<6>n - 0.0 

0038 


GO TO 311 

0039 

310 

CONTINUE 

0040 


READ(2> HH(4>1) iHH(l >1 ) >HH(7f 1 > tHH( 6 1 1 > rHH( 8 » 1) fDNY> 


C 

1 DMYf DNYvONYf DMYfDNY 

0041 

C 

HH(2>1) > 0.0 

0042 


HH<3»1> - 0.0 

0043 


HH<5fl) - 0.0 

0044 

C 

HH<9»1) • 0.0 

0045 

C 

311 

CONTINUE 

0046 


DO 50 J»l»9 

0047 


RDATA(IfJ) ■ HK(Jft) 

0048 

SO 

C 

CONTINUE 


C 

C 

C 


0049 


TYPE 60 

0050 

60 

FORMAT!' TYPE IN NAME OF DATA FILE WITH ESTIMATED DATA'/) 

0051 


ACCEPT 31 • (NAHE(I)»I>1>14) 

0052 


0PEN(UNlT>3>NAME-NAMEtTYPE-'0LD' >ACCESS>'SEQU£NTIAL' f 
1 READONLYrFORM-'UNFORMATTED' ) 

0053 


DC 81 K«ltl4 

0034 


NAME! I) ■' ' 

0055 

81 

CONTINUE 

0056 

70 

F0RMAT(12X>9E12.4) 

0057 


DO 80 K«lfN-2 

0058 


READ(Z) (HH(Lfl)>L»lf7) 

0059 


DO 80 J-lr7 

0060 


EDATA(K62fJ) > HH(Jtl) 

0061 

80 

C 

CONTINUE 

0062 

C 

DO 399 J»li7 

0063 


EOATAdfJ) ■ 0.0 
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0044 

0045 

0044 

0047 

0041 

004f 

0070 

0071 

0072 


0073 

0074 

0075 

0077 

0078 


0079 

0080 
0082 
0083 


0084 

0085 
0084 

0087 

0088 

0089 

0090 

0091 

0092 

0094 

0095 
0094 

0097 

0098 
0100 
0101 
0102 
0103 


EDATA(2»J) • 0.0 
399 CONTINUE 
C 

DO 500 I>ltN 
C0AT4(lr8>«0.0 
CDATA(I»9)>C.O 
500 CONTINUE 

lARRAY(^) » (V 
IARRAY(3> - 0 
113 CONTINUE 

C _ 

C 

TYPE 100 

100 FORNAT<' « » t SELECT UARIADLE8 » « «'///) 

C LONG I TUC I ANAL CARE 

IF(AN81.EQ.DIR) 00 TO 102 
TYPE 90 

90 FORMAT</.' VARIABLE RANGE 

1' 1. PITCH RATE - 50 DEO/SEC' /» 

2' 2. AIRSPEED - 20 FT/3EC'/» 

3' 3. ANGLE OF ATTACK - 20 DEO '/f 

4' 4. PITCH ATTITUDE - 30 DEO'/t 

5' 5. PITCH RATE ACCEL. - 50 DEG/8EC442 ' / . 

6' 4. longitudinal ACCEL. - .5 G 
7' 7. NORMAL ACCEL. - 2 0 Vr 

8' 8. elevator PSN. - 20 DEO '/> 

9' 9. * * » (BLANK)'//) 

C 

C LATERAL DIRECTIONAL 

C 

102 CONTINUE 
IF(AN81.EQ.L00) GO TO 103 
TYPE 91 

91 FORHATC/' VARIABLE RANGE +/-'//. 

1' 1. ROLL KATE - 25 DEG/SEC'/. 

2' 2. YAU RATE - 25 DEG/SEC'/. 

3' 3. SIDESLIP ANGLE - 20 DEO'/» 

4' 4. BANK ANGLE - 40 DEG'/. 

5' 5. ROLL RATE ACCEL. - SO DEQ/SEC««2 ' / . 

4' 4. YAU RATE ACCEL. - SO DE0/SECX42 ' / . 

7' 7. LONGITUDINAL ACCEL. - .5 Q'/r 
8' 8. AILERON DEFLECTION - 20 DEG'/t 
9' 9. RUDDER DEFLECTION - 20 DEG'//) 

C 

C 

C 

103 CONTINUE 
C 

TYPE 11 

11 FORHATC' INDICATE VARIABLE NO. FOR TOP PLOT'/) 
ACCEPT «.KT 

TYPE 12 

12 FORMAT(' INDICATE VARIABLE NO. FOR BOTTOM PLOT'/) 
ACCEPT «>KB 

105 CONTINUE 

IF<ANS1 .'EG.LOO) GO TO 411 
DO 421 1>1<9 
OAIN(I) •> GAINl(I) 

421 CONTINUE 
411 CONTINUE 

IF<ANS1.EQ.DIR) 30 TO 402 
DO 422 I-1.9 
GAIN(I) « 6A1N2(1) 

422 CONTINUE 
402 CONTINUE 
C 

C 


167 


c 

c 

c 

0104 

0105 
0104 

0107 

0108 110 

0109 

0110 

C 

C 

C 

C 

c 

0111 

0112 

0113 

0114 130 

0115 

0116 

0117 

0118 

0119 

0120 
0121 
0122 

0123 140 

0124 

0125 

0126 

C 

C 

012 ? 

0128 

0129 

C 

C 

C 

0130 

C 

0132 

C 

C 

C 

0133 

0134 
0136 
0138 
0140 
0142 
0144 
0146 
0148 

C 

0150 610 

0151 

C 

0152 

0153 
0155 
0157 
0159 
0161 
0163 
0165 
0167 

C 

0169 601 

0170 


CLEAR CRT AND PORN ORID FOR RLOTTINO 
CALL INIT 

CALL PL0T55( 2» 1424324444 128 ISTAT) 

DO no K>l>235f50 

CALL PL0T55(4»liK-ni8TAT> 

CONTINUE 

CALL PL0T5S<4rn229f I8TAT) 

CALL PL0T55<5*0nf I8TAT) 


FORN THE MEASURED AND ESTIMATED RESPN8ES 
DO 130 I«lfN 

IARRAT(2«I) • RDATA( IfKT)«0AIN(KT)4150 
1ARRAV(2«I-1) > R0ATA<IfK8)*GAIN<KB>450 
CONTINUE 

CALL PL0r55<9i20f2iISTAT) 

CALL PL0TS5<12f t'8 « « TINE HISTORIES » « « '>ISTAT> 
CALL PL0T55(9>50>4iI8TAT) 

CALL PL0TS5(12t>' MEASURED DATA '>I8TAT) 

CALL 0RAPH(2«NtIARRAY> 

DO 140 I-lfN 

IARRAY(2«I) > EDATA(nKT)«8AIN(KT)4150 
IARRAY(2«I-1> ■ CDATA(lrKB)«0AIN<KB)450 
CONTINUE 

CALL PL0TS5(9*50>4f ISTAT) 

CALL PL0T55(12>f' ESTIMATED DATA'iISTAT) 

CALL QRAPH(2«N»IARRAY> 


CALL PL0T55<9»50»4f I3TAT) 

CALL PL0T55<12.r' 'fISTAT) 

KFLAG > 1 


IF(ANS1 .EQ.DIR) 00 TO 699 

KFLA61 - 0 

LONOITUDINAL LABELS 

CALL PL0T5S<9f50>6iXSTAT) 
IF(KT.EQ.l) 00 TO 601 
IF(KT.EQ.2) 00 TO 602 
IF(KT.EQ.3> O'! TO 603 
IF(KT.EQ.4) 00 TO 604 
IFCKT.EQ.S) 00 TO 605 
IF(KT.EQ.6) OC TO 606 
IF(KT.EQ.7) 60 TO 607 
IF(KT.EQ.O) 00 TO 608 


COt'^TINUE 
KFLAGl - 1 

CALL PL0T55<9»50f 16tISTAT) 


IF 

(KB 

.EQ. 

1) 

(*0 

TO 

601 

IF 

(KG 

.EQ. 

2) 

00 

TO 

602 

IF 

(KB 

.EC. 

3) 

00 

TO 

603 

IF 

(KB 

.EU. 

4) 

00 

TO 

604 

IF 

(KB 

.EQ. 

5) 

00 

TO 

60S 

IF 

(KB 

.EQ. 

6> 

00 

T9 

606 

IF 

(KB 

.EQ. 

7) 

00 

TO 

607 

IF 

(KB 

.EQ. 

8> 

00 

TO 

608 

CALL 

PL0T55I 

(12i 

1 1 ' 

Q 

GO 

TO 

640 





25 DEO/SEC 'fISTAT) 
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0171 

0172 

602 

CALL ALQT99(12pp' 
00 TO 640 

V ♦/- 20 FEET/8EC 

0173 

0174 

603 

CALL PL0T95<12pp' 
00 TO 640 

ALPHA 4/- 20 DEO 

0179 

0176 

604 

CALL fL0T59<12pp' 
00 TO 640 

THETA 4/- 30 DEO 

0177 

0178 

609 

CALL PL0T99(12pp' 
00 TO 640 

Q DOT 4/- 50 DE0/SEC««2 

0178 

0180 

606 

CALL PL0T39<12pp' 
00 TO 640 

AX 4/- .5 0 

0181 

0182 

607 

CALL PL0T99(12pp' 
00 TO 640 

AN 4/- 2 6 

0183 

608 

CALL PL0T59<12pp' 

DE 4/- 20 DEO 

0184 

640 

CONTINUE 


0189 


IF(KFLAOl.EQ.O) 00 TO 610 

0187 

C 

00 TO 751 

• 

0188 

c 

699 

C 

CONTINUE 


0189 

C 

c 

c 

LATERAL DIRECTIONAL LABELS 
KFLA02 - 0 

0190 


CALL PL0T55<9p50p< 

>pISTAT) 

0191 


IF(KT.EQ.l) 00 TO 

701 

0193 


IF(KT.E0.2) 00 TO 

702 

0199 


IF(KT.EQ.3) 00 TO 

703 

0197 


IF(KT.EQ.4) 00 TO 

704 

0199 


1F(KT.EQ.5) 00 TO 

705 

0201 


IF(KT.E0.6) 00 TO 

706 

0203 


IF<KT.EQ.7) GO TO 

707 

0209 


IF(KT.EQ.8) 00 TO 

708 

0207 


IF(KT.EQ.9) 00 TO 

709 

0209 

0210 

w 

710 

C 

CONTINUE 
KFLA02 ' 1 


0211 


CALL PL0T55(9p50p16pISTAT) 

0212 


IF(KB.EQ.l) 00 TO 

701 

0214 


IF(KB.EQ.2) 00 TO 

702 

0216 


IF(KB.EQ.3) GO TO 

703 

0213 


IF<KB.EQ.4) 60 TO 

704 

0220 


IF(KB.EQ.S) 60 TO 

705 

0222 


IFCKB.EQ.6) GO TO 

706 

0224 


IF(KB.EQ.7) 00 TO 

707 

0226 


IF(KB.EQ.a) 00 TO 

708 

0223 

c 

r 

IF(KB.E0.9) 00 TO 

709 

0230 

0231 

c 

701 

CALL PL0T55(12pp' 
00 TO 740 

P 4/- 25 DEO/SEC 

0232 

0233 

702 

CALL PL0T55(12pp' 
GO TO 740 

R 4/- 2S DEO/SEC 

0234 

0235 

703 

CALL PL0T5S<12pp' 
00 TO 740 

BETA 4/- 20 DEG 

0236 

0237 

704 

CALL PL0T55<12pp' 
GO TO 740 

PHI 4/- 60 DEO 

0238 

0239 

705 

CALL PL0T55<12pp' 
60 TO 740 

P DOT 4/- 50 D£6/SEC««2 

0240 

024'. 

706 

CALL PL0735<12pp' 
00 TO 740 

R DOT 4/- 50 DEG/SEC»»2 

024)! 

0243 

707 

CALL PL0T5S(12pp' 
00 TO 740 

AY 4/- .5 0 

0244 

0245 

708 

CALL Pl.0T55< 12 p p ' 
00 TO 740 

DA 4/> 20 DEG 

0246 

709 

CALL PL0T55<12pp' 

DR 4/- 20 DEO 


' »I8TAT) 

'f I8TAT) 

' pISTAT) 

pISTAT> 

pISTAT) 

pISTAT) 

pISTAT) 


' pISTAT) 
'pISTAT) 
'pISTAT) 
'pISTAT) 
'pISTAT) 
'pISTAT) 
'pISTAT) 
' p ISTAT) 
' pISTIT) 


169 


0247 

0248 


0280 

0281 

0282 

0283 


0284 

0289 

0284 

0287 

0288 
0240 

0242 

0243 
0248 

0247 

0248 

0249 

0270 

0271 

0272 

0273 

0274 
0274 

0277 

0278 

0279 

0280 
0281 

0283 

0284 


0288 

0284 

0287 

0288 


0001 

0002 

0003 
OOC4 
0008 

0004 

0007 

0008 


0001 

0002 

0003 

0004 
0008 
0004 

0007 

0008 

0009 

0010 
0011 


c 

c 

c 

740 CONTINUE 

IF<KFLA02.EQ.0) 00 TO 710 
C 

C LOGIC FOR OENERATINO NEW PLOTS END TERHINATINO 
781 CONTINUE 

READ(Sil80) KR 
180 F0RNAT<I2) 

CALL INIT 
C 
C 

TYPE 210 

210 FORNATl' DO VOU WANT TO REPLACE TOP PLOTT (Y/N)'/) 

READ(8f220> ANS2 
220 FORMAT (Al> 

IF(AN82.EQ.N0) KFLAO » 0 
IF(ANS2.EQ.N0) 00 TO 230 
C 

TYPE 100 

IF<AN81.EQ.L00> TYPE 90 

IF(ANSl.EQ.DIR) TYPE 91 

TYPE 240 

240 FORMAT<' INDICATE NEU OARIADLE NUMBER '/> 

ACCEPT *»KT 
230 CONTINUE 
TYPE 309 

309 FORMAr<' DO YOU WANT TO REPLACE BOTTOM PLOTT (Y/N)'/) 

READ(Sf220) AN82 
IF(ANS2.EQ.N0> 00 TO 400 
TYPE 320 

320 FORMAT(' INDICATE NEU UARIABLE NUMBER'/) 

ACCEPT «fKB 
60 TO 108 
400 CONTINUE 

IF(KFLA6.E0.0. > 60 TO 410 
00 TO 108 
410 CONTINUE 
C 

C 

C 

CALL PL0T88(2>812Flf244T321^44f ISTAT) 

CALL PL0T88(0>-lr0f ISTAT) 

RETURN 

END 


SUBROUTINE INIT 
C0MM0N/STATUS/ISTAT(14) 

DATA 1STAT/14X0/ 

CALL PL0TS8(13.72»f IS7AT) 
CALL PL0TSS(13r74»»ISTAT) 
CALL PL0TSS<2»l+Sl2»tISTAT) 
RETURN 
END 


SUBROUTINE ORAPH(Nt lARRAY > 

C0MH0N/STATUS/I3TAT( 14) 

DIMENSION IARRAY(812> 

NUNBER-ISTAT<8)/8 
CALL PLOT58(7fO»Of ISTAT) 

CALL PL0TS8<8fS12f0tICTAT) 

CALL PL0TS8<2i 14(NUMBERT1 )*2r (NUMBERTl )«iOr ISTAT) 
CALL PL0TS3<3»-Nf IARRAV fISTAT) 

CALL PL0T88<lf l-NUNBERrrlSTAT> 

CALL PL0TSS(9.10rl. ISTAT) 

END 
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A. 9) CESSNA PROGrUMS 


This appendix contains listings of the programs used In the 
Cessna spin test program. 

A. 9.1) DATA ACQUISITION 


Description: This Is an assembly language program for the 

AIM 65. The program Is essentially the same as the one of Appendix 
A.l, the differences being 

- no start and end data Is taken; 

- Channels 0-14 are sampled continuously every 0.1 secs 
when the "RUN/STBY** switch Is on RUN; 

- data Is output to the TEAC tape drive every 0.5 secs. 


Program listing: 

;DATA ACQUISITION 


RNCNT-0 

*-80400 

OEX 

BLKCNT-2 

LDA #$92 

BNE CLOSEl 

BUFCNT-4 

STA MDRO 

JMP START 

IBUF-5 

LOA n 

RECORD 

OBUF-7 

STA RNCNT 

LDA #0 

CNT-9 

LDA #0 

STA IBUF 

BUFl-$200 

STA RNCNT+l 

STA OBUF 

BUF2-$300 

LDA #$FF 

LDA #>BUF1 

K0DRA2-$A481 

STA KDDRA2 

STA IBUF+1 

KODRB2-$A483 

LDA #0 

LDA #>BUF2 

KDRA2-$A480 

STA KDDRB2 

STA OBUF+1 

KSRB2-$A482 

STA KDRA2 

LDA #>INT 

DBR-S9008 

LDA |)$C0 

STA SA405 

WDC-$9009 

STA UACR 

LDA #<INT 

CDR-S900A 

START 

STA SA404 

MDR0-$900B 

LDA #$12 

LDA #$C0 

CSR-$900C 

STA MDRO 

STA UIER 

ESR-$900D 

LDA #P,EW 

LDA #<$C34E 

ISR-S900E 

JSR COMD 

STA UTIL 

MDR1-$900F 

LOA #REU 

LDA #>C34E 

WRT-$C1 

JSR COMD 

STA UTICH 

vmi-$c 2 

MAIN 

cl; 

ERA-SC3 

JSR OKEY 

LDA 10 

SLE-SC9 

CM? 

SXA BLKCMT 

REU-SCA 

BN*? V,.!?; 

STA BLKCNT4-1 

MRDY='$10 

LDA 

RECl 

FFT-$04 

JSR 

JSR SWAP 

DA-$20 

LDA #SLE 

JSH WRITE 

DBRE«$40 

JSR COMD 

REC2 

CCE-$80 

MAIN2 

JSR GKEY 

UDKB-$AOOO 

JSR GKEY 

CMP (fRECK 

UACR-$A0GB 

CMP #RECK 

BNE RECX 

UIER-SAOOE 

BEQ RECORD 

LOA BUFCNT 

UT1L-$A004 

CMP #CLOSEK 

CMP #150 

UT1CH-$A005 

BNE MAIN2 

BNE REC2 

UT2L-SA008 

CLOSE 

INC BLKCNT 

UT'H-$A009 

LDA #WTM 

BNE RECl 

UIFR-$A00D 

JSR COMD 

INC BLKCNT+1 

BIT5-$20 

LDA #WTM 

JMP RECl 

LOADK-$EF 

JSR COMD 

RECX 

RECK-$8F 

LDX #12 

LDA BUFCNT 

closek-$df 

CLOSF.1 

CMP #150 

TIMElH-$27 

1 LDA #ERA 

BNE RECX 

TIMEIL-SIO 

1 JSR COMD 

LDA #$40 
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STA UIER 
JSR SWAP 
U)A #SFF 
STA BLKCNT 
STA BLKCtm-1 
JSR WRITE 
INC RNCNT 
BNE RECX2 
INC RNCNT+1 
RECX2 
JMP MAIN2 


GREY 

LDA K0RB2 
PtlA 

LDA #TIME1L 
STA UT2L 
LOA ITIMElH 
STA UT2H 
CKEYl 
LDA UIFR 
AND IBITS 
BEQ GKEYl 
PU 

CMP KDRB2 
BNE GREY 
RTS 


COMD 

PHA 

LOA ESR 
COMDl 
LDA CSR 
AND #NRDY 
BNE COMDl 
LOA CSR 
AND //FPT 
BNE COMDl 
PLA 

STA CDR 
COMD2 
LDA I SR 
AND /ICCE 
BEQ COMD2 
RTS 


WAIT 

JSR WAITX 

WAITX 

RTS 


WRITE 
LDA ESR 
LOA #1M 
STA KDC 
LOA IWRT 
STA CDR 
LOA RNCNT 
JSR WWORD 
LOA RNCNT+1 
JSR WWDRD 
LOA BLKCNT 
JSR WWORO 
LOA BLKCNT+1 
JSR WWORO 
LOY #0 
WRITEl 
LDA (OBUF).Y 
JSR WWORD 
INY 

CPY #150 
BNE WRITEl 
WRITE2 
JMP C0MD2 


SWAP 

LDA OBUF+l 
PHA 

LDA IBUF+1 
STA OBUF+1 
PLA 

STA IBUF+1 
LDA #0 
STA BUFCNT 
RTS 


A. 9.2) jOATA READBACK 

Description: This program is used to read Che 

AIM 65 system’s tape drive to be sent to the ground 
The -urogram is similar to the one of Appendix A. 2, 
being 

“ no error checking is done by the AIM 65; 

- all 12 BITS of recorded data ere transferred 
have been recorded. 


WWORD 
PHA 
WUOROl 
LDA ISR 
AMO #0BRE 
BEQ WUORDl 
FLA 

STA DBR 
RTS 


INT 

FHA 

LDA UDRB 
BPL IN' EX 
TYA 
PHA 

LDY BUFCNT 
LOA #15 
STA CNT 
LDA $8000 
JQR WAIT 

roop 

LOA $8002 
STA (IBUF),Y 
LOA $8001 
INY 

LOA $8003 
STA (IBUF),Y 
INY 

DEC CNT 
BNE I LOOP 
STY BUFCNT 
PLA 
TAY 
INTEX 
LDA UTIL 
PLA 
RTI 
END 


data off the 
-based system, 
the dlfferer.ces 


just as they 
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Program listing 


;DATA RECOVERY 

RNCNT-0 

BLKCNT-2 

VRUN-4 

VBLK-3 

CNT-6 

BUFI-$200 

0BR*$9008 

W)C-$9009 

CDR-S900A 

MOR0«$900B 

CSR>$900C 

CSR-$900D 

ISR-$900E 

MDR1-$900F 

SI.P-SC8 

RDL-$C4 

REW-BCA 

NRDY-$10 

TORE- $02 

«R"$0D 

SCR-$9006 

SDR“$9007 

LOADC-$4C 

READC-$52 

CLOSEC-$43 

INALL-$E993 

NUMA-$EA46 

READM-$E93C 

0UTALL-$E9BC 

OUTPUT- $E97 A 

UTIL-$A004 

UTICH»$A005 

UACR-$A00B 


*-$300 

CCE 

.BYTE $80 
DA 

.BYTE $20 
MO 

.BYTE CR.'TAPE ERROR *,$A0 
MRUN 

.BYTE CR, 'WHICH RUN NUMBER', $BF 
MBLK 

.BYTE CR.'HOW MANY BLOCKS ',$BF 
MEND 

.BYTE CR.'LAST BLOCK THIS RU' ,$CE 
MINV 

.BYTE OR, 'INVALID COMMAN',$C4 
MERRl 

.BYTE CR, 'FILE MARK FOUN',$C4 
MRNCNT 

.BYTE CR, 'RUN NUMBER',$A0 


STA UTIL 
LDA #0 


STA UTICH 
LDA #$11 
EOR #$FF 
STA SCR 

MAIN 

JSR GCOM 
CMP Ir'LOADC 
BEQ MAIN2 
JSR INVAL 
JMP MAIN 
MAIN2 
LDA #$12 
STA MDRO 
LDA #REW 
JSR COMDA 
LDA #SLP 
JSR COMDA 
MAIN3 
JSR GCOM 
CMP #READC 
BEQ READ 
CMP CLOSEC 
BEQ CLOSE 
JSR INVAL 
JMP MA1N3 


CLOSE 
LDA #REU 
JSR COMDA 
LDA #REW 
JSR COMDA 
JMP MAIN 

RE<\0 

LDY #MRUN-MO 
JSR MESS 
JSR GCNT 
CMP #0 
BEQ CLOSE 
STA VRUN 
REAOl 
JSR RBLK 
BCS CLOSE 
LDA BLKCNT 
ORA BLKCNT+1 
BNE READ3 
LDY #MRNCNT-MO 
JSR MESS 


*-$400 
RESETB 
LDA //$92 
STA MDRO 
LDA #$C0 
STA UACR 
LDA #$68 

:$68-300BAUD 


-,$34-600 

;$1A-1200 

;$0D-2400 


LDA RNCNT 
JSR NUMA 
READ3 
LDA RNCNT 
CMP VRUN 
BCC READI 
BEQ READ2 
LDA #REW 
JSR COMDA 
LDA //SLP 
JSR COMDA 


JMP READI 
READ2 

LDY #MBL<-MO 
JSR MESS 
JSR GCNT 
CMP #0 
F£Q CLOSE 
SiA VBLK 


SENDS 

JMP 

SENDBl 

LDA 

RNCNT 

JSR 

SEND 

LDA 

RNCNT+l 

JSR 

SEND 

LDA 

BLKCNT 

JSR 

SEND 

LDA 

BLKCNT-H 

JSR 

SEND 

SENDBl 

LDY 

#0 

SENDB2 

LDA 

BUFl ,Y 

JSR 

SEND 

INY 


CPY 

#160 

BNE 

SENDB2 

LDA 

BLKCNT 

CMP 

#$FF 

BNE 

SENDB,^ 

CMP 

BLKCNT+; 

BEQ 

END 

SENDB3 

DEC 

VBLK 

BNE 

SENDBS 

SENDB4 

LDY 

#MBLK-MO 

JSR 

MESS 

JSR 

GCNT 

CMP 

#0 

BEQ 

CLOSEl 

STA 

VBLK 

SENDB5 

JSR 

RBLK 

BCC 

SENDB 

CLOSEl 

JMP 

CLOSE 

END 


LDY 

#MEND-MO 

JSR 

MESS 

JMP 

MAIN3 


; CHANGE TO NOPS 
;T0 TRANSMIT COUNTS 
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MESS 

LDA MO,Y 
PHA 

AND #$7F 
JSR OUTPUT 
IMY 
PU 

BPL MESS 
RTS 


CCOM 

JSR REACH 
JSR OUTAU 
RTS 

INVAL 

LDY #MINV-MO 
JSR MESS 
RTS 


GCHT 
LOA 10 
STA CUT 
GCMTl 
JSR INALL 
JSR DPACK 
see GCNTl 
LOA CST 
RTS 


DPACK 
CMP #'0' 
see RSPAC 
CMP <»S3A 
BCS RSPAC 
AND #$0F 
PHA 

LOA CNT 
ASL A 
ASL A 
CLC 

ADC CNT 
ASL A 
STA CNT 
PU 
CLC 

ADC CNT 

STA CNT 

CLC 

RTS 

RSPAC 

SEC 

RTS 


SEND 
PHA 
SENDl 
LDA SCR 
AND MTDRE 
BNE SENDl 
PU 

EOR l$FF 
STA SDR 
RTS 

RBLK 

LDA 

CSR 

AND 

#NHDY 

BKE 

RBLK 

LOA 

1164 

STA 

UDC 

LOA 

IRDL 

STA 

CDR 

JSR 

RWORD 

BCS 

RBLK2 

STA 

RNCNT 

JSR 

RWORD 

BCS 

RBLK2 

STA 

RNCNT+1 

JSR 

RWORD 

BCS 

RBLK2 

STA 

BLKCNT 

JSR 

RWORD 

BCS 

RBLK2 

STA 

BLKCNT+1 

LDY 

#0 

RBLKl 

JSR 

RWORD 

BCS 

RBLK2 

STA 

BUFl.Y 

IN': 


CPY 

#160 

BNE 

RBLKl 

JMP 

C0MDA2 

RBLK2 

JMP 

C0MDA4 


COMDA 

PHA 

LOA ESR 
COMDAJ. 

LOA CSR 
AND PNROY 
BNE COMDAl 
PU 

STA COR 
C0M0A2 
LOA ISR 
AND CCE 
BEQ C0MDA2 
C0MDA4 
LOA CSR 
PHA 

AND M2 

BEQ C0MDA5 

LDY MMERRl-MO 

JSR MESS 

PU 

SEC 

RTS 

C0MDA5 

PLA 

AND #$81 
BNE C0MDA3 
CLC 
RTS 

C0MDA3 
LDY #MO-MO 
PHA 

JSR MESS 
PU 

JSR NUMA 
LDA ESR 
JSR NUMA 

CLC 

RTS 


RWORD 
LDA ISR 
BIT CCE 
BNE RW0RD2 
BIT DA 
BEQ RWORD 
LOA DBR 

ac 

RTS 

RW0RD2 

SEC 

RTS 

END 
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A. 9. 3) DATA RECEIVE 

Description: This program Is written on the Hewlett Packard 

9825 of Cessna Aircraft Company. The program receives data from 
its RS232 port. The program is the same as the one of Appendix A. 3 
of Reference 2. 


Program listing: 


0 1 
1 1 

3i 

4t 

Si 

Ai 

7l 

8i 

9i 

to 

11 

12 

13 

M 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 


*tOMP. TERHINAL WITH CONTROL KEYS AND AUTO FUNCTIONS trIiOjfiUS' 

fnt l,eltz;l'»I 

tfO 2 

11-»Q 

wtc Q«1 

wtb Q,37 

0-HfHK 

wtc Q,0 

din L«(106] ,Ct0t255] ,D«1300 ,30] ,A«(80 ] ,q«C 1 00 1 
if H»lj*"-tL0 
qmb 


onl 

■Ir 

buP 

tfp 

pn 


8, “in' 

•in",L*,l 

Q,"in“,l 

"KEY" 


“IS'lif flq7)tnk 
if flq6)for L"1 


Ojefq 7j«fq 
to 203}for N> 


to 70)nvn(0«(L,N,N])4p|wtb 8,P 


if flq6>chPr(P)-«A«[Z«t4Z,Zl 
if flq6(dtp AAtnaxi 1 , l»n< AA)-31 > ,n«x<32, l*n ( AS) > ] jnvxt NsO-tZjwtb Q,13 
if flg6;n»xt L)Cfq 6{Sfq 7 
qto 16 

“KEY"tk*y-*C}if C»0|kP«t 
if C»66 or C"194j«fq 9j0-*Cjkr»t 

if flq? and C(C]>64 «nd C(C}<90|wtb Q,CtCl-64 ;cfq 9;kr«t 
if flg9jd»p "NOT A UALID CONTROL KEY" ,CIC1 jcfq 9>0'*C>kr»t 
if Z-79}13-»C 

if CICl»1008)4fq 7)0-»Cjkr«t 

if CtC]»1001 )Sfq 2>0-*C>0-»J}i-iK jkrwt 

if C(C]»1 000 |cfq 2!0-*C)kp^t 

wtb Q,C[C]->Psch«p(P)-»AS(Z>l-»Z,Zl 

dsp ASCnax < 1 ,l«n<AS>-*31 ) ,n«x <32, lwn< AS) ) ] 

if C[Cl«13j" "-tASjO-tZjdsp AS 

krwt 

"in"i 

if flg2}LS(l,lI-*QSCK,K];K>l-tK}diip <Jt>l)/5 

if flg2 «nd K-»31 ilSt 1 , 30 ]«»OS( J+! ] > J+1 -tJ { 1 -iK ; " "4QS 

buf "in" 

if C"7 or Cal35}qub "brsok" 

•ir Q,4 

tfr Q,"in*',l 

ir*t 

"string"*, 
for I»1 to 58 

i-»cm 

n»xt I 
for I"7B to 
I-30-*CtIl 
n«xt I 
for I ■98 to 

32-»cm 
nvxt I 
for 1«97 to 

i-»cm 

n«x t I 

for I-123 to 


87 


96 


122 


175 


175 


u»r‘ t i 

97 1 nait X 

sm X-17* t« 18S 
99i I>l44-»Ctn 
41 i n«it I 

All f«r I>18A tt 224 
42 1 3S4Cm 
43 t next I 

44 t 1-22S t« 2S0 
49i nwn<ch«p< I^lAI) > ^Cdl 
44 1 n«xt X 

47t 13-»CII4U(10'*CCI38] 

68i #«P X-214 tt 214 
49 1 X*lS84Cm 
7C I n«xt I 

71 1 40tCI 1721 )123-*Ct 1231 )94-*C(941 t32-tC(71 

72i 17-*CC47MCI 17Sl(40-»Ct4nHC( 1 48) )41 tC(41 ]4C( 149 1 }92-»C(2221 1 42-»C( 174) 
73 t 10140196)40(224) tl254Ctl2S)-*C (253) }434Ct 43) 4Ct 171) {454Ct4S)4Ct 1731 
74i 2740(42)40(1931)6340(43) 

72t 6440(1831)9140(1841)9340(1851 

74i 3940(1761 )04O(2O)4C( 1481)6140(61 140(189) 

77 t 5940(59)40(1871)5840(191) ) 4640(83)40(2161 )740(7140( 1351 
78 t 100140(76)40(220) )100040(73]40C2191 ) 101 040(74 140(218) 

79 t 100940(73)40(2181)100840(72)40(2171 
80 1 4440(217)40(89) 

81 1 pat 

B2i "bpath-ibuf "ln'')alp 8,4 

83i wtc q,l)wtb Q,8)w«it 200)wtb Q,37)utc Q,0 

84i wtc Q,l)wtb Q,37)wtc 8,0 

85 t pat 

86 1 and 

87i ftp I"! tt 5 

88i wpt 706,D«(1,),301,I 

89 t naxt I 

90: ftp X"1 to 200 

91; dap D9[I))w«it 50)naxt I)stp 
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A. 9.4) DATA PLOTTING 

Description: This program Is used by the Hewlett Packard 

9825 computer of Cessna Aircraft Company to convert and plot out 
the flight test results. (Sample plots are presented as Figures 
7.6.) 

Program listing; 


0i “FNG UNITS OON (quick lock dttt pit) tpkl)fila 6*;gto "START" ) 04Y4Z 
li "OON"iZ+l4Z 

2t for K«1 to 29 by 2)nun<Di( I ,K 1 )4U)nun<D«( I ,K4l ] )4U 
3 1 (K-l>/2+14H 

4 1 thf (U,-4>4U)btnd(V,15)4U}ibP<U,U)4U 

5: if bltdl ,U>)ioP<U,-4096)4U 

6t .0652B-.002427983U4H(H1 

7: naxt K)p*t 

8: "START ";7064R 

9i din F(21,0«(300,301,M(151 

10; din A( 15,501 ,A4(201 

11: ant "FIRST POSTIOM OF DYNAMIC DATA? ",p5, "LAST POSITION OF DYNAMIC", r6 

12i 10r5/24l4F(l ])10p6/24FI21)f*d 1 

13; ant "TRK * TEMP DATA ",T)if flql3)qto 44 

14; ant "FILE 9?",F 

15; dap "tapa CONTINUE" )atp 

16; tPk T)ldf F,D* 

17; “0" : 04D4r84B)Cfq 1,2 
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18 

If 

20 

21 

22 

23 

24 

25 
2fc 

27 

28 
2f 

30 

31 

32 

33 

34 

35 

36 

37 

38 
3f 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 
5? 

53 

54 

55 

da 

57 

58 

59 

60 
61 
62 

63 

64 

65 
66 1 
67 
68 1 

69 

70 I 

71 
72 1 


ant "STMT TIHF X AXI8*,r0(ant “COMTINUt AT TIMC*,r8 
if (<Fi2]-Fm>/lOar8trin>30)a«p *T00 8XC DYNAHIC* .rl 1 ; wait 70f0(Qta 11 
*Haant i aiia papar and CONTINUC” )ttp 
fnt jwrt 705,"Of*)rad 705,r2,r3,i*4,r8|r4-r2-*r2)r5-r3^r3 
if r2>lS650 ar r2<15558 ar r3>9650 ar r3< 9550 tbtpp t dap *eh«ck fl-f2*)atp 
pac 70S 

if F(21-F(n>50(F(2}<»D)m It49-*FC2] 

«ta a2 

if F(2]>D»D-tF(21 

far I-l ta 15jfar J»l ta 50$0 -»aI 1 ,31 )nait J;naat X 
fnt l,<,f6.1»fnt 2,2f4.t ,f4.0,f6.1 

fnt 3,* AL-R 8TA-L AL-L q r F da dr da*>i 

fnt 4,“ TAB-1. TA8-R At Ap Hp 8TA-R- 

"wrt Rt>.3|wrt R*.4*i 

far I>Fm to F(21 tl-FCllal >rl 

gab ‘CON* 

*caX Alpha - RH*t-.9656a22.5313M(l Hhti ,rl 1 

’cal Bata - LM’’i .4759-12. 03HC2HAI2,rll 

“cal Alpha - LH“ i-l .377-25. 7185MC3ltpC3,rll 

”cal q'‘tl0.04Sa72.0222H(4) »At4,rn 

•«al r"!! .8175-39. 968eN(5l4AC5,rll 

*cal F*i-3. 0776-57. 2165HI6l4At6,rl 1 

■cal da" il .575944. 13373MC71-. 151 lMt71-*ht7,r 11 

■cal dr"i-.4SI2S4l0.a90SMie}->A(8,rll 

“cal da ■ I -7 . 1 293247 . 35987M 1 9 1 4A 1 9 , r 11 

■sal TA8 - LH“i-2.318264n5.913M(10l46tl0,rl J 

■cal TA8 - RH"i-. 957554119. 69424aba<MMl l)4Atll ,rl 1 

if NIIIKO and nor t Igl ) . Ir 1 4r04r84B)afg 2 

■cal Ai“il0HC12l4A112,rll 

■cal Ay*ilOMI13)4A(13,rl] 

■cal Hp“t-10004S000Mtl41 AC14,rll 

■cal Beta - 3H" i -1 . 03424I2.4306MI l5l4AnS,r 1 3 

■for J"l to lljwrt R4, 1 ,ACJ,rl 1 (naxt J“i 

■wrt R4.2,Atl2,rl],At I3,pl l,All4,rll,ACl5,rl)"i 

if flg2(afg I 

dap I/10}n«xt X 

pan 

far I«1 to 15sJnp I 

»ci -.oeu ,ouM igxo *I3 

•cl -4e,30.}-250,230(qta 4ia 

acl -4,72,-480,480;qto 413 

•cl -4e,30,-2200,200(qto 412 

•cl -4,72,-3a00,1000igto 411 

•cl -4,72, -4400 , 400 iqto 410 

•cl -4,72,-20 ,460 )ijto 49 

•cl -4,72,-140,340fgto 48 

•cl -4,72,-a0,400jqto 47 

•Cl -46,30 ,-660 ,300 }gto 46 

•cl -48,30 , -520 , 440 >gt 0 45 

•cl -4,72,-30,18jgto 44 

•cl -48,30 ,-7.5, 16. Sjgto 43 

•cl -48,30,-2000,46000sqto 43 

•cl -48,30,-200,280)gto 4l 

5047 (far J<*S0 to l by -l}dap A(I,J1,J 

if AtI,J]-0)J-l4Zinaxt J 


73 1 pan# 2 

74i for J»1 to Zjplt . 1 J4rO4r04r6,Al 1 , J1 

75i if r6*Bjpen# 3fwrt 705 , "UC99, 0 ,8, 0 , -16, 0 ,8,-99" }pen jcpl t -l,0)pan6 2 
76> next Jipanjdap J 
77 1 naxt I 
76 t pan* 

79i if D>0iF[2]4l4F(lliif D>F[21 jFtl 14494F(2] >S4r84r8igto 26 

80i 04Y4Z>gto 11 

022683 
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APPENDIX B 


TRANSFORKATION OF AXES SYSTEMS 


This appendix shows the correletlon between several exes 
systems. 

Much Information contained In this section Is taken directly 
from Reference 28, which deals in depth with the problem of the 
different axes systems used in airplane analysis. 


PRECEDING PAGE BLANK NOT FILMED 
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There ere primarily five exes syetemc ueed In airplane analysle. 
These are described here. 

1) Body Axes 

"The orthogonal body-axes system Is fixed within the vehicle 
with the X-axls along the longitudinal center line of the body, the 
y-axls normal to the plane of symmetry, and the Z-axls In the plane 
of symmetry. This Is the axes system about which aircraft Instruments 
are usually mounted. Its main advantage In motion calculations Is 
that vehicle moments of Inertia about the axes are constant, so that 
the I terms can be omitted from the equations of motion. It Is the 
logical system to which to refer velocities, accelerations, and 
stability and control parameters In the study of aircraft handling 
qualities because the pilot's orientation with respect to this frame 
is fixed."* iTkli ,th& ax&4 oied -cn tkl6 Aepo/tt. ) 

2) Principal Axes 

"The principal axes are an orthogom^'.l body-fixed system for 

which the products of Inertia are z y, *he X and Z principal axes 

lie in the plane of symmetry; ' »' gle between the X body axis and 

the X principal axes is usual ^11 so that in many cases the body 

* 

axes can be assumed to coincide with the principal axes. 


A 

From Reference 28. 
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3) Flight Stability Axes 


"The flight itablllty axes (sometimes referred to as vehicle 
stability axes) are an orthogonal body-axes system fixed to the 
vehicle, the X-axes of which Is alined with the relative wind vector 
when the vehicle is in a steady-state trim condition but then rotates 
with the vehicle after a disturbance as the vehicle changes angle of 
attack. This system is preferred in many stability studies because, 
as with other body-fixed axes, the moments of inertia about the axes 
remain constant and also because the motions defined are prime-lly 

it 

those about the flight path rather than about body reference lines." 

{TkU -C4 oxe4 6yAtm uAzd ixi Re.^eAmc.t 22 . } 

4) Wind-Tunnel Stability Axes 

"The wind-tunnel stability axes are the system about which most 
wind-tunnel data are obtained. For this system the X-axis is in the 
same horizontal plane as the relative wind at all times .... The 
angle a between the X-axis of this system and the X-body axes is vari- 
able. (It is a constant for the flight stability axes.) This 
means that vehicle moments of inertia about the X-axis change. It 
also means that additional terms are required in the transformation 
equations for static-stability derivatives and for u,v,w derivatives 
when data are transferred to or from the wind axes or the wind-tunnel 

it 

stability axes." 


From Reference 
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5) Wind Ax«t 


"Th« wind uc«s ar« th« tyitan g«nnr«lly ua«d In calculating 
motions of tha vahicla as a point mass. Tha X-axis for this system 
is allnsd tha ralativs wind at all timss so that vahicla momants 

of Inartla «.bout this axis change. As with tha wind-tunnel stability 
axes, additional terms . « . are raqulrad in tha transformation to or 
from the wind axes and either the body, principal, or flight stability 
axes, since the angle . . . between the X wind axis and the X-axis o'. 
either of these systems is variable. Also, since the lateral a^'gJe . . . 
between the X-exes is variable, there are additional terms ... squired 
in the transformations for some of the lateral derivatives between the 

II* 

wind £xes and either of the other axes systems.' 

The correlation between these axes systems is perhaps best sum- 
marized by Table B.l. 


Table B.l Designation of Force and Moment Coefficients 
for Different Axes Systems* 


Component 

Coefficients for axes system - 1 

Body or 
principal 

Flight 

stability 

Wind-tunnel 

stability 

Wind 

X-axis force 

-‘=A 

S,s 



Y-axis force 

s 

S,s 


=c 

Z-axis force 

‘^z " -S 

^Z.s 



X-ixis moment (roll) 


^X.,s 

^2.,Wt 


Y-axis moment (pitch) 

^m 

^m,s 

*^m,wt 

^m,w 

Z-axis moment (yaw) 

^n 

^n,s 

C 

n,wt 

C 

n,w 


From Reference 28. 
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Transformation from the flight atabillty axea (as used In Ref- 
erence 22) to the body axes used in thla report Involves accounting 
for the eteady-state angle of attack (a^). The following equation 
takes care of this by correcting the Inertias. This is the only 
change required. 



NOTE: "s" denotes stability axes; no subscript denotes body axes. 

NASA Langley (Reference 30) and NASA Dryden (References 12-16) 
both use the body axes system. They both, however, use different 
designations. NASA Langley uses the X, Y, Z, il, m, n designation; 

NASA Dryden, the A, Y, N, i, m, n designation. The parameters will 
be presented in the X, Y, Z, 2,, m, n systum in this report. Table 
B.2 shows the correlation between both these systems. 

The symbols (i.e., Z^', etc.) in the definition column of Table 
B.2 are those as predicted by the MMLE "BONES'* program. For conversion 
from normal stability parameters (as per Reference 22) to these state 
vector derivatives, the reader is referred back to Tables 5.2 and 5.3. 

For rigorous conversion between the various axes systems, the 
reader is referred to Reference 28. 
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! 

Tabic B.2 Conparlson of Non>nim«nnlon«l Dcrlvaclvca * 


LONGITUDINAL ' LATERAL 



104 



















